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Astronomy of the University of Maryland collaborated on a graduate level 
course, “Recognition of Compact Astrophysical Objects.” This course was 
given at the University of Maryland in the Summer of 1975, and notes were 
distributed as the basic text. This publication is an edited version of those 
notes. 
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CHAPTER I 

STELLAR EVOLUTION FROM MAIN SEQUENCE STAR 
TO WHITE DWARF, NEUTRON STAR OR BLACK HOLE 

William K. Rose 

University of Maryland 
College Park , Maryland 


“ . • it is reasonable to hope that in a not too distant future we 
shall be competent to understand so simple a thing as a star. ” 

A. S. Eddington ( 1926) 


INTRODUCTION 

Even a cursory observation of the night sky reveals that the apparent bright- 
nesses and colors of individual stars are not the same. This basic observation is 
explained by the circumstance that the distances and intrinsic physical prop- 
erties of individual stars vary widely. In spite of the conspicuous richness of 
the starry sky, the notion that stars are simple objects persisted until obser- 
vations showed that their physical properties are very complex. 

Early astronomy was strongly influenced by Euclidean geometry. The de- 
velopment of geometry led to the belief that because the sphere was the most 
perfect shape, the Earth, sun and stars were also spherical. The trajectories 
of heavenly bodies were believed to represent the changeless order of the 
heavens. Until Tycho Brahe and Kepler discovered supernovae, celestial 
bodies were pictured as perfect objects that did not change or evolve. The 
discovery of supernovae represented one of the most important steps in the 
development of our present view of the heavens. 

At the present time, there is little doubt that the same laws of physics that 
explain natural phenomena on Earth can also explain the physical properties 
of stars. This circumstance leads inevitably to the conclusion that stars are 
not static but change with time. A major goal of modem astronomy is to 
explain the observations of stars in such a manner that their complex behavior 
can be understood as a natural consequence of the laws of physics. The de- 
velopment of satisfactory theories of pulsars and compact x-ray stars now 
represents perhaps the most fundamental problems in stellar astrophysics. 

The remarkable discovery of these objects has shown that stellar evolution is 
far more complex and interesting than many astronomers, such as Eddington, 
had believed. 
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Nuclear and gravitational energy are the two principal energy sources for the 
radiation that we observe from stars. For normal stars such as the sun, nuclear 
energy is the predominant source: Even if the sun were to contract and its 
radius decrease by a factor of two, the gravitational energy release would 
only be ^ 10~ 3 times the energy available by the fusion of protons to form 
helium. However, the predominance of nuclear energy sources is correct only 
for stars with normal stellar radii. For compact stars, gravitational energy re- 
lease may predominate. To compare the relative importance of gravitational 
and nuclear energy we compare the nuclear energy available from hydrogen 
6 X 10^ erg s' 1 g -1 ) with the gravitational energy of matter on the sur- 
face of the sun, of a white dwarf, and of a neutron star (see Table 1-1). It can 
readily be seen that if hydrogen is accreted onto the surface of a white dwarf, 
the energy release due to hydrogen burning would be ^ 50 times greater than 
the kinetic energy of the infalling matter. On the other hand, for the case of 
matter that is accreted onto the surface of a neutron star, the gravitational 
energy release can exceed the nuclear energy release by more than a factor 
of 10. 

Table 1-1. 


Nuclear Energy Release 
(hydrogen burning) 

Gravitation Energy Release 
GM/R 

White dwarf 

Neutron Star 

^ 6 X 10 18 erg g“* 

=^2X 10 17 erg g' 1 

^ 10 20 erg g' 1 


A black hole is a concentration of mass that is undergoing perpetual collapse. 

An observer who falls into a black hole will — in a finite time, as measured by 
liis own clock — reach a point in space that is arbitrarily close to the center 
of the black hole. The event horizon of a spherical black hole is the 
Schwarzschild radius, R s = 2GM/c 2 . An observer inside this radius cannot 
communicate with an observer outside the black hole. As measured by the 
clock of an outside observer, it will take the infalling observer an infinite 
time to reach R s because radiation emitted just outside R s is strongly red- 
shifted as it emerges from the vicinity of the black hole. The total energy 
that can be emitted by a mass m that falls into a black hole is always less than 
E = me 2 . This upper limit on the amount of energy that can be emitted from col- 
lapsing matter follows from the physical requirement that no region in which the 
total energy density is negative can exist in nature. Although in principle a 
collapsing star of mass M could emit as much energy as E = Me 2 , it is likely 
that the efficiency of energy release is such that no more than 10% of this 
available energy is ever radiated away by a collapsing star (or by matter that 
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is failing into a black hole). It is widely believed that the x-rays observed from 
most compact x-ray sources come from matter that is falling onto the surface 
of a neutron star, or falling into a black hole from an accretion disk. Matter 
falling onto the surface of a white dwarf could also produce significant x-ray 
emission. 

The discoveries of pulsars and compact x-ray sources demonstrate that very 
strong gravitational fields exist in certain regions of space. In such a physical 
environment the Newtonian theory of gravitation is no longer correct, and 
the general theory of relativity must be used to describe the physical state 
of matter. Although there is as yet no direct experimental evidence that 
general relativity correctly describes the behavior of matter in very strong 
gravitational fields, there is great confidence in the theory because it enables 
us to predict the motions of planets about the sun with very great precision. 


MAIN-SEQUENCE STARS 

If the luminosities of stars are plotted versus their surface temperatures, most 
stars are found to occupy one of three principal regions of this diagram, which 
is called the Hertzsprung-Russell diagram. These three groups of stars are 
called main-sequence stars, red giants, and white dwarfs. 

Main-sequence stars are stars of approximately uniform chemical composition 
that are burning hydrogen in their cores. The main-sequence is a mass sequence. 
Stars of large mass are luminous and hot, while stars of lower mass are rela- 
tively less luminous and have lower surface temperatures. It can be shown 
from the equations of stellar interiors that main -sequence stars satisfy a 
luminosity-mass relation which can be expressed in the form 


L = KM X 


where x varies from ^ 4.5 for IMq stars to about ^ 3 for 15 Mq stars. The 
mass-luminosity relation implies that high mass stars have much shorter life- 
times than stars of lower mass. The main-sequence lifetime of the sun is 
^ 10 10 years while that of a 5 Mq star is ^ 10® years. Moreover, low mass 
stars are formed more frequently than stars of higher mass. It follows that 
massive stars are much less common in nature than stars of lower mass. This 
circumstance is of great importance from the standpoint of the formation 
of pulsars and compact x-ray sources. The small number of such objects re- 
sults in part because they are the end products of the evolution of relatively 
massive stars 5M@). Therefore their rate of formation is much less than 
that of white dwarfs (the end product of the evolution of less massive stars). 
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A star begins to evolve off of the main sequence after approximately 1 0% of 
the mass of its core has been converted from hydrogen to helium. At this 
point the core of the star contracts and a thick hydrogen-burning shell is 
formed. The thickness of this shell source decreases as the radius of the star 
expands until the star reaches the red giant branch. 


RED GIANTS 

Red giants are stars with relatively high luminosities and large radii. Stellar 
interiors calculations show that their luminosities are primarily determined 
by the masses of their hydrogen-deficient cores. As more hydrogen-deficient 
mass is added interior to their hydrogen-burning shells their luminosities 
increase. The evolution up the red giant branch continues until the central 
temperature is high enough 10 8 K) for He 4 to burn to C 12 and O 16 by 
means of the triple alpha reactions. The luminosity of a red giant decreases 
after the ignition of core helium burning and remains relatively constant until 
helium is exhausted from the core and shell helium burning begins. At this 
point the star evolves up the red giant branch for the second time. 

Physical conditions in the cores of red giants are very similar to those in the 
interiors of very hot white dwarfs. In the cores of red giants the density is 
sufficiently high for matter to be degenerate, and the temperatures are 
sufficiently high for neutrino emission processes to play a major role in 
determining the temperature of the core. Stars lose significant amounts of 
their mass during the red giant phase of their evolution. Moreover, many 
elements are synthesized and ejected into interstellar space by red giants. 

See Figure Ma, b for most significant properties of a typical red giant. 


WHITE DWARFS 

Astronomers use the Hertzsprung - Russell diagram to classify stars according 
to their luminosities and surface temperatures. When stars are classified in 
this manner, it is found that there are numerous stars, called white dwarfs, 
with temperatures greater than or equal to that of the sun but with much 
lower luminosities. Since the luminosity and surface temperature of a star 
are related by the equation 


L = 4tt R 2 o T 4 
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Density 



Energy 

Flow 




Temperature 

(a) Red Giant Envelope 


Mass 


Figure I-la. Some important physical properties of the interior of a red giant are 
shown schematically (based on calculations by I. I. Iben). (source: Iben, 1. 1., 1967, 
Ann. Rev. Ay. Ap. 5, 571). 
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(b) Core of Red Giant 


Figure Mb. Some important physical properties of the interior of a red giant are 
shown schematically (based on calculations by I. 1. Iben). (source: Iben, 1. 1., 1967, 
Ann. Rev. Ay. Ap. 5, 571). 
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where R is the radius of the stars, T is the effective temperature and o is the 
radiation constant, it follows that these objectives must have much smaller 
radii than the sun. The radii of white dwarfs are typically one percent of the 
solar radius. 


Before the nature of white dwarfs can be understood, it is necessary to de- 
termine their mass. Fortunately, several white dwarfs that are relatively close 
to the sun are in binary systems, including the white dwarf companion of 
Sirius A. It is possible to determine their masses directly. 

The discovery that the masses of white dwarfs were comparable to the mass 
of the sun posed a dilemma to astronomers. Since the surface gravitational 
field is 


GM 



it follows that white dwarfs must have very strong gravitational fields. Since 
they are exceedingly stable, the inward force of gravity must be balanced 
by some outward pressure. What is responsible for this pressure? If the 
state of matter in the interior of a white dwarf were similar to that in the 
interior of the sun, very high temperatures would be required to produce the 
pressure gradients necessary to balance gravity. Such high interior tempera- 
tures would produce a very large outflow of radiation, and consequently 
white dwarfs would be very luminous. Since their luminosities are generally 
very low, the state of matter in their interiors cannot be that of a non- 
degenerate gas. 

It is well known that electrons (and other spin l /i particles) obey fermi 
statistics and consequently only a limited number of electrons with energies 
less than some specified value can occupy a given volume of space. It follows 
that if matter is compressed to the high densities that must exist in the 
interior of a white dwarf (p 21 10 6 g cm' 3 ) most of the electrons will have 
high kinetic energies even if the interior temperature is relatively low. 

In a white dwarf, the pressure gradients that are necessary to balance the 
strong gravitational field are produced by the zero point kinetic energy of a 
degenerate electron gas. Because the pressure in the interior depends almost 
entirely on the density and is nearly independent of density, it is possible 
to predict a relation between the mass and the radius of a white dwarf. Ob- 
servations of white dwarfs show that this predicted mass-radius relation shown 
in Figure 1-2 is nearly correct. 
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Figure 1-2. The mass-radius relation is 
shown for a white dwarf. The dashed 
line would apply if inverse beta decay 
could be neglected (Rose, W. K., 1973, 
Astrophysics , Holt, Rinehart, and Win- 
ston! New York). 


Although a degenerate electron gas can balance the gravitational field of a 
solar mass white dwarf, it is impossible to form white dwarfs with masses 
greater than approximately 1 .4 M@. This fundamental result indicates that 
a star exceeding this limit after it has exhausted its nuclear energy sources 
must collapse and become either a neutron star or a black hole. The existence 
of an upper limit to the mass of a white dwarf arises because a degenerate 
electron gas becomes relativistically degenerate at high densities (p ^ 10 6 
g cm' 3 ). The pressure-density relation for a relativistically degenerate electron 
gas is P = Kp 4 ^ 3 . If we use this equation of state and the equations of hydro- 
static equilibrium, 

dP GM(r) 

dr " P r 2 ’ 


and mass conservation, 


dM(r) 

^T" 


47T r ^ p. 


we find the dimensional relation 


KM 4 / 3 _GM 2 
R 5 “ R 5 


This relation tells us that for sufficiently large M, gravitation will dominate 
pressure for arbitrarily small values of R. 
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PLANETARY NEBULAE AND THE FINAL STATE OF THE SUN 

A planetary nebula is a mass shell that expands away from the hot central 
star at a velocity that is typically 20 km sec -1 . Most of the visual radiation 
emitted by planetary nebulae comes from emission lines. The excitation of 
this radiation is caused by ultraviolet photons radiated by the central star. 
Planetary nebulae are hydrogen-rich and have masses that range from ~ . IMq 
to ~ I.OMq. If we can infer its diameter from its measured distance and 
angular diameter, we can use the measured velocity of expansion to estimate 
the time at which a nebula was ejected from the central star. It is found that 
the central stars of planetary nebulae follow an evolutionary path (see Figure 3) 
in which their positions in the Hertzsprung-Russell diagram are determined 
by the age of the nebulae. Young nebulae have luminosities that are 

10 4 Lq while those older than cz 2 X 10 4 years have luminosities <, 10 2 Lq. 
(See Figure 1-3.) 



Figure 1-3. The evolution of a star from the red giant 
branch into the white dwarf state is shown in the 
Hertzsprung-RusseU diagram. It takes approximately 
20,000 years for a star to evolve from the position de- 
noted 0 to 20. (source: Rose, W. K., 1973, Astro - 
physics , Holt, Rinehart, and Minston: New York). 
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From their inferred lifetimes and total number in the galaxy , it can be esti- 
mated that ~ 3 planetary nebulae are produced each year in our galaxy. This 
production rate is comparable to the rate of production of white dwarfs and 
it is likely that most stars (including the sun) will pass through a planetery 
nebula phase before they become white dwarfs. 

The surface temperatures of the central stars of planetary nebulae are typically 
> 50,000 K. Such hot stars must be hydrogen deficient. It is natural to ask 
what kind of star we would produce if we constructed a stellar model from 
the hot, hydrogen-deficient central star and the hydrogen-rich mass shell. 

The only stellar model that can be constructed with these components is a 
luminous red giant. It follows that red giants must be the progenitors of 
planetary nebulae and that planetary nebulae represent a phase of significant 
mass loss that occurs in the history of a star soon before it becomes a white 
dwarf. 

The existence of planetary nebulae and other evidence for mass loss from 
red giants make it clear that stars lose very significant amounts of mass before 
they evolve into their final states. It is likely that even main sequence stars 
as massive as — 5 M 0 will lose sufficient mass as red giants to become white 
dwarfs (i.e. final mass < 1 . 4 M 0 ), rather than neutron stars or black holes. 

The theory of stellar evolution provides an explanation for the formation of 
planetary nebulae. Stellar model calculations show that when a red giant 
evolves to the evolutionary phase in which it has both a hydrogen and a 
helium burning shell surrounding a degenerate core, the helium burning 
shell is characterized by a series of thermal pulses (i.e. thermonuclear 
“runaways”). The violence of these thermonuclear runaways increase as the 
star evolves. Eventually they lead to luminosities at the base of the red giant 
envelope that are sufficient to eject much (in some cases all) of the red giant 
envelope. Observations of the planetary nebula FG Sagittae provide strong 
evidence that thermal pulses arise in evolved stars. 


PRODUCTION OF MAGNETIC FIELDS IN STELLAR INTERIORS 

Although some stellar magnetic fields are probably fossil fields (i.e. magnetic 
fields that were trapped in the star at the time it was formed), most of the 
strong magnetic fields associated with stars were probably formed by means 
of hydromagnetic dynamos that were produced during the stars’ evolution. 

As will be argued below, such magnetic fields may be produced if a region of 
a star is non-uniformly rotating and thermally convecting over several density 
scale heights. (See Figure 1-4.) 



STELLAR EVOLUTION 


11 



Figure 1-4. An idealized model consisting of a differentially rotating convecting core 
surrounded by a stationary medium with infinite electrical conductivity (i.e. rj = 0). 
The resistivity of the convecting core is 2 *15 V t (source: Levy, E. H. and Rose, 
W. K. 1974, Ap. J., 193 , 419). 


A dimensionless parameter called the dynamo number determines whether 
there are solutions to the dynamo equations, which describe the generation 
of stellar magnetic fields. The dynamo number N is defined as 



(\ V cyA 

V > 
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where rj t ~ .15 V t C (V t = turbulent convective velocity, £ = density scale 
height) is the turbulent magnetic diffusivity; R is the radius of the convective 
core or thickness of the convective envelope; y Q ~ V 0 /R is the rate of non- 
uniform rotation; and is the cyclonic velocity of the turbulently convecting 
region. If |N| becomes sufficiently large, the production of magnetic fields 
by means of a hydromagnetic dynamo is predicted. 

The maximum magnetic field strength that can be produced by a hydro- 
magnetic dynamo can be estimated once it is recognized that Coriolis force 
on convecting fluid elements is responsible for the efficient production of a 
dynamo magnetic field. For this reason, it is plausible to assume that when 
magnetic stress becomes comparable to the Coriolis force and thereby inhibits 
the motions of the fluid, the maximum magnetic field will be attained. We 
write 


<B p B t ) 
47rR 


"pV t £2, 


where B p and B t are the poloidal and toroidal components of the magnetic 
field respectively, p is the mean mass density, V t is the turbulent velocity 
of the convecting fluid and £1 is the angular velocity of the convective region. 

Where are the most likely sites for the production of stellar dynamo fields? 
Stellar magnetic fields are produced in the convective envelopes of stars 
(e.g., in the convective envelope of the sun). Convective regions produced as 
a result of thermonuclear runaways in helium burning shells and carbon 
burning cores are plausible sites for the production of those magnetic fields 
that are observed in magnetic white dwarfs (~ 1 0^ — 1 0 7 gauss) and in pulsars 
10 12 gauss). 

Main sequence stars more massive than about 1 .5 Mq have surface angular 
velocities that are typically 30 times greater than that of the sun. If the 
angular momentum per unit mass were conserved during the evolution of 
such stars, the angular velocities of their interiors would become very high 
when their interior densities had become comparable to those of white 
dwarfs (p ^ 10 6 g cm" 3 ). However, the observations of white dwarfs that 
have evolved from massive stars show that they have relatively low angular 
velocities. Stellar model calculations make is plausible to assume that 
magnetic fields in stellar interiors are responsible for the transport of angular 
momentum that is required to explain the observed low velocities of white 
dwarfs. 



STELLAR EVOLUTION 


13 


We can estimate the spin-down time scale of a rotating star in the following 
manner. Suppose that a stellar core of radius R is rotating with angular 
velocity with respect to some external region. The magnetic field exerts 
a torque <B p B t > R 3 which acts to reduce the relative angular velocity between 
the core and the external region. Since in ~ p R 5 n, we can write 


dl£2 _«cdf2 - 

P R5 -dT-<BpB t >R 3 . 


dt 


From this equation, it follows that the spin-down time for the rotating stellar 
core is 


R 2 n/<B p B t >. 

STELLAR COLLAPSE TO NEUTRON STAR OR BLACK HOLE 

Supernovae occur when the core of a star collapses to form a neutron star or 
black hole. Such collapse can occur when the core mass of a red giant exceeds 
- 1 .4 Mq the upper mass limit for a white dwarf. In our previous discussion 
we have emphasized that the relative infrequency of supernovae explosions 
as compared to the rate of production of planetary nebulae is primarily the 
result of mass loss from red giants. We have argued that the cores of red giants 
probably do not rotate as rapidly as they would if angular momentum per 
unit mass were conserved during the evolution of a star. If this were not the 
case, the differentially rotating cores of red giants could exceed 1.4 Mq by 
a significanf factor before collapse. 

A neutron star (see Figure 5) is a degenerate star whose interior density is 
comparable with or greater than the density inside the atomic nucleus (p ~ 

2.5 X 10 14 g cm" 3 ). If matter is compressed to such high densities it will 
form a degenerate (mostly) neutron gas. Neutrons, unlike electrons, strongly 
interact and therefore their equation of state is not that of a free fermion gas. 
In order to compute the equation of state of such a gas one must understand 
the forces between neutrons. If these forces are primarily repulsive, the 
equation of state will be stiffer than that of a free neutron gas. The reverse 
will be true if the forces are primarily attractive. 

The best available equations of state for neutron matter show that attractive 
interactions dominate for neutron densities < 3 X 10 14 g cm" 3 , but at that 
at somewhat higher densities repulsive interactions dominate and consequently 
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Figure 1-5. Model of a neutron star (source: Hewish, A,, 1975, Rev. Mod, Phys. 47 , 

570). 

the equation of state becomes P = K p** with 7 > 5/3. This result implies 
that neutron stars more massive than that allowable for a free neutron gas 
(i.e. ~ .7 M©) can be formed. This is an extremely important result since 
the collapsing cores of red giants are predicted to have masses that are £ 1 .4 M©. 
It follows that neutron stars, and therefore pulsars, would probably not be 
formed at all if the pressure of matter did not become greater than that of a 
free neutron gas at high density. 

Because neutron stars with relatively high mass (M > 1 .4 Mq) can probably 
exist, it follows that it is difficult to form black holes by means of stellar 
collapse. Black holes are likely to be formed during the collapse of a very 
massive star or if mass is accreted onto the surface of a neutron star that is 
a member of a close binary system. It is interesting to speculate that the black 
hole candidate Cygnus X-l might have been formed by the latter mechanism. 
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Greenbelt, Maryland 

and 
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INTRODUCTION 

The definitive history of pulsars will be a truly fascinating account of the 
visionaries, Zwicky and Landau, who postulated the existence of neutron 
stars only one or two years after Chadwick discovered the first neutral atomic 
particle, of the theoreticians, Oppenheimer and his students, who came to 
understand neutron - star structure in the years before they built the atomic 
bomb, of the astrophysicists, Woltjer, Hoyle, and Pacini, who worked on the 
nonobservable (!) consequences of neutron stars in the Fifties and early 
Sixties; of the discoverers, Hewish and his colleagues, whose persistence in 
completely exploring scintillation of radio sources, led them to detect the 
pulsating radio sources; and finally, of the many scientists who, in February 
1968 and the months thereafter, rushed to the telescopes around the globe 
or to their copies of Landau and Lifshitz Classical Theory of Fields , often 
treading on the less hardy, to explore the details of the pulsar phenomenon. 

This history has not been written. It will not be given here. There are how- 
ever many progress reports on the pulsars in the open literature, some of which 
are tabulated chronologically in Table H-l. 

In this account it is assumed that a pulsar is a collapsed remnant at the very 
end of its stellar evolution path (see the discussion by Rose in Chapter I). 
with the angular momentum and magnetic flux of a common star, 10 4143 
kg m^ s * and 10^2 24 respectively. The collapse of a common star to 
neutron star dimensions, 10 km in radius, with angular momentum and 
magnetic flux conservation, can produce spin frequencies well over 1 Hz and 
surface magnetic fields of 10^ gauss, values which lead to many consistency 
arguments concerning fundamental observations: 


‘NAS/NRC Resident Research Associate 
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Table II I 

Brief Chronology of Pulsar Investigations 


REVIEW PAPERS 

Hewish, A., et al ., 1968, Nature , 217, 709 (the discovery). 

Maran, S. P., and Cameron, A. G. W., 1968, Physics Today , 21,41. 
Hewish, A., 1969. Endeavor , 28, 55. 

Radhakrishnan. V.. 1969. Proc. ASA , 1,254. 

Ginzburg, V. I., Zheleznyakov, V. V., and Zaitsev, V. V., 1969. 

Ap. Sp. Sci , 4, 464. 

Hewish, A., 1960, Ann. Rev. Ay. Ap. , 8, 265. 

Drake, F. D., 1971 , McGraw-Hill Encyclopedia , “Pulsar”, 95-97. 
Tilson, S., 1970, IEEE Spectrum , 7, 43. 

Ostriker, J. P., 1971, Scientific American , 224, No. 1, 48. 

Ruderman, M., 1971 , Scientific American , 224, No. 2, 24. 

Smith, F. G., 1971, Rep. Prog. Phys. , 35, 399. 
ter Haar, D„ 1972, Physics Reports (Sec.C) , 3, 57. 

Ruderman, M., 1972, Ann. Rev. Ay. Ap. , 10, 427. 

Manchester, R. N., 1973, Proc. IEEE , 61, 1205. 

Hewish, A., 1975, Science , 188, 1079 (Nobel lecture), 
ter Haar, D., 1975, Contemp. Physics , 16, 243. 

Ginzburg, V. L., and Zhelezniakov, V. V., 1975, Ann. Rev. Ay. Ap. , 
13,511. 


MAJOR CONFERENCES 

1968 New York: Maran, S.P., 1968, Sky Tel ., 36, 4. 

1969 Flagstaff: Misc.. 1969, Publ, Astr, Soc. Pacific , 82, 375. 

1969 Rome: unpublished (?). 

1970 Jodrell Bank: Davies, R. D., and Smith, F. G., 1971 , IAU 
Symposium No. 46 . 

1970 Brighton: Ginzburg, V. L„ 1971, Highlights of Astronomy , 
2,20. Hewish, A., 1971 , Highlights of Astronomy, 2, 3. 
1974 Stanford: Manchester, R. N., 1974, Science , 186, 66. 
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(1) the range of spin frequencies (12); 

(2) the decrease of spin frequencies with time due to rotational energy 
loss; 

(3) the Poynting flux of the magnetic field through the “light cylinder” 
(a cylinder aligned with the rotation axis and with radius r c = c/12) 
which satisfies the energy source needed to maintain the Crab 
nebula. 

(4) an oblique dipole moment of the field which provides a natural 
source of asymmetry for the low duty-cycle pulsar emission, as 
well as a source of high-energy particles through unipolar induction. 

The following sections will cover many of the areas of observational knowledge 
and occasionally elaborate on the simple physical picture outlined above. 

They will stay closer to what we know, than to what we would like to know 
or what we presently conjecture about pulsating radio stars. 


SPIN PROPERTIES 

A narrow bandwidth sample of pulsar radiation is, to first order, a sequence 
of delta functions, not unlike the clock pulses from a good atomic standard. 

In an hour one might count 5000 pulses and know the arrival times of the 
first and last pulses referred to an atomic standard to a few ms. This leads 
to a period determination accurate to about 1/rs, or 1 part in 10 6 . Repeating 
this over several months would lead one, as it led Hewish, to discover a 
Doppler pattern in the succession of period determinations due to the 
Earth’s velocity about the sun. To investigate the nature of an extra-solar 
system periodic signal, it is necessary to reduce one’s observations to an inertial 
frame of reference. This is done by correcting topocentric pulse arrival 
time (t a ) at the telescope to barycentric arrival time (t b ) at the gravitational 
center of the solar system, a point often just outside the solar surface. The 
correction is 


t b = t a + (R.(t a ) + Ro(t a ))-S(t a ) + K(t a ), (i) 

where R e is the position of the earth center of mass with respect to the solar 
system barycenter, R 0 is the position of the observatory on - or, for terrestrial 
satellites above — the earth, S is the position of the source on the celestial 
sphere as viewed from the barycenter, and K is a correction for the error in 
the terrestrial clock as it moves in and out of the gravitational potential well 
of the solar system. Since errors in t a may be as small as 100 jus, one needs 
errors in R to be of order lOjus, or 3 km. 
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To measure t a it is necessary to average many pulses summed synchronously 
with an accurate period of date both to obtain a stable waveform and to have 
high sensitivity. With a set of bary centric arrival times distributed over a year 
or more for any pulsar (tg) one can investigate the stability of the spin by 
comparing the values tg to the best estimate for the spin frequency (£2) and 
its derivative (£2) to obtain a residual “phase” 

0i = £2(tg-tg) + j£2(tg -tg) 2 -N*, (2) 

where N* is the integer which reduces <t> 1 to a value less than 0.5. The phase 
residuals are then fit in a least -squares sense to errors in S (in Eq. (1), 12, 12, 
and perhaps 12. The standard equation of condition is 

V = R j ■ [AS 0 +J/(t‘ b - tg)] + A£2 (tg - tg) + \ A£2 (tg - tg) 2 , (3) 

where AS 0 is the vector position error at tg, (Aa, A6), and n is the apparent 
proper motion of the pulsar, (M a ,M 6 )- The effect of annual parallax is below 
the accuracy of present experiments. The final residual phase after the fitting 
procedure reflects any systematic variations of pulsar spin from the model. 
Observed values are 1 .8 <£ £2 <£ 1 80 s"' and 0" > £2/£2 2 > -4 X 10 33 . Examples 
of systematic variations are (A) discontinuities of the spin, most frequently 
observed as a sudden fractional increase of £2 by 1 0"^ to 10"® and a decrease 
of ft, and (B) slower wanderings of the residual phase by small fractions of a 
pulse period over time scales of months. These have been interpreted as 
sudden and slow changes in the moment of inertia of the star . 

Much has been written about the spin properties of pulsars. The rotational 
energy of a pulsar is 

E = il£2 2 , (4) 

where I is the moment of inertia obtained from neutron star structure calcu- 
lations, roughly 10 38 kg m 2 . The rotational energy loss rate is 

E = - 1£2£2 , (4) 

and varies between 10 22 and 10 32 J s -1 . Theories of the pulsars predict a 
spin dependence of E « £2 m with 3 < m < 5. Measurements of ft of a pulsar, 
available now only for the Crab nebula pulsar, lead to a value of 
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£ 2£2 

m= 1 + ~3.6 

£2 


( 6 ) 


A correlation of observed values of £2 and £2 (Figure 1) leads to a similar 
value for m, but the scatter is large and it is not clear that all objects belong 
to a single statistical sample. 

The ratio of £2 and £2 gives a time scale for evolution of the spin frequency: 

t e = £2/ £2 , (7) 

which is twice the age of a pulsar if m = 4 and if the magnetic dipole field 
does not decay with time. The evolution time is as small as 10 4 yr, but a 
more typical value is 10 6 yr. We view the relation of t e and age with caution 
today because of the irregular variations of £2 and £2 mentioned above and of 
the existence of an object with t e < 4 X 10 9 yr. 

Finally we mention that, in models of pulsars which spin down as a conse- 
quence of dipole radiation from an oblique magnetic field, the perpendicular 
component of the dipole field at the surface can be related to the spin 
properties by 


« (nn) 1 / 2 . 


Much labor has been spent relating one pulsar property to another within the 
context of the models we have described. Very few attempts have been made 
to account for the lack of independence of quantities such as E, t e , and B s j^ 
due to their derivation from the same observables, £2 and £2 . An important 
contribution in this regard has recently appeared (Lyne et al., 1975). They 
find no correlation between Pj(= 2 tt/£ 2) and Pj ( = 2w£2/£2 2 ), their observables, 
but they do find correlations with 408-MHz luminosity, equivalent width of 
the average pulse and galactic Z distance of 


i^os a *i u n 0 - 88 

Woe p-0.46 p-0.23 
Z a Pj -43 p+0.36 



( 9 ) 


respectively. 


(correlation - 0.43 ± 0.05), 
(-0.57 ± 0.07), 
(0.17 ±0.10), 
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Figure II-l . Spin properties of pulsars. A scatter plot of pulsar spin frequency (fl) against 
the time derivative of the spin frequency 
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AVERAGE WAVEFORM 

In fact pulsar radiation is not a sequence of delta functions. An intensity- 
time series of samples of the total power output of a narrow-band signal received 
by a radio telescope, when folded modulo the apparent period, displays a 
waveform with a typical width of 0.04 Pj . There are a number of properties 
of the total-power waveforms of the ensemble of pulsars: 

(1) Taken as an ensemble, there is no time asymmetry of the waveforms. 

(2) The waveforms indicate the presence of separate “components” 
either by resolved peaks or by inflection points which suggest un- 
resolved “components”. 

(3) In the shortest-period objects there is a high incidence for com- 
ponents at two points in the pulse period separated by 0.50 ±0.15 
periods - a main pulse and an interpulse - and in two cases for a 
low level “bridge” of radiation extending throughout the narrower 
gap between the two regions of emission. 

(4) The component-to-component peak-intensity ratios can vary with 
radio frequency to such a degree that one may effectively disappear. 

(5) The widths and separations of components frequently decrease 
monotonically with radio frequency. 

(6) The two-component objects do not appear to be the sum of two 
independent components, but rather are single double-peaked beams. 

These properties, excepting (3), are consistent with a model for the average 
pulsar beam which is, most frequently, an oblique, rotating hollow-cone 
of outer dimension 10°. An observer with a polar angle within 5° of the 
polar angle of the cone axis views a single cut through the cone — a single- 
peaked waveform when far from the axis and a double-peaked one when near 
the axis. When the cone axis and the observer are near the rotation axis, a 
broad pulse is observed with identical morphology. This scaling of morphology 
is exhibited in figures 2 and 3, which include a category of triple-peaked 
waveforms that appear to arise from a hollow cone with a central, and some- 
what unrelated, pencil beam added to the more common double-peaked 
waveform. Two beams spaced at 180° in the equatorial plane are required to 
explain interpulse components. 

Pulsar radiation is strongly polarized. A complete average waveform then 
requires four Stokes parameters: the total intensity I discussed above; the 
linear polarization parameters Q and U, which are differences of power re- 
ceived by probes of the radiation field at respective position angles of 0° 
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files with longitude expected for random orientation of observer and hollow-cone beam axes with respect to the 
rotation axis (source: Backer, 1976). 
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Figure II-3. Average waveform morphology. (See comments to Figure II-2.) 


and 90 and 45 and 135° ; and the circular polarization parameter V which 
is the power difference of probes sensitive to left and right circular polariza- 
tion It is common to display waveforms of the linear polarization P = (q2 + 
u2 ) , the position angle of the linear polarization A = 1/2 tan* 1 (U/Q) and V. 
The same quantities are sometimes given with intensive quantities normalized 
by the total intensity. 
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Average values of P/I range from 0 to very near 100% in various parts of various 
objects. The angle often rotates monotonically through the waveform. Values 
of V/I are generally less than 40%. Some polarization waveforms are given 
in Figure 4 — from observations at 430 MHz made with the Arecibo telescope. 
The polarization angle rotation is often consistent with the rotation expected 
for emission in the plane defined by the observer and cone axis vectors - 
which is the original concept that led to hollow-cone model. Examples of the 
polarization angle rotation for this model are given in Figure 5. Radiation 
from relativistic particles being accelerated along the axial field lines of an 
oblique magnetic dipole could provide the required beam geometry for both 
total intensity and the polarization angle. 


RADIO SPECTRUM 

The height of the total-intensity waveform averaged over the entire period 
is the flux-density of the equivalent continuum source. The low average flux 
densities of pulsars, a few Janskys at meter wavelengths, is the reason that 
these radio stars were not catalogued and investigated earlier. The one excep- 
tion is the Crab nebula pulsar, which was discovered and studied, especially 
by Hewish and his colleagues, in pre-pulsar years. However, at the long wave- 
lengths where the Crab pulsar is detectable, it does not pulse regularly because 
of a smearing of the pulse by multipath propagation effects (see the section 
on Interstellar Medium). The flux density spectra have steep, power-law de- 
creases at cm wavelengths, with indices between -1 and -3 and a range of m- 
wavelength morphologies from continued power laws to a second, flatter power 
law, to sharp cutoffs. Figure 6 displays several examples. The radio lumi- 
nosities of pulsars, calculated with distances estimated from the thermal plasma 
column densities (see the section on Interstellar Medium) and based on the 
conical beam model, range from L r = 10 18 to 10 23 W, or 10 16 * 21 Joules per 
pulse. The ratio L r /E is always less than 10 -3 . By comparison, laboratory 
lasers have only recently attained the megajoule level for pulses of nanosecond 
duration: 10 15 W. 

At a single frequency, radiation brightness temperatures inferred from the 
peak flux-density in the average waveform and the 10° beamwidth are near 
20 25 K. This emission, with hi; » 10* 26 J and with kT b » 10 2 J on average 
(!), constitutes the greatest thermodynamic imbalance known to man. For 
contrast, the radiation from compact extragalactic sources has hv ~ 10" 23 J 
with kT b < 10“® J. The radiation mechanism must involve the coherent 
action of a large number of charges with no perturbing influences that would 
open up channels with hv closer to kT b . Coherent radiation processes are in 









RECOGNITION OF COMPACT ASTROPHYSICAL OBJECTS 



ion axis plane - for a number ot observer polar angles. Cl. 1 
(Source: Backer, Rankin ana Campbell, 197b). 
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Figure II-6(b). Flux-density spectra. Observations of two pulsars at different epochs 
summarized indicating the existing of stable spectra (source: Backer and Fisher, 


axe 
1976 ). 
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principle, narrow bandwidth since many particles must transfer their free 
energy to radiation in unison. Perhaps as much as a decade bandwidth would 
be possible for processes like coherent magnetic bremsstrahlung radiation. 
However, the spectra of pulsars extend over two decades at comparable levels; 
in particular, the power-law form of the spectra suggests a distribution of 
parameters governing the center frequency of the coherent process in the 
emission which extends the spectrum over the observed bandwidth. 

FLUCTUATION WAVEFORM 

The individual pulses are seldom similar to the average waveform, either in 
total intensity or in polarization state. We deal at first with the total intensity 

fluctuations, and approach the matter as through a microscope, with increasing 
time resolution. 

As we have described in the preceding section, repeated measurements of the 
flux density have shown 2 to 1 variations at meter wavelengths with a charac- 
teristic time 


W - 40 da ys. (10) 

Present evidence suggests that this modulation is more severe at longer wave- 
lengths, which could imply a propagation effect or could relate to the higher 
degree of coherence required for the stronger radiation at long wavelengths 
(the degree of coherence should scale with the radiation brightness tempera- 
ture which increases with A 3 ' 5 for pulsars). At present no source of memory 
for the slow variations has been identified. 

We now focus on variations over roughly 100-1000 pulses, an increase of 
resolution by several orders of magnitude, but we are still concerned with the 
average flux density, or perhaps the average waveform. Above a few hundred 
MHz we have also hurdled over a domain of intensity variations imposed by 
variable refraction in the intervening medium. (See the section on Inter- 
stellar Medium.) In the 100 - 1000 Pj domain, many objects have a “bursty” 
intensity -time sequence, and a few exhibit sudden changes in the average 
waveform between one of two “modes”. In relation to the latter effects, 
there are also abrupt changes in emission patterns within sequences of in- 
dividual pulses, again between one of several modes. Since these variations 
are much longer than the pulse-to-pulse changes, we can label them according 
to their frequency character: 


100 Pj <t red < 1000 Pj. 
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A sequence of pulse energies unaffected by propagation effects has a fluctua- 
tion spectrum composed of several elements: 

(a) the “red” fluctuations mentioned above, 

(c) a uniform level of “white” fluctuations, and 
(b,d) one or more features indicating a quasi-periodic process in the data. 
(See Figure 7.) 



Figure II-7. Fluctuation waveform. Schematic of the typical fluctuation spectrum of 
pulsar data from a single longitude point. The components are (a) a low-frequency rise 
due to interstellar scintillation; (c) a “white” noise level arising from random subpulse/ 
micropulse variations; and (b,d) features indicating the presence of a quasi periodic 
process of a quasi periodic process in the data (source. Backer, 1975). 

The presence of a quasi periodic process leads us to increase our resolution 
another several orders of magnitude, from lPjtolO^Pj. We find that 
individual pulses are often narrower than the average waveform and very 
erratic from one pulse to the next. To the extent that there are well- 
defined peaks in individual pulses, we can identify “subpulses which many 
identify as a cut through the fundamental instantaneous beam of the emission 
from one or a few coherent groups of particles: 

tsubpulse- 10 ' 2 P 1 (12) 

The random statistics of the subpulse emission is responsible for the “white” 
fluctuations introduced above. 
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In most objects there is a definite correlation of the subpulse emission from 
one pulse to the next. This may take the form of short bursts of pulses; of a 
loosely defined pattern of subpulses on and off in different regions of the pulse 
that is replicated with a pattern period which one detects as a quasi periodicity 
in fluctuation spectra; or of a highly organized pattern of regularly spaced 
subpulses which, in the simplest case, appear to drift through the average 
waveform with a drift rate D, a separation within the pulse P 2 , and a pattern 
repetition period P 3 = D P 2 (see Figure II-8). A high degree of random be- 
havior is expected in pulsars, since the level of radiation is a very sensitive 
function of the number of particles and of their arrangement along the line 
of sight. It is less easy to understand the organized behavior. There is a need 
for a short-term memory in the emission mechanism which allows it to 
generate for several periods at a time in the most random objects, and to create 
a periodic array of beams which slip slowly out of view over tens of periods 
in the organized objects. 

One attractive explanation for these effects is the creation of an array of 
subpulse beams periodically distributed in a circle around the axis of the 
hollow-cone average beam discussed earlier in the Average Waveform section 
(Figure II-9). If these beams rotate slowly around the axis because of electro- 
magnetic body forces on the current producing the radiation, then many of the 
correlated subpulse phenomena outlined above could be modeled; e.g., sub- 
pulses which drift in a nonlinear manner, subpulses in double-component 
objects which disappear in the mid-pulse region, and sudden changes of D 
and P 3 . 

In calling this section Fluctuation Waveform we refer not just to the waveform 
of an individual pulse which deviates from the average, but also to summaries 
of how statistical properties of the second moment of the intensity data 
vary across the pulse. We are identifying the data as I jk = I (if,-, t k ) where 
<t>j is the relative position within the average waveform, or relative longitude 
of the sub-earth point of the underlying star, and t k ~ kPj . This may be trans- 
formed to fluctuation spectra: 

p (M>*i^(M> p /Vk> (12) 

where d k j and Ij k are a Fourier transform pair and (lj k ) k is the average intensity 
at longitude <j>y An example of P k j for PSR 1237 + 25 is given in Figure 10 
which displays a remarkable symmetry about the midpulse point. A histo- 
gram analysis can give an idea of the probability-density function for in- 
tensities 


P(I k l0j) = h(I k l0j)/N, 


(13) 
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(b) INDIVIDUAL 


Figure II-9. (a) Hollow-cone beam model, (b) One two-dimensional beam pattern 
which is consistent with many of the observational details discussed in this section 
and the section on average waveforms. 
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Figure IMO. Fluctuation waveform. A sequence of fixed-longitude fluctu- 
ation spectra (right) obtained from the sequence of longitudes whose 
average intensities are plotted to the left. Five components are noted in 
Roman numerals. The data display a remarkable mirror symmetry about 
the midpulse point (source: Backer, 1973). 
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where h is the frequency of occurrence histogram of intensities at longitude 
0j and N is the number of pulses analyzed. An example of p k j is given in 
Figure 1 1 for PSR 0950 + 08. From these two summaries of intensity fluctu- 
ations, one can form simple waveforms of the rms intensity fluctuation about 
the mean: 


<7(0j) = 


O(0j) = 



-ifc 


(!k-<Ijk> k ) 2 P(I k l0j) 


-i* 


p ( f k^j) 


(14) 


(15) 


In addition one can form “fractional” rms measures from P k j by summing over 
only one of the several elements of fluctuation mentioned above; e.g. a red , 

° white> ^feature- Finally, it is possible to use the phase information in the 
transform cS(f k j 0j) near f k ~ 1/P 3 to quantize the most regular drifting- 
subpulse patterns. 

Pulsar radiation has strong elliptical polarization. The polarization state of 
individual pulses — Q(0j , t k ), U( 0 j , t k ), and V( 0 j , t k ) — is not identical to 
the average state multiplied by the variable total intensity waveform. Varia- 
tions in linear polarization A(0j, t k ) = (Q 2 (0j, t k ) + U 2 (0j, t k ))*/ 2 range from a 
minimum of about 30% about the mean to a total range of 0% to 100%. The 
polarization angles 


. t k ) = x tan ' 1 


/U( 0 j ,t k )\ 

Wi’W 


(16) 


often fluctuate 30 about the most probable values, with, however, a small 
probability of being 90° ± 20° away. There are objects, or regions of objects, 
which have a nearly random angle distribution. And the circular polarization 
may fluctuate by 30% about the mean, which is most frequently less than 30%. 
Figure 1 2 is an extension of the histogram analysis mentioned above to the 
polarization state for PSR 0525 + 21 . Moving across the displays one finds j 
0j with 360° = Pj , the total intensity <I> = < 1(0 tk)), ; the vector-averaged 
linear polarization ? L = «Q( 0 j , t k )> 2 + <U( 0 j ,t k ) k > 2 )^; the average linear 
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Figure H-12. Fluctuation waveform. Probability density "map” of the polarization 
data from PSR 0525 + 21. 


polarization P2 = <(Q 2 (0 j, t k ) + U 2 (^, t k ))*> k , an intensity map of the 
probability of linear polarization with an overlay of the vector-averaged 
quantities; the analagous quantities for the polarization angle and the circular 
polarization; and ten times the fraction of data at which was above a 
threshold suitable for including the data in the three probability maps. 

The most striking property of the individual pulse polarization data is the 
double-valued angle state at angles differing by nearly 90°. It is the relative 
probabilities of these two states which often lead to the complex polarization 
angle waveforms mentioned in the Average Waveform section. By taking 
these two states into account, we find that the monotonic angle rotation 
across the pulse is common to most objects — giving further support to the 
hollow-cone emission beam model. In addition, we find that (a) in many 
cases the occurrence of the two angle states in the data has a pattern bearing 
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a fixed relationship to the correlated subpulse phenomena mentioned above, 
(b) the transition between angle states in an individual pulse is very rapid and 
at a point of low linear polarization; and (c) the two states often occur with 
opposite senses of circular polarization. 


There is no physical explanation for this phenomenon at present. The data 
suggest a subpulse beam with positive and negative lobes of linear and circular 
polarization, with the linear polarization solely in one Stokes parameter (see 
the discussion of Stokes parameters in Average Waveform section). If this 
beam drifts as described above and shown in Figure 13 for the drifters and 
objects with periodic subpulse patterns, and has a fixed orientation for other 
objects, then we can satisfactorily model the polarization data. It is curious 
that the hypothetical beam is precisely the beam of an ultra-relativistic 
charge accelerated along a circle with the symmetric linear Stokes parameter 
attenuated. While Stokes parameters are a convenient parametric expression 
of the state of polarization of radiation for calculations of radiative transfer 
(Chandrasekhar, 1963), they do not constitute orthogonal modes of the 
radiation. However, guided by the observational requirement of a single 
linear-polarization Stokes parameter, one needs to look for radiative processes 
that emit or, very importantly, perhaps amplify, one or the other orthogonal 
components of the electric field depending on the orientation within the 
total-intensity beam. Finally, we note that the polarization transition de- 
scribed above is a rapid reflection of the polarization state through the origin 
of the Poincare sphere (whose Cartesian axes are Q, U and V). 


Now for one final increase in time resolution — by about two orders of mag- 
nitude — to 10' 5 Pj . In a number of objects one finds another regime of 
deep fluctuations, micropulses, with a typical time scale, determined from 
autocorrelation analysis, of 100 /is: 

W»u« = l» 4p i- (17) 

Statistical analyses of this form of data have found it to be consistent with 
nearly random Gaussian noise amplitude-modulated by the subpulse pattern. 
The “nearly” used above is a qualification needed because in PSR 2016 + 28 
there are often periodic micropulses with P 5 ^ 1 ms (P 4 - DPj = P 3 P 1 /P 2 
is the time for drifting subpulses to cycle around the star if that is their orbit 
and if ?2 and P 3 are constant along the orbit). The simple picture for micro- 
pulses is either individual coherent emissions or random interference between 
independent coherent emissions. There is little information at present on 
the polarization properties of micropulses. 
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INTERSTELLAR MEDIUM 

The intersteUar medium contains three components which affect pulsar 
radiation: thermal electrons, a weak magnetic field, and neutral hydrogen. 

In a few objects there is or may be a circumstellar cloud of electrons and 
field but in general such an effect is not important and it will not be dis- 
cussed further here. The thermal electrons alter the group velocity of the 
radiation, yielding a dispersed signal. The arrival time of the pulse is given by 

t a ~ 0.4 DM v'g (18) 

relative to infinite frequency where is the signal frequency in units of 
10 8 Hz and DM is thermal electron column density in units of electrons 
pc cm"3 

( PSR 

n e d£. 09) 

earth 


This signal dispersion necessitates either the use of narrow bandwidths to 
obtain 10‘ 3 Pj resolution for the studies described in previous sections, 

B ~ 0.12 Pj DM' 1 MHz, (20) 

or a multi-bandwidth “antidispersion” device. Dispersion sets a time resolu- 
tion limit on direct observations of micropulses of 


dt. 


B 0 - 1/Bo, 


( 21 ) 


Bq ~2.4 X 10" 4 I'g DM' 1 MHz 2 . 

Micropulse studies mentioned in the previous section used antidispersion 
techniques to explore the radiation down to 10fis (B 0 ~ 100 KHz) at »> 8 ~ 1 
The column densities measured by pulse dispersion have led to a revision of 
the astronomer’s understanding of the interstellar medium. 

In a dilute plasma threaded by a weak magnetic field the orthogonal propa- 
gation modes of electromagnetic radiation are circularly polarized. The phase 
velocities of the two modes differ by minute amounts and, as a result, the 
polarization angle rotates with decreasing frequency. The rate of rotation is 
called the rotation measure 
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— = RM = 0.81 
dA 2 


/•PSR 

dearth 


n e d6. 


( 22 ) 


where B M is the component of the interstellar magnetic field directed toward 
the observer in microgauss. The ensemble of pulsar rotation measures have 
been used to investigate the large scale structure of the galactic magnetic field. 


It has been known for many years that the interstellar gas is not homogeneous. 
There are arms (spiral or otherwise), isolated clouds of dust which absorb 
optical radiation at Doppler-shifted frequencies of atomic lines, and clouds 
of neutral hydrogen which emit and absorb at frequencies near 1420 MHz. 
Within a year after pulsars were discovered it was clear that the electron gas 
in the interstellar medium was perturbed by parts in 10 2 on a scale of 10 11 cm, 
ten orders of magnitude smaller than the neutral clouds known previously. 

The observation that led to this understanding was the existence of deep 
modulation of the pulsar emission over a scale of many kHz at meter 
wavelengths which persisted for a thousand periods. Recall that this modu- 
lation was neglected in the preceding section. Without any further informa- 
tion, this observation led one to think of diffraction effects from scattered 
radiation rather than properties of the coherent emission mechanism. The 
straightforward dimensional analysis of the phenomenon led to the prediction 
of strong frequency dependences of the spectral and temporal scale lengths 
of the modulation, as well as a number of associated phenomena. 

Envision a layer of thickness Z filled with minute perturbations of size a and 
electron density deviation 5n. From a random walk argument the rms 
perturbation of the phase of a wavefront will be 



At meter wavelengths 50 is many radians. This will scatter the radiation by 
the angle 


50 a (24) 

where 8^ is the typical distance over which the phase differs by a radian, 
roughly a/60, and produce a Fraunhofer diffraction pattern at distances 

z >y 50 = 8|/X. 


(25) 
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The delay between direct and scattered radiation is 

r b ocZ60 2 =ZX2^ 2 , (26) 

which can become observable at meter wavelengths. Since 6^ » 1 , 2 0 scales 
as X" 1 and hence r b « X 4 . Recall the mention in the introduction that the 
Crab nebula pulsar was known as a compact continuum source at decametric 
wavelengths, due to scattering, in the years before the pulsar discovery. Adding 
random phase data together with delays of r b will result in destructive inter- 
ference with a scale 


B s <* l/r b . (27) 

Now if we allow the configuration of pulsar, irregular medium, and observer 
to change slowly because of relative velocities, the diffraction viewed by an 
observer with a narrow bandwidth (B < B s ) receiver will change in a time 

T i*V v f (28) 

where v t is the transverse component of the velocity of the line of sight rela- 
tive to the scattering medium. As mennoned above, 2^ scales as X' 1 and 
hence t s « X' 1 . The frequency dependences of 80, r b , B s , and t s from these 
simple dimensional arguments have been observed. If the neutral electrons 
are uniformly distributed in the interstellar medium and if Sn/n is constant 
everywhere, then the scattering observables should scale with dispersion 
measure, Z« DM. However, the observed dependences of r s and B s (or r b ) 
are stronger functions of dispersion measure than are predicted. 

To a DM not DM 1 / 2 , 

(29) 

B s oc DM 2-1 * not DM 2 

It is suggested that Z is proportional to DM, but that 6n grows with Z because 
of the growing probability of encountering ever stronger perturbations of 
electron density as a line of sight extends into the interstellar medium. The 
quantity (Z/a)^ 5n is a measure of the scattering. 

At meter wavelengths experiments with observatories spaced by 10 8 " 9 cm 
have been performed to attempt a measurement of the drift velocity v { oi 
the diffraction pattern since 2^ is near 10 9 " 19 cm. These experiments have 
not been very successful owing to the short length of recordings (measured 
in units of r ) and to instabilities in the pattern indicating that the medium 
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rearranges itself within a few t s . The experiments do indicate velocity com- 
ponents v t as high as 300 km s' 1 , although a typical value is 50 km s' 1 . These 
velocities are consistent with the hypothesis that pulsars receive a momentum 
impulse at birth from an asymmetric explosion, a binary disruption, or 
radiation reaction. 

At frequencies near 1420 MHz pulsar signals are absorbed by neutral hydro- 
gen clouds along the line of sight. A comparison of absorption spectrum 
features to emission spectrum features, observed when the pulsed signal is 
zero, leads to distance estimates for a pulsar based on a velocity model for 
the galaxy. Distance estimates from neutral hydrogen measurements — and, 
for the Crab nebula pulsar, from properties of the nebula — provide a crude 
calibration of the dispersion measure-distance relation. Further precision 
in the DM-Z relation must await parallax measurements using long-baseline 
interferometry. 


By using the present cude DM-Z relation we can explore the space distribution 
of pulsars. This is given in Figures 14 and 15, scatter plots in cartesian 
galactic coordinates X and Y and X and Z. One finds a distribution in Z of 
±150 pc and a general confinement to the inner 12 kpc of the galaxy. 

CONCLUSION 

We have covered most of the territory mapped out by radio astronomers 
since 1968. We have omitted, unfortunately, the large amount of theoretical 
work done on the magnetosphere electrodynamics and internal solid state 
physics of rapidly rotating, highly magnetizeu compact stars. At the present 
level of mathematical sophistication, much of this work has been unable to 
relate to the details of pulsar radiation. The results of more speculative 
theoretical work have been included in the preceding sections. 

Where are we going? In the coming years we will have firmer knowledge of 
the space distribution and motion of pulsars. Sensitive surveys are being 
conducted to find new objects, and long-baseline-interferometry studies 
have begun to determine proper motions and parallaxes. These results as 
well as complementary studies of X-ray stars (see the discussion in following 
chapters) will tie pulsars into the overall picture of stellar evolution in the 
galaxy. While every form of statistical analysis has not been attempted on 
pulsar radiation, it is not clear that the specification of more details of sub- 
pulse phenomena is in order, with the important exception of curious polari- 
zation effects. Certainly the micropulse phenomenon is worth further probes 
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to assess its “true” character; of particular importance in this regard is an 
investigation of the polarization properties of micropulses. In at least one 
object, PSR 2016 + 28, there is an indication of a polarization signature to 
the micropulses within a subpulse — a fact which suggests that the micro- 
pulse is a beam of radiation, and not merely random amplitude modulation 
within the subpulse envelope. By dealing with the properties of individual 
beams, it seems that the observers and theoreticians may be able to make 
progress in communicating with each other. 

It is clear also that in the coming years we will have a better understanding of 
the interstellar medium through “column measure” observations of pulsars, 
combined with parallax determinations, high energy photon observations of 
the interstellar medium (Hu, UV etc), and theoretical work. With pulsar 
proper motions and parallaxes in hand, we can return to doing scattering 
diffraction pattern studies and hope to isolate medium effects from pulsar 
transverse velocity effects. 

Following is a short bibliography; further references may be found in the 
review articles in Table II-l or in the NASA/STAR volumes. Appendix A 
of this chapter is a listing of data on 147 pulsars kindly provided by Y. 
Terzian. 


SHORT TOPICAL BIBLIOGRAPHY 


Introduction 


(see review articles in Table 1) 

(also two volumes of early Nature letters and articles entitled Pulsars ) 

Spin Properties 

Boynton, P. E., Groth, E. J., IIP, Partridge, R. B., and Wilkinson, D. T., 

1969, Ap.J. 157, L197 (Crab). 

Hunt G.C.. 1971. MNRAS. 153, 119. 

Manchester, R. N., and Peters, W. L., 1972, Ap. J. , 173, 221. 

Manchester, R. N., and Taylor, J. H., 1974, Ap. J. , 191, L63. 

Papaliolios, C., Carleton, N. P., and Horowitz, P., 1970, Nature, 228, 445, 
Reichley, P. E„ Downs, G. S„ and Morris, G. A., 1970, Ap. J. (Letters), 159 , 
L35. 

Reichley, P. E., and Downs, G. S., 1971, Nature P. S., 234, 48 (Vela). 
Richards, D. W., Rankin, J. M., and Zeissig, G. A., 1974, Nature , 251, 37 . 



PULSARS 


47 


Average Waveform 

Backer, D. C., 1976, to be pub. in Ap. J. , November, 1976. 

Backer, D. C., Boriakoff, V., and Manchester, R. N., 1973, Nature P.S., 243 77 
Komesaroff, M. M., Morris, D„ and Cooke, D. J., 1970, Ap. Letters. 5. 37. ’ 
Lyne, A. G„ Smith, F. G., and Graham, D. A., 1971 , MNRAS 153 ’ 337 
Manchester, R. N., 1971, Ap. J. (Suppl.) 23.283. 

Radio Spectrum 


Backer, D. C., and Fisher, J. R., 1974, Ap. J., 189, 137. 
Backer, D. C., and Fisher, J. R., 1976, in preparation. 
McLean, A. I. 0., 1973, MNRAS . 165. 133. 

Sieber, W., 1973, Astron. Astrophys., 28, 237. 


Fluctuation Waveform 


Backer, D.C., 1973, Ap. J„ 182,245. 

Backer, D. C., Rankin, J. M., and Campbell, D. B„ 1976, Natur e. 263 202 
Hankins, T. H., 1972, Ap. J. (Letters), 177, LI 1. 

Huguenin, G. R„ Taylor, J. H„ and Helfand, D. J., 1973, ApJ„ 181, LI 39. 
Manchester, R. N., Taylor, J. H., and Huquenin, G. R., 1975, Ap. J., 196, 83. 
Rankin, J. M., Campbell, D. B., and Backer, D. C., 1974, Ap. J., 188, 609 
Taylor, J. H., and Huguenin, G. R., 1971, Ap. J. , 167, 

Taylor, J. H., Manchester, R. H., and Huguenin, G. R., 1975, Ap. J., 195, 513 

Interstellar Medium 


Backer, D. C., 1975, Astron. Astrophys. . 43. 395 CScatterinpt 
Gomez-Gonzalez, J., and Guelin, M., 1974, Astron. Astr ophys., 32, 441 (HI). 
Gordon, K. J., and Gordon, C. P., 1973, Astron. Astroph ys.727.~l~l 9 (HI). 
Gordon, K. J., and Gordon, C. P., 1975, Astron. Astrophys ., 40, 27 (HI). 
Komesaroff, M. M., Hamiton, P. A., and Abies, J. G. 1972 Aust J Phvs 
25, 759, (Scattering). ’ " 

Lang, K. R., 1971, Ap. J ., 164, 249 (Scattering). 

Lazareff, B., 1975, Astron. Astrophys., 42, 25 (HI). 

Manchester, R. N„ 1972, Ap. J„ 172, 43 (RM, DM). 

Mutel, R. L„ Broderick, J. J„ Carr, T. D„ Lynch, M„ Desch.M., Warnock 
W. W„ and Klemperer, W. K„ 1974, Ap. J„ 193, 279 (Scattering). 
Rickett, B. J., 1970, MNRAS, 150, 67 (Sc^ttohig). 

Salpeter , E. E., 1969, Nature, 221, 31 (Scattering). 

Scheuer, P. A. G., 1968, Nature, 218, 920, (Scattering). 


48 


RECOGNITION OF COMPACT ASTROPHYSICAL OBJECTS 


Slee 0 B Abies J. G., Batchelor, R. A., Krishna-Mohan, S., Venugopal, V. R., 
’ and Swarup, G., 1974. MNRAS, 167,31 (Scattering). 

Sutton, J.M., 1971, MNRAS, 155,51 (Scattering). 


Evolution and Models 


Backer, D.C., and Sramek.R. 'A., 1976, AJ., 81,430. 

Ginzburg, V. L„ and Zhelezniakov, V. V., 1970, Comments Ap. Sp^Sci., 2, 


1 67 

Huguenin, G. R., Manchester, R. N., and Taylor, J. H., 1971, Ap. L , 169, 9/. 
Hulse, R„ and Taylor, J. H., 1975. Ap. J ■ (Letters), 195, L51. 

Prentice, A. J. R., and ter Haar, D„ 1969, MNRA S, 146, 423. 

Manchester. R. N., Taylor, J. H„ and Van, Y. Y., 1974, Ap. J. (Letters), 189, 


LI 19. 



APPENDIX ll-A 

LISTING OF DATA ON 147 PULSARS 


49 


50 


RECOGNITION OF COMPACT ASTROPHYSICAL OBJECTS 




PULSARS 


51 


* • - - l- • 

• r- m « - ^ 


• • • • . • 

<'* • «< «o «* • * « _ _ 


0.600997 



52 


RECOGNITION OF COMPACT ASTROPHYSICAL OBJECTS 




PULSARS 


53 





CHAPTER in 

BINARY STARS AND COMPACT X-RAY SOURCES 

William K. Rose 

University of Maryland 
College Park, Maryland 


BINARY STARS 

More than half the stars in the sky are probably members of binary sys- 
tems. Kepler’s third law implies that the orbital period P and distance 
between stars are related by 


(^Ja^-CHMj+M. (1) 

where a is the distance between the centers of the components, and M. and 
M 2 are their masses. 1 

Visual binaries are systems whose components are usually sufficienuy distant 
from one another (i.e. 100 times the diameters of individual stars) that the 

evolution of one component is practically unaffected by the presence of the 
other. On the other hand, in close binary systems the distances between com- 
ponents are comparable to the radii of individual components and the evolu- 
tion of the components is affected by the transfer of mass between them. 

Close binary systems are also spectroscopic binaries, since the orbital motions 
of the components have velocities — 100 kms"l and consequently their motions 
can be inferred from the Doppler effect. Eclipses are observed whenever the 
orbital plane of the binary system forms a sufficiently small angle with respect 
to the line of sight. Observations show that the outer layers of the components 
are distorted because of tidal interactions. 

In many binary systems (Cygnus X-l is a very important example) the spec- 
tral lines of only one component can be observed. For such systems the 
masses of the components cannot be determined without additional infor- 
mation and all we determine is the mass function 


M? sin^ i 

f(M) = — 

(Mj+M 2 ) 2 


( 2 ) 


ss 
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where M. is the mass of the primary and i is the complement of the angle be- 
tween the line of sight and the orbital plane. For an eclipsing binary system 
the angle i is close to 90 degrees. In order to find the masses of the stars in a 
single line binary system it is necessary to have additional information such 
as knowledge that one component is a main sequence star or supergiant of 
known mass. 


Calculations of the evolution of stars in close binary systems are almost always 
based on the assumption that mass transfer occurs when the radius of one com- 
ponent exceeds some critical value, which is generally taken to be the radius 
of the Roche lobe. For this reason, it is important to derive an expression for 
this critical radius in any given binary system. 

The orbital angular velocity is 


GfMj + M 2 ) 


( 3 ) 


We define a coordinate system that is rotating with angular velocity a> and 
whose origin is at the center of mass of the system. The equation of motion 


for a small point mass m is this coordinate system is 


m d 2 r = Fj + F 2 - row X (wXr) -2mu X r (4) 

dt 2 

where F. and F, are the gravitational forces of and M 2 on m. If we de- 
fine the direction of w to be along the z axis and let the x axis lie along the 
line separating the two masses, the x and y components of the above equa- 
tions become 

••= GM 2 (x-x 2 ) ' G(M t +M 2 )x 

X [(x-x^+y 2 ] 3 / 2 [(x-x 2 ) 2 +y2]3/2 a 3 

GM iy GM 2 y GQ^M^y^ (6) 

Y [(x - Xj) 2 + y 2 ] 3 / 2 [(x-x 2 ) 2 +y 2 ] 3 / 2 a 3 


where 


Mo 


-Mi 


Xi = 


, a; x 2 =• 

1 M x +M 2 Mj+M 2 


a, 


are the positions of M^ and M 2 along the x axis. 
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In the rotating frame of reference the total energy is a constant of the 
motion. To show this we note that the centrifugal force can be derived 
from the potential 


v c ~ “ \ mco 2 (x 2 + y 2 ) ( 7 ) 

and the Coriolis force, since it is perpendicular to the velocity can do no work 
in the rotating frame. It follows that in the rotating frame the total energy 
can be written 


E = ^ m (x 2 + y 2 ) + V (8) 

with 

_ -GmMj GmM 2 Gm(Mj + M 2 )(x 2 + y 2 ) 

[(x - Xl ) 2 + y 2 ] 1 ' 2 ” [(x - x 2 ) 2 + y 2]i /2 ^3 

It is convenient to define 





V 


y 

a’ 


*i 


m 2 

Mj +M 2 ’ 


$2 ~ 


Mi 

Mi + M2 


In terms of these quantities, V can be expressed as 


V 


GmCMj +M 2 ) 
a 


' ^2 + -ti 

m - «i ) 2 + v 2 ) 112 ] 112 + 1« - 1 2 ) 2 + v 2 ) 1 12 


-\(P + r, 2 ) 


( 9 ) 


A plot of the .unction V(£ ,0) is shown in Figure HI-1. We note that there 
are three critical points (called Lagrangian points) along the £ axis at which 
the potential V(£,0) has maxima. Two additional Lagrangian points are 
located off the £ axis. It can be shown by expanding the right hand mem- 
bers of the force equations in powers of the displacement about the La- 
grangian points that motions in the vicinity of these points are unstable. 
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Equipotential surfaces close to either Mj or M 2 are almost spherical. More 
distant surfaces become more distorted until they intersect at the inner 
Lagrangian point (LI) for some critical value of the potential. This equi- 
potential surface (see Figure II1-2) is called the critical Roche potential, and 
the lobes surrounding the two stars are called the Roche lobes. At larger 
values of the potential the two stars become surrounded by a common equi- 
potential surface. 



An approximate value of the radius r of a Roche lobe can be obtained by 
setting 4/3rrr^ equal to the volume with the Roche Lobe. When this is done 
the following useful expressions for the radii of the Roche lobes are found: 


M, M 1 

.8 <—<20, — =08 +.2 log— , 

M 2 a M 2 

Ml r l _ 2 / M 1 \ 1/3 

° < M^ < - 8, "a ~34/3 \Mj +M^| 

The orbital angular momentum J of a binary system is 

GCMjM^a 

Mi + M 2 


( 10 ) 


(ID 


( 12 ) 


If the total mass and angular momentum of the system are conserved during 
mass transfer, the distance from the centers of Mj and M 2 are obtained from 
the expression 




Figure III-2. A plot of the equipotential surfaces. The critical Roche equipotential 
surface intersects the inner Lagrangian point (L j ). 


a = MJ 2 _ Constant 
G(MjM 2 ) 2 M 2 (M -Mj) 2 

As a star evolves it may expand and overflow its Roche lobe. Under these 
circumstances it is possible to use the above equations to calculate how much 
mass is transferred and then to calculate the effect of this mass transfer on 
the binary system. It is important to emphasize that it is not possible to de- 
fine a critical potential surface if the rotation of the stars is nonsynchronous 

STELLAR EVOLUTION IN CLOSE BINARY SYSTEMS 

Previously, we discussed the evolution of single stars. Since one known 
pulsar and the most interesting compact x-ray sources are members of 
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close binary systems, it is of considerable interest to understand how stars 
evolve in such systems. 

The Roche equipotential surface is the critical potential surface that deter- 
mines whether a point mass is gravitationally bound to a component of a 
binary system. If the radius of one star expands so as to fill its Roche lobe, 
mass will be transferred to the other star through its first Lagrangian point. 

There are three evolutionary phases in which a star will expand and fill its 
Roche lobe. In some binary systems the components are sufficiently close 
that mass transfer will arise before the primary (i.e., and the more massive 
component) has completed hydrogen core burning. Binaries of this general 
class have orbital periods of approximately 1-2 days. In other systems mass 
transfer will occur after hydrogen core exhaustion but before helium core 
ignition. The orbital periods of the latter binary systems will range from 
~ 2 to -200 days. In still more widely separated binaries mass transfer will 
begin only after the primary has evolved up the asymptotic giant branch 
(i.e., hydrogen and helium burning shells). These relatively widely separated 
binary systems have orbital periods of years. 

For the first class of binary system discussed above, mass transfer from mas- 
sive primary to the less massive secondary will cause the Roche lobe to shrink 
and mass transfer will continue on a thermal time scale until the mass ratios 
are reversed. Thereafter, mass transfer will continue on a nuclear time scale. 

For the second class of binary, mass transfer will continue until almost all 
the mass of the primary is transferred to the secondary. Stellar model calcu- 
lations indicate that the mass of the primary will exceed the limiting white 
dwarf mass (1 .4 M® ) after mass exchange only if the initial primary mass is 
^16 M®. Mass transfer is expected to occur very quickly for widely separate 
binary systems that exchange mass when the primary is on the asymptotic 
giant branch. For such systems we expect that significant amounts of mass 
will be lost from the system. 

It is of interest to ask which binary systems are likely to be the progenitors 
of x-ray binaries. Stellar model calculations indicate that in order to produce 
an x-ray binary containing a neutron star or a black hole, the initial primary 
mass must be £ 16 M® unless a second phase of mass exchange occurs m which 
mass from the initially less massive star is added to the white dwarf remnant of 
the initial primary. 

The masses and orbital periods of such x-ray binaries as Cyg X-l and Cen X-3 
can be understood on the basis of the above picture for binary evolution. 
However, in order to produce these x-ray binaries it is likely that con- 
siderable mass and angular momentum be lost from the binary system. It is 
important to note that mass exchange before mass loss does not reduce the 
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angular momentum per unit mass of a binary system. However, a net loss 
of angular momentum per unit mass may arise during uie evolution ot stars in 
a binary system. Suppose that the masses of the primary and secondary are 
almost the same. In this case the secondary may fill its Roche lobe as the 
primary is still transferring mass, and mass from both stars will be transferred 
through the second and third Lagrangian points. Such mass transfer will cause 
a considerable loss of angular momentum from the binary system. The mass 
transferred through the second and third Lagrangian points may flow away 
very quickly or it may form an envelope around the binary stars. In the latter 
case the binary stars will experience a drag force as they orbit inside this en- 
velope. Such a drag force may cause the members of the binary system to 
move closer together. 

NOVAE 

Novae are stars that are observed to brighten by many orders of magnitude 
in a short period of time (~1 day). They are typically 10 5 times more lumi- 
nous than the sun at maximum light, and occur quite frequently. Approxi- 
mately 25 - 50 novae occur every year in our galaxy. It is of interest to 
compare the frequency of nova outbursts with the production rate of plan- 
etary nebulae and supemovae. In our previous discussion of planetary nebulae 
we have noted that their rate of formation is approximately equal to the rate 
of formation of white dwarfs (=3 per year per galaxy). Since post-novae stars 
and presumably also pre-nova stars are hot white dwarfs (also called ultra- 
violet dwarfs) it may appear surprising that their frequency of production 
should exceed that of white dwarfs. However, as we shall discuss below, 
novae must recur many times in a single stellar system. (See Figures III-2). 

Supemovae are stellar explosions that occur much less frequently than novae 
(— one supemovae per 70 years per galaxy) but are much more spectacular 
events. The amount of energy released in a nova outburst (= 10 46 ergs) is 
approximately 10-’ times less than is observed in a supernova outburst. 

Since the binding energy of an evolved star is 10 50 ergs, it is clear that the 
amount of energy released during a nova outburst is much too small to greatly 
affect the entire star. Only the outer envelope of the pre-nova star is ejected 
into space. On the other hand, the energy released during a supernova out- 
burst is sufficient to produce a remnant star that is fundamentally different 
than the pre-supernova star. 

During a nova outburst a shell of matter that is approximately 10"^ M 0 is 
ejected into space at a velocity that is ==2000 km s -1 . The kinematics of the 
expanding shell can be inferred by the appearance of highly broadened emis- 
sion lines and blueshifted absorption lines. 
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The absorption lines are blueshifted because they arise in the outer layers of 
the matter that is moving towards the observer. Some novae occur sufficiently 
close to the sun that we can measure the change in the angular diameter of the 
expanding gas shell. By measuring the velocity of expansion of the nova mass 
shell and also how rapidly its angular diameter changes with time, we can de- 
termine the distances to these novae. Since we can measure the apparent 
brightness of a nova, a determination of distance implies that we can infer its 
luminosity. Although the visual luminosities of novae at maximum light vary 
by a considerable factor, there is an empirical relation between the rate of de- 
cline from maximum light and the absolute brightness at light maximum. The 
more luminous novae tend to decline more rapidly than the less luminous. 

This circumstance implies that one can estimate the distance to a nova by 
measuring its apparent brightness at light maximum and the rate of decline 
of its apparent brightness. Since novae are quite luminous, they can be usee 
to estimate the distances to galaxies. The first observational evidence for ex- 
ternal galaxies came from observing the light curves of novae in the Andro- 
meda galaxy (M31). 

Perhaps the most important clue concerning the origin of novae comes from 
the observational evidence that all novae arise in binary systems. The typical 
nova outburst occurs in a close binary system in which one component is a 
white dwarf and the other a star of normal stellar diameter (most likely a sub- 
giant or main sequence star). The observations indicate that mass is flowing 
from the larger star onto the surface of the white dwarf. If a nova outburst 
occurs each time the amount of accreted mass equals the typical amount of 
mass ejected (=* 10 -5 M®) then many thousands of nova outbursts can arise 
from a single stellar system before the secondary (i.e., the larger star) has 
exchanged much of its mass. 

So far we have not discussed what effect the accreted matter will have on 
the white dwarf. Calculations show that as hydrogen-rich matter is accreted 
onto the surface of a white dwarf, the outer envelope will contract and 
thereby cause the interior to heat up. For plausible mass accretion rates 
(10" 8 - 10 -10 M®/year) temperatures sufficiently high to ignite hydrogen 
burning at the base of the hydrogen-rich envelope are attained after the mass 
of the hydrogen-rich envelope exceeds approximately 10 -4 M®. (See Figure 
III-4 for temperature distributions in hot white dwarf models). 

The ignition of hydrogen burning in the hydrogen-rich envelope of a white 
dwarf will produce a thermonuclear runaway that can lead to the mass loss 
that is associated with the outburst. Several mechanisms for mass loss have 
been proposed. If the thermonuclear runaway takes place under conditions 
of relatively high density (p > 10 3 gm cm -3 ) and high CN abundance, very 
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f*™* 6 ! 111 " 4 ’ 1116 temperature tT) is shown as a function of radius (R) for 0.75 M„ and 
1.395 M® models of hot white dwarfs. The 1.395 Mq model has a very high effective 
temperature (T e ) and a very steep rise in temperature near the surface (source: Rose, 
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high rates of nuclear energy generation may result (e CN ^ 10 16 erg g -1 s 1 ). 
Such high burning rates can lead to an ejection of the envelope in a dynami- 
cal time scale. Pulsations provide another plausible mechanism for mass loss 
since the white dwarf will be pulsationally unstable shortly after the ignition 
of hydrogen burning. Surface magnetic fields of 10 6 - 10 7 gauss have been 
inferred from the observed polarization of the light from the old nova DQ 
Her. Such high magnetic fields will suppress convection in the outer layers 
of the white dwarf envelope and therefore lead to mass loss when the stellar 
luminosity exceeds the Eddington limit 

L crit ~ 1038 ergS S_1 ' (H) 


Although mass accretion onto white dwarfs does not usually lead to strong 
x-ray emission, white dwarfs can be x-ray sources. This circumstance can be 
most readily seen if we equate the thermal energy density to the gravitational 
energy density on the surface of a white dwarf. We find 


T f 

1 surface r 


(15) 


MASS ACCRETION ONTO NEUTRON STARS 

The compact x-ray sources, Her X-l and Cen X-3, are believed to be neutron 
stars that are members of close binary systems in which mass accretion is re- 
sponsible for the observed x-ray emission. These systems are extremely inter- 
esting because they provide an opportunity to determine the masses of neu- 
tron stars. 

Her X-l which is one of the most carefully observed x-ray sources, has three 
periodicities associated with it. The period of the x-ray pulsations (1 .24 s) is 
the shortest period. This is believed to be the rotation period of an obliquely 
magnetized neutron star. The 1 .7 day orbital period is the next longest period. 
Finally, the x-ray source is observed to turn on and off with a period of ap- 
proximately 35 days. 

The shapes of the x-ray pulses are observed to vary in time. Sometimes there 
is a single main pulse; at other times the main pulse has two peaks. The star 
Hz Her is the optically identified companion of Her X-l . Most of the optical 
variability of Hz Her is probably caused by x-ray heating of its exposed hemi- 
sphere. The optical light curve near minimum intensity is much narrower than 
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the x-ray eclipse. Optical emission from an accretion disk surrounding Her 
X-l may explain this interesting effect. Such an accretion disk is predicted if 
matter flows from Hz Her to Her X-l through the inner Lagrangian point. 
Observed x-ray absorption dips provide additional evidence for an accretion 
disk. 

Measurements of the orbital Doppler variations of the pulsation period allow 
one to measure the projected orbital radius (3.95 X 10 11 cm), velocity (169 
km s' 1 ), and mass function 


f(M) = 


M Hz sin3 1 
(M H z + M x )2 


.85 M®. 


(16) 


The optical light curve shows that the angle of inclination is >85°. Although 
x-ray heating complicates the optical determination of the mass of Hz Her, it 
lies in the range 1 .5 — 2.5 M®. The corresponding range in x-ray mass is .7 — 
1.7 M®. If it is assumed that Hz Her fills its Roche lobe, it is possible to relate 
the mass ratio to the radius of the Roche lobe and the diameter of the occult- 
ing region. In this case the masses of the binary components become M Hz — 
2.2 Mo and M x ^ 1 .3 M®. 

Cen X-3 is another pulsed binary x-ray source. The x-ray pulses have a 
periodicity of 4.8 seconds and contain more than 70% of the x-ray emission. 
X-ray eclipses are observed, and from the Doppler shift of the period of the 
pulsations it is possible to infer that the source is a member of a massive 
binary system whose orbital period is 2.09 days. Changes in the pulsation 
period are P/P ~ 10" 4 year' 1 but are not monotonic. The x-ray spectrum, 
which has a possible low energy cutoff at 2-4 keV, is flat up to ^25 keV. 

The x-ray spectrum is harder near the peak intensity of the pulse. 

The binary companion of Cen X-3 is a BO supergiant that is at a distance of 
approximately 10 kpc. The projected orbital radius of the binary is 1.1 X 10 12 
cm and the mass functions 15.4 M®. The measured mass function implies 
that the optical companion has a mass greater than 15.4 M®. Mass estimates 
for Cen X-3 are not as precise as for Her X-l but do indicate that the mass of 
Cen X-3 is comparable to that of Her X-l . 

Although Her X-l and Cen X-3 were formed as a result of stellar collapse, 
they are both members of close binary systems. It is natural to ask under 
what physical conditions a supernova explosion can occur in a binary system 
without causing its disruption. We consider a binary system in which the 
components move in circular orbits about their center of mass. In the center 
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of mass frame velocities of Mj and M 2 are 


v 1= m 2 

V 2 =M! 


a(Mj + M2) 


1/2 


ll/2 


a(Mj + M 2 ) 


(17) 

(18) 


respectively. The supernova outburst is assumed to be spherically symmetric. 
Let q Mj equal the mass of the stellar remnant and (1 - q) M equal the ejected 
mass. Because the time scale for the supernova outburst is short as compared 
to the orbital period, the change in the potential and kinetic energy of the 
binary system is sudden. The acceleration must remain finite and consequently 
the velocities of Mj and M 2 will remain constant at the time of the explosion. 

Immediately after the explosion the kinetic and potential energies in the 
original center of mass frame are 


, -gm,m 2 

T = — [qM, + M,] 

2a(M 1 +M 2 ) 14 1 21 

(19) 

U = q — 

(20) 


respectively. It follows that the total energy in the original center of mass 
frame is 


GMiM 2 

(21) 

The kinetic energy associated with the motion of the center of mass is 


2(q Mj + M 2 ) 

where 

p = M 2 V 2 -qM 1 V 1 (22) 

is the momentum in the original center of mass frame. It follows that the 
total energy in the new center of mass frame is 
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q 

and therefore the binary will become unbound if 


(23) 


M 1 -M 2 >2qM 1 . (24) 

In most binary systems the more massive star will undergo stellar collapse 
before its companion has left the main sequence. Under such conditions the 
above equation shows that the binary system will probably become unbound 
since the mass of the remnant is generally a small fraction of the mass of the 
primary. Moreover, asymmetries in the supernova outburst itself will usually 
act to make a binary system unbound. These physical arguments explain why 
only one known pulsar is a member of a binary system. On the other hand, 
in a close binary system where mass exchange is important, the exploding 
star can be less massive than its companion at the time of a supernova ex- 
plosion and consequently the binary system can readily remain bound. 


After a supernova outburst, the binary orbit must be eccentric. It follows 
that the highly circular orbits of Her X-l and Cen X-3 must have been cir- 
cularized if a supernova outburst occurred in these systems. Although the 
physical mechanism for this circularization is not understood, mass exchange 
may act to circularize the orbit of the binary. It is interesting that the orbit 
of the binary pulsar PSR 1913 + 16, in which mass exchange is not significant, 
is highly eccentric. 

The x-ray luminosities of Her X-l and Cen X-3 are 10 37 ergs s _1 . The 
luminosity available from accretion onto the surface of a neutron star is 


GM dM 
R dt 


(25) 


and consequently the required rate of mass accretion is dM/dt a* Kr 9 M® year. 

We define the Alfvdn surface to be the surface below which the plasma sur- 
rounding a neutron star corotates with the star. The location of the Alfven 
surface is determined by the condition that the kinetic energy density of the 
infalling matter (1/2 pv 2 ) should be equal to the energy density of the stellar 
magnetic field (B 2 /8 tt). The Alfven surface is located at r A ^ 10 8 cm, if the 
surface magnetic field of the neutron star is B = 10 12 gauss. The accreting 
matter may have sufficient angular momentum for centrifugal force to halt 
its infall before it reaches the Alfvdn surface. In this case an accretion disk 
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will be formed. Viscous forces within the disk can transport angular mo- 
mentum and thereby allow an infall of plasma. 

Infalling plasma will induce surface currents above the Alfvdn surface. Such 
surface currents will shield matter above the Alfven surface from the influ- 
ence of the stellar magnetic field. Magnetic fields are likely to be frozen 
into the infalling matter. This circumstance implies that the radial compo- 
nent of the field should increase as 1/r 2 during infall, while the perpendicu- 
lar component of the field will grow as 1/r 2 . Amplification of the infalling 
magnetic field will continue until dissipative processes cause it to decay. 
Such dissipation will prevent the magnetic energy density from exceeding 
the kinetic energy density of the infalling matter; therefore magnetic fields 
will not stop infall. 


Since the velocity of the infalling matter is likely to be supersonic with re- 
spect to the Alfven surface, a shock wave will probably form in front of 
the surface. The thickness of this shock wave will be comparable to the 
distance over which the plasma is deflected, i.e. the ion Larmor radius 


12 eB 
m c c 

where v is the free fall velocity of the plasma. We note that this distance is 
less than the collisional mean free path and consequently the shock is colli- 
sionless. 

Magnetic forces of the form J X B will inhibit the flow of plasma across mag- 
netic field lines and the plasma will flow toward the polar regions of the neu- 
tron star. The magnetosphere is likely to be subject to Rayleigh-Taylor in- 
stability of some other macroscopic instability in the vicinity of the magnetic 
poles. For these reasons, the plasma may flow to the surface near the poles. 
The anisotropic radiation patterns necessary to explain the x-ray pulses are, 
therefore, produced by the anisotropic infall of plasma. 


MASS ACCRETION ONTO BLACK HOLES 

Among the known compact x-ray sources, Cyg X-l is considered most likely 
to be a black hole. Cir X-l is another plausible candidate. The x-ray fluxes 
from both these objects show variability over very short time scales. In the 
case of Cyg X-l the fluctuations extend down to times less than 1 ms. 
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Cyg X-l is known to have at least two x-ray states (an upper and a 
lower state). In the upper state the source has an excess of low energy 
(i.e. < 6 keV) x-rays. In 1971 a radio source appeared at the position of 
yg X-l at the same time that the x-ray emission changed from the upper 
to the lower state. During this transition the 2-6 keV flux decreased by a 
factor of 4 and the 10-20 keV flux increased by a factor of 2. Recent (1975) 
x-ray observations show that Cyg X-l has returned to its high state. 

Cyg X-l is a member of a 5.6 day spectroscopic binary. Its companion is the 
9th magnitude BOI supergiant HDE 226888. Recent observations indicate 
that the distance to HDE 226888 is >2 kpc. This distance determination is 
extremely important. If HDE 226888 is more distant than 1 kpc, it must be 
a massive supergiant (i.e. M > 20 M®). A lower limit can be placed on the 
mass ratio M(Cyg X-1)/M(HDE 226888). It follows from this measured lower 
limit that Cyg X-l must have a mass <J4M®. This mass is probably too high for 
a neutron star or white dwarf, and consequently Cyg X-l is probably a black hole. 
Attempts to interrupt the x-ray emission have assumed that the x-rays are coming 
from an accretion disk that surrounds a black hole (M >4 M®). A mass accre- 
tion rate of 10~ 9 M®/yr may be sufficient to produce the required x-ray emis- 
sion. It is important to point out that spherical accretion is a very inefficient 
means of producing radiation, since only a small fraction of the available 
energy is emitted during infall. Moreover, spherical accretion does not pro- 
duce strong 1-10 keV x-ray sources. For this reason, accretion disks are the 
only known plausible models for intense 1 -1 0 keV x-ray emission from black holes. 

Two basic types of mass accretion can arise in close binary systems. For the 
case of Roche lobe overflow, the infalling matter has a great deal of angular 
momentum and consequently an extensive accretion disk will be formed. In 
the second type of mass accretion, mass transfer occurs as a result of a stellar 
wind. In this case it is touch and go whether or not the infalling matter has 
sufficient angular momentum to form an accretion disk. 


The above discussion makes it clear that before we can understand x-ray emis- 
sion from black holes we must understand the physical properties of accretion 
disks. Such disks must have outer and inner boundaries. The location of the 
outer boundary is determined by the specific angular momentum of the accret- 
ing matter. The location of the inner boundary is determined by the circum- 
stance that, by general relativity, it is impossible for stable particle orbits to 
exist closer than =* 3r s (r s = 2GM/c 2 ) to a black hole. In discussing the struc- 
tures of accretion disks we will assume that mass accretion is steady and that 
matter in an accretion disk is in Keplerian orbits (i.e., centrifugal force balances 
gravity in the radial direction). 
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We assume that dM/dt is the rate at which gas is added to an accretion 
disk. The specific angular momentum at a radius r from the center of a 
black hole is 


l ,(2 L) 1,2 ,= ( cm,)‘'2. « 

An element of mass must lose almost all of its angular momentum before it 
reaches the inner radius. For this reason, angular momentum must be con- 
tinuously transported away from the outer boundary of the accretion dis 
by infalling gas under steady state conditions. The rate of removal of angu- 
lar momentum is 




(27) 


where r Q is the radius of the outer boundary. 

The surface density of matter in the disk is 

2 = 2hp ( 28 ) 


where h is 
two sides, 
stress is 


the thickness of the disk and p is the density half-way between its 
The radial velocity is v < 0 and v^ = u> r = (GM/r) 1 . The viscous 


2l 7°0r 


(29) 


where r? is the dynamic viscosity and the shear of the Keplerian orbits is 


3 3/GMX 1 / 2 

From the conservation of mass, it follows that 


(30) 


dM _ rv 

— — = 27rrv r 2. 
dt 


(31) 


At any particular time angular momentum is transported across any radius r 
by the infall of mass. In order to conserve angular momentum outwards. 
Angular momentum conservation implies 
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— (GMr) 1 < 2 = 2m 2h^ r r + J c . (32) 

The left hand side of the above equation represents the rate at which angular 
momentum is transported inward through a radius r and 


j c ^— (GMr,) 1 / 2 (33) 

(where r, is the inner boundary of the disk) is the rate at which angular mo- 
mentum is deposited on the black hole. The remaining term represents the 
rate at which viscous stresses transport angular momentum outwards. 

The above equations imply that the viscous stress is determined by the 
mass accretion rate and the mass of the black hole. Solving for the viscous 
stress, we find 


dM 

2h V = ~ 2 [(GMr) 1 / 2 - (GMr,) 1 / 2 ]. (34) 


If at any r, 2hT,0 r is smaller than given in the above equation, the viscous 
stresses will not b^ sufficiently strong to allow the mass flow dM/dt and con- 
sequently mass accumulation will take place. Mass accumulation in the outer 
regions of the disk will lead to an increase in its radius. 


Fhe rate of generation of heat per unit area in the disk is 


2he = -4h - 


3 dM GM 

<t>i °4>i ~ 4jrr 2 dt T ’ < r>>r i)> 


(35) 


and consequently the rate of generation of heat between r, and r 2 is 



( 36 ) 


r 2 ^ r l » r ] 
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As shown below, heat is produced in the disk by viscous stresses as well as 
gravitational energy release. Gravitational energy is released between r L and 
r 2 at the rate 

dM /GM GM\ (37) 

It follows from the virial theorem that 1/2 of this energy goes into heat and 
1/2 into orbital kinetic energy ; therefore viscous forces which transport angu- 
lar momentum outwards deposit heat between rj and r 2 at the rate that is 
2/3 the total rate of generation of heat given in equation 36. 

The rate at which energy is transported across a radius r is 


3 dM GM 

E = ui ^iTt T 


( 38 ) 


The energy flux from the top and bottom of the disk is 2F, where 

3_ dM GM (3 , 

dt r ’ 

We note that F is independent of the viscosity and depends on r, the mass 
accretion rate, and the mass of the black hole. The total luminosity from 
the disk is obtained from the expression 


L = 



2F 27rrdr. 


(40) 


The minimum disk temperature T min (r) at a given radius r is achieved if the 
disk radiates as a blackbody; i.e., 


-aT 4 . =F(r). 

4 mm v ' 

On the other hand, the maximum value for T(r) is determined by 

GM m P _ m P c2 

T max- r, k' 6 


(41) 

the relation 

(42) 
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The x-ray observations of Cygnus X-l indicate that the hottest electron tem 
perature T e in the disk surrounding Cygnus X-l is T e 10 9 K, which is be- 

tween T min and T max . 

The vertical structure of an accretion disk is much less well determined than 
the radial structure, since the determination depends on a knowledge of the 
viscous stresses in the disk. It is plausible to assume that 

V r — 

If a were known, then the vertical structure could be determined by the con- 
dition of hydrostatic equilibrium 


the equation of state 



pk(Tj + T e ) 
nip 


(43) 

(44) 


some assumption about the relative magnitudes of Tj and T e , and angular 
momentum conservation. It has been suggested that Tj » T e in the inner 
portion of the accretion disk surrounding Cygnus X-l. In this picture the 
inner part of the disk bulges because Tj > 10 n K but T e is kept < 10 9 K as 
a result of Comptonization of relatively soft x-rays (<6 keV) emitted by the 
outer regions of the disk. This theory explains the spectrum of the hard 
x-rays (> 10 keV) from Cygnus X-l . 

According to general relativity the motion of a particle moving about a 
spherically symmetric black hole (i.e., moving in a Schwarzschild field) can 
be written in the form 


where 





(45) 


and 


r S “ 


2GM 


h = r 2 0 


= constant 


is the specific angular momentum and r is the proper time. The above expres- 
sion differs from the corresponding Newtonian result in that as a result of the 
last two terms on the right hand side of the equation, the effective potential 
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energy becomes large and negative for sufficiently small r. Thus, in general 
relativity gravitational capture is possible in the absence of tidal forces, and 
stable orbits about a black hole are impossible for sufficiently small r (or suf- 
ficiently high E and fixed h). 
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CHAPTER IV 

X-RAY SIGNATURES OF COMPACT OBJECTS 


Elihu Boldt 

Goddard Space Flight Center 
Greenbelt, Maryland 


INTRODUCTION 

Although the concept of an extremely compact stellar object (i.e., neutron 
star or black hole) has existed as a theoretical curiosity for a long time, the 
fundamental observational consequence of detectable X-radiation was some- 
how overlooked prior to the advent of experimental X-ray astronomy. We 
recall what a surprise it was to find that X-ray astronomy is even possible, in 
the sense that there are any detectable signals at all from beyond the solar 
system. Then there was the greater surprise that over a hundred stellar sources 
within our own galaxy each radiate X-rays at more than 10 36 ergs/s (i.e. the 
luminosity of a thousand suns, in X-radiation alone). Now, after looking at 
the temporal and spectral signatures of several of the brightest sources, we 
suspect that most of the strongest sources, if not all, are associated with com- 
pact objects. Hence, we have come to realize that the main reason X-ray 
astronomy is flourishing is that compact stellar objects do in fact exist and 
that intense X-radiation is a direct consequence of energetic processes funda- 
mental to them. 

We can now go back and do the homework we should have done about twenty 
years ago. That is, we can see how far we could have gone in discussing the 
potential X-ray characteristics of compact stellar objects without knowing 
anything other than the most basic theoretical parameters of compact objects 
and some elementary laws of physics. Of course, given the advantage of know- 
ing something about the actual X-ray data, this procedure is obviously artifi- 
cial. Nevertheless, it may help to give us the kind of background needed for 
interpreting these data in a discriminating way and extracting the genuinely 
new information they must carry. We will try this out in practice later by 
examining some specific data which from the work of the x-ray astronomy 
group at Goddard Space Flight Center; provide us with a good challenge 
for three of the most interesting objects yet discovered, viz: Cye X-l Cyc X-3 
and Her X-l . ’ 
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II. X-RAY SCALES FOR COMPACT OBJECTS 

As was indicated by Dr. Rose in Chapter I, a neutron star is expected to have 
a mass approximately that of the sun and a radius of about 10 6 cm. Hence, 
we introduce the dimensionless parameters Af and R^ defined by: 


M = M/M 0 


R 6 = R/(10 6 cm). 

The velocity (v ff ) of freely falling matter, as measured at R, is always given 
by 

v ff =(2GM/R) 1 / 2 (0 

Expressed in our dimensionless parameters, this becomes 

v ff = 0.542 c(A//R 6 ) 1 ^ 2 . (2) 

In other words, for a typical neutron star, the free fall velocity is about half 
the velocity of light. Hence, we note immediately that collisions with am- 
bient matter will give rise to bremsstrahlung radiation in the X-ray band. 
Specifically, the endpoint for the bremsstrahlung spectrum of photons ob- 
served at infinity, emerging from the collisions of matter falling to the surface 
of a neutron star, is given (with m 0 the electron rest mass) as 

(hf)« = m 0 c 2 {[l -(v/c) 2 ]* 1 / 2 - l} {l -(Vff/c) 2 } 1/2 . (3) 

For Af = R^ and v = v f f this gives 82 keV as the endpoint, placing the bulk of the 
emission in the usual hard X-ray region. On the other hand, for M > 3.4 R^ , 
v f f = c and Equation (3) shows that (ht') 00 = 0. Hence, if R^ — is funda- 
mental for a neutron star, then the mass cannot exceed 3 M® without 
the object becoming a black hole. X-radiation from near the surface of a 
compact object is then quite natural for a neutron star but impossible for a 
black hole. For a non-spinning black hole (Schwarschild metric), the orbital 
velocity for the innermost stable circular Keplerian orbit (Kepler radius 
R^ = 0.88 Af) is precisely half the velocity of light and we again have condi- 
tions natural for collisional production of X-rays, but well above the surface 
of the star (Schwarschild radius R^ = 0.29 Af). 

The other basic parameter we consider is the magnetic field (H). A field of 
10 6 - 10 7 G for a white dwarf or the core of a red giant would, with flux 
conservation, become a field on the order of 10 12 - 10 13 G when the object 
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collapses to a neutron star. Again we introduce a dimensionless parameter 

H 12 =H/(10 12 G) 

The quantum of the magnetic field (i.e., the energy of cyclotron photons) is 
then 


* H = H 12 (11.6 KeV). (4) 

This is obviously in the X-ray band and implies important emission and ab- 
sorption effects, as we shall describe later. Curiously enough, this huge mag- 
netic field disappears as stellar collapse approaches the Schwarzschild radius 
and these magnetic effects are strictly absent for black holes. 


III. ELECTROMAGNETIC RADIATION LIMITS 

The general classical expression for the rate of energy radiated by an electron 
of velocity (v) in a field (E, H) is given by 


where 



( E *7 XH ) l -i< E 



7 2 = [1 -( v / c ) 2 ]” 1 and a 0 = j ir e 2 /mc 2 ) 2 


( 5 ) 


is the well known Thomson cross-section. We consider an electromagnetic 
environment of photons to obtain the results appropriate to inverse Comp- 
ton scattering. We specify an ambient photon propagation direction by a unit 
wave vector (k) such that 


H = k X E, k *E = 0. (6) 

Equation 5 may now be written 

W r =o 0 U EH c [ 1 “0 * v/c)] 2 /[l -(v/c) 2 ] (7) 

where U EH = E 2 /87r is the energy density of such photons. In order to in- 
vestigate the “force” exerted by these ambient photons we note that the rate 
of radiated momentum is W r v/c 2 and that of absorbed momentum is W a k/c, 
where W a is the energy absorbed. Using energy and momentum conservation, 
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we find that the force in the k direction is (for v 2 « c 2 ), given by 

= °o ^EH ' v / c )l’ 00 

Averaging over photon propagation directions we get the net force arising 
from all ambient photons, as follows 


<F > = * 0 U EH [<k> - 2 k^>]. (9) 


For particles falling onto a compact object the second term is always negative; 
hence 

( F ) ^ o 0 Ugjj <k >, (10) 

and the net total force inward (gravitational minus this electromagnetic radia- 
tion force) is always 


net 


l< 


°o L 1 


tJMm 
_ R 2 4?rR 2 c 


( 11 ) 


where L is the luminosity, given by 

L = 47rU EH c<k -R>R. (12) 

This net total force vanishes when 

L = (47rGMmc)/ a 0 = L E , (13) 

where L E is so defined (i.e., the Eddington luminosity limit). 

The net total force may then be rewritten as 



That is, in effect, the gravitational source term (GM) is reduced by the factor 
(1-L/L E ). 

Because of the charge neutrality of the infalling plasma, the total mass per 
electron is essentially one atomic mass unit (amu); the protons couple to the 
gravitational field, the electrons couple to the ambient radiation field, and 
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Coulomb attraction couples the electrons and protons. Hence the mass (m) 
in Equation 13 refers to a mass of 1 amu, and the Eddington limit so eval- 
uated is 

L e = 1.3 X 10 38 M erg/s. ( 15 ) 

It is convenient, therefore, to measure luminosity in units of 10 38 ergs/s and 
define a dimensionless luminosity parameter (Lg) as 

Lg = L/10 38 ergs/s) 

If this limiting luminosity (i.e., Equation 15) arises from blackbody emission, 
the temperature is given by 

T = 2.0Ml/4(R 6 )l/2 X 10?K, (16) 

providing a spectrum that resides well within the X-ray region. Measuring 
temperature in units of 10 7 K, we define 

T 7 = T/(10 7 K) 

and obtain the general blackbody expression for temperature in terms of 
luminosity and size as 


T 7 = 1.9l4 /4 R^ 2 . (1 

For optically thin emission, the Eddington limit places restrictions on the 
population of radiating particles. For an ambient magnetic field (H=£ 0, 

E = 0), Equation 5 gives (with j3 = v/c) 




(18) 


Since the number of radiating particles is 

N = L/W r (19) 


the Eddington limit (i.e., Equation 15) and Equation 18 give 


N < 4 X 10 39 M(H 12 )-2 r 2 


( 20 ) 
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where equation 4 implies that the electron velocity (0) in Equation 20 needs 
to exceed a minimum value, given by 

0 2 > 0.045 H|2* ( 21 ) 

Hence, 

N^XIO^MH"^ ( 22 ) 

The upper limit on the kinetic energy content W e of electrons corresponding 
to Equation 19 is 

w e < 2 X 10 33 M H ‘ 2 2 ergs. (23) 

This upper limit will be particularly important when we examine the energy 
content of flare events associated with some of the compact sources. 

For bremsstrahliing radiation, we use Equation 5 to give the radiation rate 
under the condition that 


<(V-EJ2> =1e2v2. (24) 

obtaining , -v 

W r = 20 c c(^l-^Ju E , (25) 

where the electric field energy density U E is given (with a — 1/1 37) by 

U E = Z 2 m 0 c 2 0n/(1377r), (26) 

n being the number density of nuclei (charge Ze). Considering an equal num- 
ber density of electrons and protons, the Eddington limit here puts a limit on 
the emission measure of the emitting region: 

/ W r n d (vol.) < L e (27) 

which, with Equations 25 and 26, may be expressed 

f n 2 d (vol .) < 1 .7 X 1 0 60 M (\ - y j fT 1 cm" 3 (28) 

where, for the X-ray region, we require 0 =;0.05 - 0.5. Since <n 2 > is always 
greater than(n) 2 we may rewrite Equation 28 as 
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N 2 /(volume) < 1.7 X 10 60 A/ ^1 J /H cm -3 . (29) 

In terms of the energy content W of radiating particles, Equation 29 yields 
W< 7 X 1029 Afl/2 1/ % 3 / 2 (volume)l/ 2 me 2 . (30) 


This will help us isolate the particle population responsible for the burst events 
seen in Cygnus X-l (a black hole candidate). 


IV. ACCRETION 

For emission driven by an accretion process with matter initially at rest, the 
luminosity is given by 


L = 


dM\ GM , 

d r -r('-wle). 


where the rate of mass accretion (dM/dt) is given by 

^ = 4ttpR 3 / 2 t/ [2GM (1 - L/L*;)] 1 / 2 


(31) 


(32) 


with r? the ratio of infall velocity to free-fall velocity, 

V = <(V • R/R)> [2GM (1 - L/L e )/R]-1/2. (33) 

From the above equations we obtain a lower limit on the density, given by 
P > 3 .7 X 1 (H> L 8 M- 3 / 2 R^ 1 ^ 12 g/cm 3 (34) 

(i.e., obtained for i? = 1 , L « L E ). 

For a plasma at a temperature , this may be used to infer a lower limit on 
the thermal kinetic energy density as 


U T =-|kT > 29 X 10 20 T 7 M-3/2 r- 1 ' 2 eV/cm 3 . 


(35) 
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For free-fall, the kinetic energy density U ff would be given by 



7 M ' 

\l/2 

. ( M \ 


U ff = pc 2 

Lr- i94 v 

) J 

- pc \ 015 r;)- 

(36) 


From Equation 36 we note thatU ff ^ 0.2 pc 2 near the surface of a neutron 
star or the innermost stable orbit of a Schwarschild metric black hole, and 
U ff ^ pc 2 near the Schwarschild radius of a black hole. 

From Equation 34, we get 

pc 2 > 21 X 10 26 L 8 j W- 3/2 R' 1/2 eV/cm 3 . (37) 

For comparison, the energy density of the magnetic field is 

U H = 25 X 10 33 H 2 2 eV/cm 3 . (38) 

Considering a dipole field (i.e. H oc R -3 ), the minimal kinetic energy density 
(Equations 36-37) would equal that in the magnetic field at a radius approx- 
imately 10-1 0 2 times that for the maximum field strength. Hence, particle 
motion in this inner region is likely to be dominated by magnetic effects, 
leading to beaming and pulsed modulation as the dipole rotates. At a radius 
approximately (10 5 /T 7 ) that of a neutron star, we expect (from Equations 
35 and 36) that the accreting plasma is dominated by thermal and/or external 
effects. The vast intermediate region is of course dominated by the gravita- 
tional field of the compact object. 

V. RADIATIVE TRANSFER 

Having established a lower limit to the density of accreting matter as a func- 
tion of radius (see Equation 34) we are now in a position to evaluate a lower 
limit to the column density of the entire accreting atmosphere through which 
we might view the X-ray emission, viz: 


pdr >7.4 Rj /2 LgM- 3 / 2 (T 1 ! 2 - l)g/cm 2 (39) 

Jr 6 

where we have assumed that this atmosphere extends to at least TR^ . First 
we notice the remarkable Feature that this immediately implies a highly ionized 
atmosphere, for otherwise photoelectric absorption would render the source 
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invisible in the usual X-ray region (i.e., below =* 20 keV). Even if the gas were 
so hot that all elements of lower Z than iron were completely ionized, the K 
shell electrons of iron at cosmic abundance alone would cause an appreciable 
absorption dip at about 9 keV. Hence, the temperature of the gas in this 
atmosphere is likely to be greater than 10® K and is almost certainly greater 
than 10 7 K. 

Using Equation 34 for density, we evaluate the optical depth for free-free 
absorption as (t v ) as 


r„ > 6X10- 5 (T 7 )- 7 / 2 £n T (40) 

where x = hi>/kT and r is the ratio of the outer radius to the inner radius de- 
fining the atmosphere considered. For T 7 > 1 , Equation 40 implies that the 
atmosphere could be extremely thin to emission and absorption. This is a 
lower limit to the optical depth, due mainly to the fact that we took 7 ? = 1 in 
evaluating the lower limit for density given in Equation 34. However, it is not 
likely that r\ would be smaller than about 10" 3 . Although this would increase 
the optical depth by about three orders of magnitude, it would still be less 
than unity (i.e., optically thin). 

The optical depth for Thomson scattering is (r s ) 

r s = 2.9Rl /2 (rl/2-l)L 8 A/-3/2 (41) 


where is the radius of the source and r(R^) is the outer radius of the 
scattering atmosphere. It is quite likely, therefore, that this atmosphere be 
optically thick to Thomson scattering. We note, however, that Thomson 
scattering may be substantially altered in the presence of a strong magnetic 
field, since the transfer of energy to electron motion normal to the field is 
quantized at intervals Ifco = X 1 1 .6 keV and thus is severely restricted. 
The modified cross section (a H ) may be expressed in terms of the Thomson 
cross section (a H ) as follows: 


( fl ) : 


(hy ± 1 1.6 keV H 12 ) + 2 Sm2fl 


(42) 


where 6 is the angle between the field direction and the direction of photon 
propagation. For hv « 1 1 .6 keV H J2 , Equation 42 shows that the cross 
section at small angles can become much smaller than the value given by the 
unperturbed Thomson cross section. This opens up the interesting possibility 
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of an optically thin direction along the dipole axis even though the atmos- 
phere in general might be optically thick to electron scattering. 

An atmosphere that is optically thin for free-free absorption but thick for 
electron scattering has a rather novel consequence. If the electrons in this 
atmosphere are at a temperature T, then photons which come into thermal 
equilibrium with them by elastic Thomson scattering will acquire a Boltzmann 
energy distribution (not Planckian) given by 

HN 

— oc v 2 exp (-hr'/kT). (43) 

Monochromatic photons passing through a plasma of temperature T transfer 
an average energy per collision 

<AE e > = (hn/mo c 2 ) (hn - 4kT). (44) 

When hn « 4kT, the photon energies get boosted by the collisions and the 
radiation is “Comptonized.” If hn > 4kT, the photons heat the gas until 
kT = hn/4 and then emerge with the Boltzmann distribution, for which 

<hy> = 3kT, <(hn) 2 > = 12 (kT) 2 , (45) 

and hence the expectation value for energy transfer (given by Equation 44) 
goes to zero. 

For a black body source 

<hi>) = 2.68 kT BB , <(hn) 2 >= 10.1 (kT BB ) 2 , (46) 

the condition of zero average energy transfer would (see Equation 44) imply 
that 


T = < (hn) 2 > / (4 < hn> = 0.94 T BB . (47) 

Finally, we compare the Boltzmann and Planck distributions at the same T: 
(dN/dJ>) BolUmann 


(dN/dn) 


•= [1 - exp (-hn/kT)] ; 


(48) 


Planck 
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This shows that they deviate from each other significantly only when 

ku<kT. 


VI. SOME OBSERVATIONS 

Dr. Rose has already discussed, in Chapter 1, two well studied X-ray pulsars 
in eclipsing binary systems,’ Her X-l and Cen X-3. Most of the X-ray data 
used in the interpretation of their binary behavior are temporal in character 
(e.g. Doppler variations of the pulsation period) and indicate that the com- 
pact member in each of these cases is a neutron star. We shall now empha- 
size the information provided by the X-ray spectra themselves. In addition 
to Her X-l, we shall describe some spectra obtained by the Goddard X-ray 
group for two compact objects which appear quite different, Cyg X-l and 
Cyg X-3. 

The data we shall describe were all obtained with the same complement of 
X-ray detectors during two rocket-borne experiments with launch dates sepa- 
rated by a year (October 4, 1973 and October 3, 1974). As is shown in Fig- 
ure IV- 1 , the first of these flights examined Her X-l , Cyg X-3, and Cyg X-l 
in that order, with exposures on the order of a minute each. We note the 
remarkable constrast in temporal behavior among the sources exhibited in 
this figure. For Her X-l we see the pronounced periodic signal of the 1.24s 
pulsar. For Cyg X-3 the variations are strictly statistical (i.e., on such a small 
time scale the 4.8 hour periodicity produces negligible systematic effects). 

For Cyg X-l , rapid temporal variations are most pronounced and chaotic on 
time scales down to a millisecond (Figure IV-2), 

The spectrum of Her X-l is extremely simple (see Figure IV-3). The photon 
number spectrum follows a rather strict power law of unit index (i.e., dN/d(hi>) 
a (h^) -1 ) from somewhat below 2 keV to somewhat above 20 keV. Observa- 
tions at higher energies indicate that the harder emission might not be pulsed. 

We pursue the suggestion that the emission from Her X-l comes from accre- 
tion onto the magnetic poles. Due to the suppression of Coulomb collision 
losses in the environment of intense magnetic fields, the particles falling onto 
the neutron star are stopped by nuclear collisions. The resulting penetration 
amounts to an optical depth for usual Thomson scattering of 10. However, 
as we have seen in the last section (see Equation 41), Thomson scattering 
along the field is also highly suppressed (i.e., along the dipole axis for regions 
near the poles) for hv « (1 1 .6 keV)H 12 • Hence, we expect the intensity 
near the dipole direction to vary as follows: 
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exposure samples of each source obtained during the same rocket-borne experiment on October 4, 1973 (GSFC 
Flight 13.010). The inset shows the overall flight profile in counts per 2-second interval versus time after launch 
(in seconds). The arrows indicate when the 15-second data samples occurred relative to the full exposure. 
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Figure IV-2. High time resolution data on Cyg X-l obtained in GSFC Flight 1 3.010. The lower portion is an 
expansion of the peak of the “flare" at about 318.5 seconds after launch. The shaded elements correspond 
to statistically significant bursts of millisecond duration. 
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ENERGY (keV) 

Figure IV-3. The average spectrum for Her X-l obtained during GSFC Flight 13.010. 
The open symbols are taken from the argon-filled counter, and the closed symbols 
from the Xenon-filled counter. The dashed curve corresponds to attenuation of a 
pencil beam in the ordinary mode whereas the solid line corresponds to the extra- 
ordinary mode (with r $ = 10 and Hj 2 = 10, in Equation 41). 
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r s / \2 

I(h,, 0) « I 0 (hr) exp J sin 2 0 - r s (49) 

We see that the exponential factor only becomes a significant attenuator for 

sin0>[2/r s ] 1 / 2 (50) 

and/or 

hr>[(11.6keV)H 12 ]/T $ 1/2 . (51) 

For r s = 10, Equation 50 implies 0 > 27° , consistent with the half-width of 
the pulsed beam. For H 12 = 10, Equation 51 implies that spectral attenua- 
tion sets in at hv > 37 keV. This leads us to believe that the high energy cut- 
off behavior for Her X-l at about 24 keV may be obtained via Equation 41, 
suitably fine tuning the parameters r s ^ 10 and H 12 *** 10. A source spectrum 
of the form Iqoc^)* 1 is not surprising for hv small compared with any of the 
possible characteristic driving energies (e.g. > 100 keV for kT, 1 1 .6 keV(H 12 ) 
or free-fall bremsstrahlung end point). As shown in Figure IV-3, we obtain a 
remarkably good fit for attenuation that is dominated by the extraordinary 
mode (i.e., minus sign for Equation 42), which is to be expected for r s (hi>)/ 

(1 1.6 keV) H 12 > 1 . 

Figure 4 shows the spectrum obtained for Cyg X-3 on October 4, 1973. Con- 
sidering data below 10 keV, it is well fit by a blackbody spectrum of kT = 1 .24 
keV (a high energy tail above 10 keV cannot be excluded). Normalized to a 
distance of 10 kpc, the apparent source luminosity is L(1 .7 - 40 keV) = 

8.8 X 10 37 ergs/s, whereas the blackbody luminosity at the inferred tem- 
perature (T 7 = 1 .4) is L = 7.2 X 10 37 ergs/s. Since 10 kpc is a likely lower 
limit to the distance, we note that the luminosity is close to the Eddington 
limit. Using Equation 15, we obtain 

M >0.7. 


From Equation 1 7 we obtain the blackbody radius as 
(1.9)2li /2 

*6— ~ = 1.84 1.6. 


Hence, the mass and radius of the blackbody are comparable to those expected 
for a neutron star. 
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Figure IV-5 shows the spectrum obtained for Cyg X-3 one year later (October 
3, 1974). The differences may be summarized as follows: 

(1) At energies of about 4 keV the new spectrum is depleted of photons by 
an order or magnitude, but at energies above about 1 5 keV it becomes 
harder than the earlier spectrum. 

(2) The second spectrum shows a strong spectral feature near the energy 
region where iron K-line emission would be expected. 

Acceptable fits may be obtained for the spectrum of Figure 1Y-5 by assum- 
ing optically thin thermal emission with kT > 17 keV and a line at 6.7 keV 
with an equivalent continuum width of 1 .2 keV (i.e. for kT = 17 keV, the 
power in the line is 7% of the total power in the continuum). The apparent 
source luminosity is L(1 .7 - 40 keV) = 4.2 X 10 37 ergs/s. However, this is 
a lower limit since the spectrum may still have an appreciable contribution 
from beyond 40 keV. Further, since this observation was taken near zero 
phase of the 4.8 hour periodicity (i.e., at an intensity minimum) it needs to 
be normalized upward for comparison with the earlier observation which 
occurred at a phase 0.8. We conclude, therefore, that both luminosities are 
about the same and are at or very near the Eddington luminosity. As we 
have seen in Section IV, the accretion rates and density profiles may change 
appreciably near the Eddington limit, thereby producing a variety of spectra 
at the same luminosity. 

To what extent could the optically thin spectrum be considered that for the 
source of the optically thick spectrum seen through an optically dense atmos- 
phere? It is interesting to note, in this connection, that the optically thick 
spectrum of Figure IV 4 could be as well fit by a Boltzmann distribution (see 
Equation 43) as by a blackbody (both at the same temperature), since there 
were in fact no observations below kT = 1 .24 keV, and even if there had been 
they would have been severely modified by interstellar absorption effects. 

As described in the last section, the atmosphere for a source near the Edding- 
ton limit is likely to be optically thick for electron scattering, but thin for 
free-free absorption. 

What sort of inputs are permissible for an emerging Boltzmann spectrum at 
kT = 1 .24 keV? From Equation 45 we note that 

<(hF) 2 )/<(h^)> = 4kT«5 keV. 

Hence, for example, a single line at 5 keV would be one possibility, but that 
is not what was seen. For a spectrum of the form 
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Figure IV-4. The spectrum obtained for Cyg X-3 on 4 October 1973 
(GSFC Flight 13.010). 
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dN 

dhi' CC 


Chi')" 1 exp (-hp/kT') 


we have 


< (h v) 2 ) 
<h^> 


= kT'. 


(52) 


This implies kT' = 5 keV, clearly less than the lower limit of 17 keV obtained 
for the spectrum of Figure IV-5. For a power-law spectrum of number index 
T and end-point E 0 w'e find 


<(hp)2) _/2-r\ 

(hr-) \3-rJ*o- 


(53) 


Taking kT for the Boltzmann Equation, then Equations 45 and 53 imply 


2- T 

3- T 


4kT 

E 0 ' 


(54) 


For kT = 1 .24 keV and Eq = 40 keV this implies V ^ 1 .9, which is steeper 
than the optically thin spectrum observed. Finally, let us consider a black- 
body as the source spectrum. From equation 46 we find that the correspond- 
ing blackbody source would be one for which kT = 1.4 keV, not much differ- 
ent from that already considered. 


Although the spectrum of Cyg X-l is notorius for variations, it is in fact the 
most stable in X-ray astronomy. From the first detection of radio emission 
in March 1971 to the apparent transition in April 1975 (i.e., ~4 years) the 
spectrum remained remarkably constant. For example, the two spectra 
(Figures IV-6 and IV-7) obtained one year apart with the same instrument 
are both statistically significant and indistinguishable for the band covered 
(—1.5- 40 keV). They both yield a power-law spectrum of number index 
1.57, presumably highly characteristic of this emission from an accretion 
disk. Normalizing to a minimum distance of 2.5 kpc, we find the luminosity to 
be L > 2.9 X 10 36 E?- 43 ergs/s, where Eq is the upper end of this power 
law spectrum. Since Eq > 40 keV, we obtain 

L> 1.4 X 10 37 ergs/s. 

Since this source is probably several solar masses (i.e., a black hole), the 
luminosity might very well be much below the Eddington limit. 
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As Dr. Rothschild will discuss in detail in Chapter 6, the bulk of this emis- 
sion may be characterized by randomly occurring events (i.e., shot noise) of 
comparable magnitude, such that we might express the average luminosity 
as 


L= WX 

where W is the mean energy per event and X <; 20 s' 1 is the average rate 
of events (at least during the two observations considered here). Hence we 
conclude that 


W> 7 X 1035 er g S ^ 

Following the arguments of Section III, we conclude that 
W e < 1 .5 X 10 33 L 8 /Hj 2 . 

Hence, the lower limit to the energy per event places an upper limit on the 
magnetic field associated with electrons radiating cyclotron radiation in this 
band. This yields Hj 2 < 2 X 10 -2 . 

For bremsstrahlung radiation we have the condition 

W< 6 X 1 0 29 x/Lg ■ (Vol.)t/2 mc 2 

(1 -/3 2 /3)!/2 

Since the duration of these events is a fraction of a second and 0 < 0.5 we 

have that the emitting volume <3 X 10 30 cm 3 . Hence 

W < 10 38 [m/m (electron) ] ergs, 

and we do not obtain a significant restriction (i.e., the shot noise events are 
well under 10 3 ® ergs each). However, the millisecond bursts have an energy 
of ^ 1 .4 X 10 35 ergs each. For millisecond bursts the maximum emitting 
volume is down by a factor of about 10“ 9 relative to “one second” events, 
and we obtain the corresponding constraint 

W< 3 X 10 33 [m/m (elecUon) ] ergs. 

This tells us that the radiating particles must be more massive than electrons 
(for example, the energy content in protons could be as high as 6 X 10 36 ergs). 
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As Dr. Rothschild will discuss in detail in Chapter 6, the bulk of this emis- 
sion may be characterized by randomly occurring events (i.e., shot noise) of 
comparable magnitude, such that we might express the average luminosity 
as 

L = WA 

where W is the mean energy per event and A ^ 20 s" 1 is the average rate 
of events (at least during the two observations considered here). Hence we 
conclude that 


W>7 X 10 35 ergs. 

Following the arguments of Section III, we conclude that 
W e < 1 .5 X 10 33 Lg/H^ 2 . 

Hence, the lower limit to the energy per event places an upper limit on the 
magnetic field associated with electrons radiating cyclotron radiation in this 
band. This yields H 12 < 2 X 10“ 2 . 

For bremsstrahlung radiation we have the condition 

W< 6X 10 29 \/L^ — (Vol.) 1 / 2 me 2 . 

(1 - 0 2 / 3 ) l / 2 

Since the duration of these events is a fraction of a second and 0 < 0.5 we 
have that the emitting volume <3 X 10 30 cm 3 . Hence 

W< 1038 [m/m (electron) ] ergs, 

and we do not obtain a significant restriction (i.e., the shot noise events are 
well under 10 33 ergs each). However, the millisecond bursts have an energy 
°f £1-4 X 10 35 ergs each. For millisecond bursts the maximum emitting 
volume is down by a factor of about 10" 9 relative to “one second” events, 
and we obtain the corresponding constraint 

W< 3 X 10 33 [m/m (electIon) ] ergs. 

This tells us that the radiating particles must be more massive than electrons 
(for example, the energy content in protons could be as high as 6 X 10 36 ergs). 
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So far we have restricted ourselves to some data below 40 keV, the region 
of the spectrum accessible from rocket-borne proportional counters. The 
larger exposures and extended response needed for examining the harder 
emission from compact X-ray sources is best supplied by satellite-borne scin- 
tillation counters and solid state detectors. Results from one such pioneering 
experiment are summarized in Figure IV -8. These spectra represent long-term 
averages and give us a rather simple picture for characterizing the gross spec- 
tral properties for the X-ray emission from sources in binary systems. Ex- 
pressed as energy spectra (i.e., keV/keV), they are remarkably flat at low 
energies (< 10 keV) and then fall off rather abruptly at a few tens of kilovolts, 
with essentially no detectable emission for energies above about 50 keV. As 
we noted in Section II, such energy scales are to be expected for compact 
objects such as neutron stars and black holes. 
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CHAPTER V 

GAMMA RAYS FROM COMPACT OBJECTS 


Hakki Ogelman* 

Goddard Space Flight Center 
Greenbelt , Maryland 


I. INTRODUCTION 

In contrast to radio and x-ray emission studies from compact objects, the 
gamma ray observations, at best, can be classified as “underdeveloped.” If 
we examine the productivity of these regions of the electromagnetic spectrum 
we see that about 150 compact objects (pulsars) have been discovered in the 
radio, and about 40 in the x-ray region (mass-accreting neutron stars). The 
gamma ray region, on the other hand, has two compact objects on its inven- 
tory. Furthermore, of these two, the Crab Nebula and Vela pulsars were 
initially discovered in the radio region. We need not be too apologetic, how- 
ever, since in the rest of the electromagnetic spectrum the Crab pulsar is the 
only object seen in other than the radio, X-ray, and gamma ray regions. 

Two of the strongest gamma ray sources which the present level of experi- 
mental sophistication allows us to see are compact objects. This is similar to 
the history of x-ray astronomy, in which the strongest sources detected ini- 
tially, are now interpreted as mass-accreting compact objects. In contrast, 
the two strongest radio sources first discovered were a supernova remnant 
and another galaxy. Today pulsars are a small set of thousands of radio 
sources, whereas compact x-ray sources constitute about 30% of all the 
sources observed in x-rays. The conclusion we can draw from this type of 
data is that as gamma ray astronomy techniques develop we can hope to see 
compact objects emerging as the main component of sources emitting in this 
energy region. The expectations of the future, then, justify a discussion of 
the implications of just two gamma-ray sources in our understanding of the 
physical processes that take place in compact objects. 

We may ask, what is special about gamma rays and compact objects? In the 
previous lectures we have heard about radio emission from pulsars and x-ray 
emission from neutron stars in close binary pairs. These emissions reflect 
quite different conditions. The radio emission from pulsars signifies a non- 
thermal process requiring strong magnetic fields and very energetic electrons. 
In such emission processes we have to have a very tenuous medium, vacant of 

*NAS-NRC Senior Postdoctoral Research Associate at GSFC, on leave from Middle East 
Technical University, Ankara, Turkey. 
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any significant material, so as to prevent the energetic electrons from ther- 
malizing in the medium. In contrast, the proposed emission mechanism for 
x-ray -emitting compact objects is thermal emission of material that gets heated 
up as it falls into the gravitational potential well of a neutron star. As we have 
seen, the rate of energy release in this well if related to the size of neutron star, 
has an effective temperature that causes emission in the x-ray binaries. We 
would then expect to encounter a material starved, magnetic field rich emis- 
sion process. Our first response would be that under these circumstances, 
nothing very special should happen in the gamma ray region. The electromag- 
netic emission spectrum in magnetic fields more or less follows the trend in the 
electron spectrum; if anything, it steepens with respect to the normal trend at 
high energies. Since, a priori , we expect the energetic electrons in pulsars to 
have some form of power law distribution that gives rise to all the emission, 
we should not expect any special enhancement in our simplistic approach. 
Indeed, the spectrum observed for the Crab Nebula pulsar does justify this 
notion; however, the Vela pulsar is not observable between the radio and 
gamma ray regions. Actually, a source is observed in the 2-10 keV range 
which appears to show pulsations occasionally. Certainly, the two objects 
behave quite differently in the region between 10” 3 to 10 7 eV. 

One process that lends special importance to gamma rays around 1 00 MeV is 
the production of n° mesons in high energy interactions of nucleons with mat- 
ter. Subsequently the tt° mesons decay into two gamma rays with ^70 MeV 
energy each in the tt° reference system. The geometry of the collisions with 
respect to the observer doppler-broaden this in the laboratory frame; however, 
the peak of the distribution is always centered around 70 MeV. Effectively, 
then, gamma rays around this energy are the characteristic quanta that high 
energy nuclei emit in material, just as thermal activity is indicated by emis- 
sion lines of the atoms. However, as we shall discuss later, the pulsars are 
believed to have very little material in their magnetospheres and the 7r° pro- 
duction mechanism does not appear to be effective. 

Turning back to the material starved, electromagnetically rich pulsar mag- 
netosphere, we can speculate on how we might get special emphasis in the 
emission around the gamma ray region. As a beginning, consider observations 
of secondary gamma rays at sea level in the earth’s atmosphere. If we per- 
formed such an experiment, we could notice an excess of gamma rays around 
100 MeV. Although the details of the individual interactions are complex we 
can guess roughly why the energetic photons cluster around 100 MeV. Ini- 
tially we start on top of the atmosphere with a very energetic gas of cosmic 
rays. The atmosphere is many interaction mean-free-paths thick both to this 
hot gas and to its secondary components. The energetic protons quickly 
interact, producing 7r -mesons which immediately decay into electrons, positrons, 



GAMMA RAYS FROM COMPACT OBJECTS 


101 


and gamma rays. The electrons predominantly lose their energy by brems- 
strahlung to photons of comparable energy, and the photons pair-produce 
and create electrons with similar energies. The electromagnetic cascade thus 
multiplies and grows. When the average energy drops to around 100 MeV, the 
electrons start losing their energy predominantly by ionization and the pho- 
tons by Compton processes. At this stage the electron and photon compo- 
nents can no longer sustain each other and the shower stops growing. Further- 
more, around 100 MeV the photons have their minimum cross section and 
they can emerge more readily at sea level. We may expect a similar situation 
in the pulsar magnetospheres. In this case, although there is no material, the 
electrons and photons are coupled to each other through the intense magnetic 
Fields. In particular, the photons anihilate with the characteristic photons of 
energy treB/mc of the magnetic fields, and produce electron-positron pairs. 

The threshold of this process is around 25 MeV on the surface where B ^ 10 13 
eV at the speed of light cylinder. In a simplistic approach then, we can sus- 
pect that energetic electrons and photons will create a shower wherein the 
emerging flux of photons would have energies in the range indicated, which 
conveniently falls into the gamma ray range; thus we have a gamma ray biased 
emission. 

With these general considerations we shall proceed to examine the gamma ray 
data in Section II, and gamma ray production in pulsars in Section III. In 
Section IV we shall try to estimate the relation of gamma rays from pulsars 
to the total galactic gamma ray emission. 

The discussion of the data will be limited to the energy region 35 MeV to 
100 MeV and we shall primarily discuss the results from the SAS-2 satellite, 
since at present it offers the most complete self-consistent survey of the 
galaxy in this energy region. 


II. OBSERVATIONS 

The First compact source discovered in gamma rays is the pulsar NP0531 in 
Crab Nebula. The energy range in which this object is detected extends from 
radio photons to 10 9 eV and possibly up to 10 12 eV. The electromagnetic 
spectrum in the range 10 5 to 10 9 eV is shown in Figure V-l. The dashed line 
represents the intensity of the Crab Nebula supernova remnant. As can be 
inferred from the Figure, while in the 1-10 keV range the pulsar emission is 
about 10% of the total nebulosity; above 10 MeV the pulsar emission domi- 
nates over the supernova remnant. Figure V-2, shows the phase distribution 
of gamma ray events above 35 MeV as measured by SAS-2 detector folded 
modulo of NP0531 radio period. Also shown on the same Figure are the 
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NP0532 

GAMMA RAYS >35 McV 



TIME IN FRACTIONS 
OF A PULSE PERIOD 


Figure V-2. The phase distribution of gamma ray events above 35 MeV as measured 
by SAS-2 detector folded modulus of NP0531 radio period. The arrival times of the 
main (M) and secondary (S) radio pulses corrected for interstellar dispersion are 
marked (Source: Kniffen et al M 1974). 
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arrival times of the main (M) and secondary (S) radio pulses, corrected for 
interstellar dispersion. The gamma ray distribution shows two peaks which 
are coincidental with the corresponding radio features. The statistical prob- 
ability of getting this phase distribution is around 10"^ . The flux of pulsed 
gamma rays above 35 MeV is (6.2 ± 2.8) X 10"** photons cm" 2 s~* , and above 
100 MeV it is (2.2 ± 0.7) X 10" 6 photons cm* 2 s" 1 . 

Since the Crab pulsar, with its 33 ms period is the fastest known pulsar, it is 
natural to go down the line and examine other pulsars. The second fastest 
pulsar, PSR1913 + 16, has a period of 59 ms. However, it is a binary pulsar 
and its period and orbit parameters are not known accurately enough to per- 
mit gamma ray pulsation analysis. The third fastest pulsar is PSR0833-45, 
the Vela pulsar, which is located four times closer than the Crab pulsar. In- 
deed, SAS-2 results of the gamma ray emission of Vela show a strongly pulsed 
flux with the 89 ms period of the radio pulsar. Figure V-3 shows the phase 
plot of these gamma rays. One strange feature of the gamma ray phase dis- 
tribution is that, unlike its radio counterpart, it shows two peaks. For the 
radio observations there is only an upper limit to a second pulse away from 
the main pulse amplitude. In the gamma ray region the second pulse is almost 
as intense as the main one. Moreover, the radio peak does not coincide with 
either of the gamma ray peaks. The pulsed flux above 35 MeV is (15.1 ± 2.4) 

X 10"* photons cm“ 2 s" 1 . In the wavelength region between radio and gamma 
rays, the features of the Vela pulsar differ from that of the Crab Nebula pulsar. 
The intrinsic luminosity ratio of Crab to Vela pulsars, Lc ra b/Ly ela , is 8 around 
100 MeV, 80-100 around X-ray energies, and >5000 at optical wavelengths. 

As far as gamma rays are concerned there is hardly any difference in the phase 
plots of these pulsars. Figure V-4 shows the pulse shape ot these two pulsars 
in various energy bands. 

The various observational features of these two pulsars are summarized in 
Table V-l. 

In addition to the positive results we can set upper limits to a variety of pul- 
sars that have accurate enough period and period derivative information. The 
preliminary results of such an effort are summarized in Table V-2. What has 
been done in the analysis is the following: Gamma rays arriving from the 
proper position of the pulsar, within about 2a of the angular accuracy of the 
measurements, were folded modulus of the pulsar period at the solar system 
barycenter. The phase distribution was then examined in bins of 0.1 period 
phase. Since we do not have, a priori , a good estimate of what pulse shape 
to expect (remember that the phase plot of gamma rays from Vela pulsar 
does not look like its radio counter part), we merely check the y} analysis 
of the phase distribution’s compatibility with random data. Since none of 
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PSR 0833-45 
GAMMA RAYS > 35 MeV 



FEB. 15-21 , 1973 


FEB. 21-27, 1973 



TIME IN FRACTIONS 
OF A PULSE PERIOD 


Figure V-3. Phase plot of gamma rays from Vela pulsar using 89 ms period 
(Source: Thompson et al., 1975). 
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NP 0532 


PSR 0833-45 




PER 100 -33 MSEC 

TIME IN FRACTIONS OFA PULSE PERIOD 


Figure V4. Pulse shape of the Crab and Vela pulsars at various energy ranges. 
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Summary ofNP0531 and PSR0833-45 Data 
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Table V-2 

Summary of Data on Pulsars with Upper Limits to 
7-Ray Flux above 35 MeV 


Pulsar 

X 2 (9 deg. 
of freedom 

3a limit of 
pulsed flux 
(cm -2 s" 1 ) 

P(sec) 

P/2P (yr) 

D (pc) 

7 Luminosity^ 
upper limits 
(photon s -1 ) 

P0525 + 21 

7.3 

1.1 X 10- 6 

3.745 

1.5 X 10 6 

1460 

3.0 X 10 32 

P0540 + 23 

13.6 

1.2 X 10- 6 

.246 

2.5 X 10 5 

1980 

5.7 X 10 37 

P0611 + 22 

5.9 

15.5 X 10- 6 

.335 

9.0 X 10 4 

2720 

1.1 X 10 38 

P0943 + 10 

16.2 

0.7 X lO' 6 

1.098 

4.9 X 10 6 

300 

7.8 X 10 35 

P0950 + 08 

4.4 

0.6 X 10- 6 

.253 

1.8 X 10 7 

80 

4.9 X 10 34 

PI 133 + 16 

17.5 

0.6 X 10-* 

1.188 

5.0 X 10 6 

138 

1.4 X 10 35 

P1237 + 25 

16.1 

0.4 X 10-* 

1.382 

2.3 X 10 7 

296 

4.5 X 10 35 

P1451 + 68 

7.9 

3.7 X 10- 6 

.263 

>1.4 X 10 6 

222 

2.2 X 10 36 

P1642 -03 

16.3 

0.6 X lO"* 

.387 

3.4 X 10 6 

125 

1.2 X 10 35 

PI 706- 16 

7.0 

1.0 X 10-* 

.653 

1.7 X 10 6 

131 

2.1 X 10 35 

P1727-47 

6.2 

1.9 X lO" 6 

.830 

7.7 X 10 4 

3960 

3.6 X 10 38 

PI 747 -46 

1 1.3 

1 

1.8 X 10- 6 

.742 

1.7 X 10 s 

560 

7.0 X 10 36 

P1749 -28 

4.6 

1.4 X lO* 6 

.562 

1.1 X 10 6 

1310 

2.8 X 10 37 

PI 8 1 8 - 04 

13.6 

0.9 X 10" 6 

.598 

1.5 X 10 6 

2590 

7.6 X 10 37 

P 1 9 1 5 + 13 

3.4 

1.5 X 10-« 

.195 

4.3 X 10 s 

2480 

1.1 X 10 38 

P1919 + 21 

7.1 

1.4 X lO” 6 

1.337 

1.6 X 10 7 

328 

1.8 X 10 36 

P1929 + 10 

10.1 

1.7 X lO" 6 

.227 

3.1 X 10 6 

89 

1.6 X 10 35 

P2020 + 28 

3.0 

2.1 X 10-* 

.343 

2.8 X 10 6 

651 

1.1 X 10 37 

P2021 + 51 

6.7 

2.3 X lO" 6 

.529 

2.8 X 10 6 

624 

6.3 X 10 36 

P2 1 1 1 + 46 

5.2 

1.4 X lO' 6 

1.015 

2.3 X 10 6 

317 

1.7 X 10 36 

P2217 + 47 

7.0 

1.1 X lO" 6 

.538 

3.1 X 10 6 

1250 

2.0 X 10 37 


^The emission solid angle of the pulsar assumed to be 0.4* steradians. 
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of the plots showed a x 2 more than 17.5 (Probability >.05 of occurring 
by random chance on a x 2 distribution with 9 degrees of freedom), we 
have listed in Table V-2 the upper limits so the gamma ray luminosity for 
these objects together with assumed distances and the actual values of x 2 . 

The upper limits were roughly calculated by assuming that if the pulsar flux 
existed it would be in one bin and the limit is 3 a of the mean per bin of the 
phase plot. In the future a more refined analysis with actual comparison of 
radio arrival times may reveal that these upper limits are actually close to 
positive measurements, but for the time being we shall only treat them as 
upper limits. Since even these upper limits - when expressed as a limit to 
the luminosity - contain large errors due to the errors in distance determina- 
tion, they should be treated cautiously. In the subsequent discussion we shall 
only use them as a general check on the extrapolated behavior of the Crab and 
Vela pulsars. 


III. PULSAR ENVIRONMENTS AND GAMMA RAY PRODUCTION 

As we have seen from the previous discussions on compact objects, they 
have the following typical parameters: M ~ 1M@, R ~ 10 6 cm, B ~ 10 12 
gauss; the rotation periods P cover the range .03 to 3 seconds. 

From considerations of electrodynamics, a rotating dipole has to set up a 
large voltage between its surface and the speed of light cylinder (where the 
magnetic field lines can no longer loop around and so close on the surface). 
This is simply a consequence of the fact that near the pulsar the plasma has 
to co-rotate with the neutron star, and far out it cannot do so. The voltage 
thus developed does not depend very much on what the dipole axis inclina- 
tion is with respect to the rotation axis. We can estimate the magnitude of 
the voltage simply by using Faraday’s law of induction: 


<t> (voltage) ~ — 


1 d f 

c dtj B ' ds ' 


If we take dt ^ I/O , B • dS at BR 2 , then 

r 0 ftBR 
C R c 


(cgs). 


(III-l) 


(III-2) 


Substituting typical neutron star parameters, we see that the voltage ranges 
between 10 12 and 10 10 volts/ cm near the surface. This potential difference 
is very attractive for particle acceleration. However, what fraction of this 



110 


RECOGNITION OF COMPACT ASTROPHYSICAL OBJECTS 


voltage is available for particle acceleration is difficult to say, since the man- 
ner in which the magnetosphere is covered with the network of currents that 
create this potential has a great bearing on this question. 

We can also estimate the magnitude of the currents and hence the density n 
of the charge-separated plasma that will provide these currents. Again let us 
use elementary electrodynamics. From Gauss’s Law, we get 

V • E = 47 rp = 47ren. (HI-3) 


Using E from Equation (III-2), we get 


UB 


n 


47rec 


(III-4) 


With typical parameters we have charge -separated densities around 10 13 cm" 3 
at the surface and 10 4 cm" 3 at the speed of light cylinder for a Crab Nebula 
type pulsar where SI =* 200 s" 1 . Although these densities are calculated for 
currents, the regions which do not partake in these currents could not have 
densities very different; at least we can use this n as a lower limit. 

It is interesting to take the surface density of 10 13 cm" 3 and integrate it to 
infinity assuming that it falls off as R" 3 . The result is a surprising sum of 
~ 5 X 10 18 particles cm" 2 . If we assume neutrality, we can multiply it by 
the proton mass to get the total column density of 10 micrograms - cm -2 ! 

There is an independent way of calculating this column density based on the 
pulsar dispersion measurements. The radio signals from pulsars show shift in 
arrival times as a function of frequency which is interpreted as due to the 
dispersion in the interstellar medium. Effectively the dispersion measure 
gives us Jn e • d£ between us and the radio emitting region on the pulsar. Most 
of the dispersion measures ranging from 10 to 500 parsec cm" 3 can be accounted 
for by travel in the interstellar medium. Let us estimate that less than 5 parsec 
cm" 3 of this could be intrinsic to the pulsar. This gives us 


I n e dC < 1.5 X 1019 electron cm" 2 
•'pulsar 

vhich is in accord with the estimated value of 5 X 10 18 cm -2 . It seems that 
pulsars have very little material above their surfaces, something around 10 to 
20 micrograms cm -2 , a factor of 10 8 less than that Dr. Boldt has estimated 
for mass accreting neutron stars. 



GAMMA RAYS FROM COMPACT OBJECTS 


111 


It is interesting to do a rough estimate on what mass accretion rates are ex- 
pected for pulsars, inasmuch as they are neutron stars and at distances large 
compared to the speed of light cylinder; they should appear no different than 
neutron stars of x-ray binaries. Let us assume that the pulsar has a velocity v 
with respect to the interstellar medium. The accretion radius r ac can be de- 
fined as that where the kinetic energy of the gas with respect to the moving 
pulsar is less than the gravitational energy at that distance, i.e., 


r 


ac 


2GM 

v 2 


(III-5) 


The mass accretion rate is then given by 


dM 

dt 


: pv 


_ tt(2GM) 2 p 


(III-6) 


where p is the density of the interstellar medium in which the pulsar moves. 
With M ^ 1M© , p c- 1 hydrogen mass cm -3 , and v ^ 50 km s _1 , we get 

« 7 X 10 9 g s" 1 = 7.4 X 10” 17 M© y r 1 . 


Compared to mass accretion rates of 10~ 6 - 10" 9 M® yr -1 for x-ray stars, 
this accretion rate is too small to be of any importance in generating energy. 
However, such a process may cause the enhancement of the material in the 
line of sight. To get an order of magnitude estimate of the amount of mate- 
rial between the speed of light cylinder r c and the accretion radius r ac let us 
assume p 0 - 1 atom cm* 3 gets compressed as (r ac /r c ) 2 and then it trickles 
on the neutron star surface. Then 


~ r ac J @0 * (HI-7) 

~ ‘ac 

With the previously assumed values and — 200 s” 1 , we obtain 10" 18 atoms 
cm” 2 . Again, this value is in accord with the other two methods of estimation. 

We have seen that the potential gradients in the magnetospheres of pulsars is 
sufficient to accelerate protons above the ir° production threshold. However, 
we are also forced to the conclusion that there is very little material above 
the surface of the neutron star that could make 7r° mesons out of these photons. 




112 


RECOGNITION OF COMPACT ASTROPHYSICAL OBJECTS 


As a matter of fact only about 10" 5 /72 of the protons would interact. To 
get 10^5 er g S S -1 j n gamma rays would require a proton energy luminosity of 
more than 10 43 ergs s*" 1 which is a factor of 10 5 above the rotational energy 
loss rate of the pulsar. 

One may speculate on the possibility of n° production by the protons of the 
magnetosphere accelerated towards the surface of the neutron star; in this 
case we do have plenty of target material. At about 1 GeV proton kinetic 
energy, the cross section times multiplicity of 7i° production is about 10' 26 
cm 2 , at higher proton energies this product increases as approximately the 
0.5 power of the proton kinetic energy. Roughly, GeV protons convert a 
few percent of their energy to gamma rays of interest. As the energy of the 
proton goes up, the efficiency of energy conversion to gamma rays around 
100 MeV drops (^as T -1 / 2 ). Even in the optimum range of proton energies 
we need 10 37 to 10 38 ergs s”* 1 , or 10 39 to 10 40 protons s“ l . This is about a 
factor 10 6 larger than the estimated primary electron flux from the pulsar. 

Let us leave protons, pions, and nuclear processes temporarily and look at the 
possibility of electromagnetically generating the gamma rays. 

What happens to energetic electrons and photons in a medium such as that 
encountered in the pulsar magnetosphere? 

Electrons can interact mainly in one of two ways yielding end products that 
are again photons. One process is synchrotron radiation, in which the energy 
of the radiated photon in eV is 

E 7 = 5 X 10- 20 B| Eg (eV). (III-8) 

The energy loss rate in eV per second, is given by 

dE e 

— = 5 X 10- 15 E? (eV)B 2 (III-9) 

dt e 

with B in gauss. Or, we can write the mean pathlength of the electron as 
cE 

\ “ OEjTt “ 6 x 1024 E e* B ' 2 cm - OH-10) 

The other process that affects the electron is curvature radiation. In this case 
the energy of the radiated photon is 
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E y = 6.5 X 10- 23 Eg (eV) R" 1 e V, (IIM 1) 

where is the radius of curvature of the trajectory in centimeters. The 
energy loss rate is 


dEg 

dt 


= 5 X 


10-20 E 4 r -2 eV 


and the mean-free-path is 


(III-l 2) 


X c ^ 6 X 10 29 E; 3 R 2 cm. (IIM3) 


We could also calculate same type oi formulas for inverse Compton type in- 
teractions. But since, at the surface, the magnetic field energy density H 2 /87 r 
is about 4 X 10 22 ergs/cm 3 and photon energy density is around 2 X 10 14 
ergs cm" 3 even if the object is radiating at the Eddington limit, for all practi- 
cal purposes all the main energy loss mechanisms are going to be due to the 
magnetic field. 

Now, without going into the details of the processes let us try to understand 
what will happen to energetic electrons trying to get out of the surface of the 
pulsar. 

Notice that if we ask what sorts of photons will a 10 12 eV electron radiate 
through synchrotron processes near the surface, where B ^ 10 12 gauss, we 
get an unrealistic value like 5 X 10 16 eV; this is much larger than the electron 
energy. The conflict is due to the fact that Equation (III-8) is a classical for- 
mula whereas we are really dealing with quantum electrodynamical regimes 
at such electron energies near the surface. The parameter T that separates 
the classical regime from the quantum region is given by 


T 


E JB_ 
me 2 ®cr 


(1IM4) 


where the critical magnetic field B cr is that at which the energy of the photon 
with Larmor frequency he B/mc exceeds the rest mass energy of the electron, 
n^c 2 : 


®ci ~ ' 


m 2 c 3 


= 4.4 X 10 13 gauss. 


eh 


(IIM 5) 
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With these parameters we can express T as: 

T = 4.5 X 10- 20 E(eV) B(gauss) (IIM6) 


If T is less than 1 we are in the classical regime, and otherwise in the quantum 
mechanical regime. In our example of E = 10 12 eV and B - 10 12 gauss, the 
parameter T ^ 4.5 X 10 4 , well into the quantum mechanical regime. Syn- 
chrotron emission in such a regime is similar to bremsstrahlung emission in 
solid material. In contrast to the photons from a bremsstrahlung spectrum 
which have equal amounts of energy per logarithmic interval up to the elec- 
tron energy, the magnetic bremsstrahlung has a broad peak at 



(III-17) 


which approaches E e at T » 1 . In the quantum regime the characteristic 
length over which electromagnetic interactions take place is of the order of 
the Compton wavelength, X 10" 11 cm. Again if we apply the condition 
of E e a* 10 12 eV and B ^ 10 12 gauss to Equation (III-10) we see that the 
mean-free-path of the electron is only 15% of the Compton wavelength. The 
implications are then that we will get cascades like 

e - + B -► e" + electrons + positrons + photons + B 


in a very small region of the order of the Compton wavelength 


Notice that the condition for pair production, 

heB 

(iii ' 18) 

is equivalent to 

^7 B 

— 3T>1, (III-19) 

M e c 2 “cr 

which is identical to the condition that electrons radiate in the quantum 
mechanical phenomenon with the same Feynman diagram as the bremsstrah- 
lung. 
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The essential import of all of this is the following: 

As long as the energies E of the photons and electrons are greater than 
2 X 10 19 /Bj^ eV, there is a cascade growth that reaches maximum in a very 
short length when the average energies of the electrons and photons drop 
down to a critical value E c = 2 X 10 19 /Bj^. For surface fields of 10 12 gauss, 

E c ~ 20 MeV. Hence any acceleration process on the surface cannot yield 
energizing electrons or photons above this energy, and most of the emerging 
photons will have energies around 20 MeV. Unless there is more acceleration, 
these shower electrons would radiate the remainder of their energy with a 
spectrum similar to the magnetic bremsstrahlung of a 20 MeV electrons in 
fields slightly less than the surface fields of the pulsar. The shower photons 
would emerge from the pulsar with their ~20 MeV energy. Indeed we thus 
have a mechanism for generating a large number of photons around the energy 
of the gamma rays of interest, just as most photons in extensive air showers 
have the critical energy appropriate for air. 

IV. CONTRIBUTION OF PULSARS TO THE GALACTIC 
GAMMA RAY EMISSION 

Generalization based on two observations may lead to erroneous conclu- 
sions, however, scientific knowledge follows the extrapolation of existing 
knowledge and we cannot resist the move. Furthermore, we do have a num- 
ber of upper limits that will aid our generalizations. 

The intention then, is to consider the Crab, Vela and other pulsars as the 
manifestations of the same kind of object at different ages, from the data to 
determine an empirical luminosity-versus-time curve, and to examine its 
implications. 

In Figures 5, 6, and 7 we have plotted the gamma ray luminosities for the 
Crab and Vela pulsars together with upper limits for the other pulsars as a 
function of the period P, apparent age P/2P, and the total rotational energy 
loss rate respectively. Those upper limits which have x 2 greater than 15 are 
circled as being more promising as possible positive fluxes. 

The most common interpretation of the relationship of P and P is in terms 
of the radiation losses of a rotating dipole. In this case, the energy radiated 
away is proportional to P -4 . If we assume that the radiated energy is sup- 
plied by the rotational energy loss of the pulsar, 
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Figure V-6. Gamma ray luminosity for Crab and Vela pulsars together with 
upper limits for the other pulsars as a function of apparent age (P/21*). 
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'igure V-7. Gamma ray luminosity for Crab and Vela pulsars together with uppe 
other pulsars as a function of the total rotational energy loss rate. 
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dE rot _ 47r 2 IP 
dt p3 


(IV-1) 


where I is the moment of inertia of the pulsar. Under these assumptions P 
should be proportional to . Furthermore, if the radio luminosity is pro- 
portional to the rotational energy loss rate, we expect L radio to be propor- 
tional to P" 4 . Experimentally, when we compare these measurements as 
determined in the radio region, we find that P shows no correlation with P, 
and radio luminosities follow aP ' 1 3 to P" 1 * 8 power law dependence. The 
lack of correlation between P and Pcan be interpreted as possibly the large 
variation 20) between initial parameters such as the moment of inertia 
and magnetic field strength for different pulsars. The relatively slow drop- 
off in radio luminosity as a function of period implies that in the earlier 
stages of the pulsar most of the dissipated rotational energy goes into other 
things, possibly electrons in the remnant. Progressively as the pulsar ages we 
get a larger fraction of the available energy out as radiation. 


In the same spirit let us examine the gamma ray data. In Figure 5 we see the 
plot of the logarithm of gamma ray luminosity L y versus log P. A straight 
line plotted between the Crab and Vela pulsars gives 

L 7 (>35 MeV) = 2.8 X 1 035 p-2.2 ± .5 photon S -1 (F/-2) 

Upper limits of L y to pulsars P0950 + 08, P1642 - 03, and P1929 + 10 fall 
more than an order of magnitude below this extrapolation. However, we 
should bear in mind that the similar radio luminosity plots exhibit a range of 
two order of magnitude for pulsars with the same period. The gamma ray 
data are also consistent with pulsars with periods larger than 0.2 second not 
emitting gamma rays. This seems like an unlikely conclusion since the volt- 
ages in the magnetosphere fall off as P" 1 and the density of material available 
for acceleration falls off as P 4 . In a rough model of gamma rays being the 
product of magnetospheric showers of energetic electrons, we would expect 
the number of gamma rays n 7 emitted at the shower maximum per energetic 
primary to be proportional to 


riy (at shower max) 



(IV-3) 


where E 0 is the energy of the shower initiating electron and E crit is the energy 
when photons and electrons decouple. The total gamma ray luminosity b y 
should then be proportional to: 
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n 0 


■^crit 


(IV-4) 


If we estimate that E 0 will fall off as P -1 (Equation III-2) and n 0 , the number 
of shower-initiating primaries, to fall off as P - * (Equation III-4), and that B, 
which determines E crit , does not change much, then L 7 ® P -2 . This type of 
luminosity dependence on period seems to be supported by the Crab and Vela 
pulsars. 

If we examine the L^, versus apparent age t = P/2Pplot of Figure V-6 we notice 
that Crab and Vela pulsars fit 


L y (>35 MeV) = 4 X 10 37 


pr ) j -°- 9± ' 2 


photons s" 1 . (IV-5) 


Although a few of the upper limits are slightly below this extrapolation, con- 
sidering the uncertainties in distances, Equation (IV-5) is not in conflict with 
the rest of the data. Pulsars PI 1 33 + 16 and PI 642 - 03 with x 2 values greater 
than 16 lie on this extrapolation. 

On Figure V-7 we show the distribution of Ly and upper limits to it as a func- 
tion of the spindown enegy loss rate assuming that each pulsar has a moment 
of inertia I = 1 O 45 g cm 2 . Also drawn on this figure is a line that indicates 
the case of all rotational energy loss going into production of gamma rays 
above 35 MeV. Notice that the Crab pulsar converts about 10 -4 and Vela 
10~ 3 of the rotational energy loss into gamma rays observed. The extension 
of Crab and Vela pulsar luminosities implies that when the gamma vary lumi- 
nosities are around l^t X 10 37 photons s* 1 , all of the available rotational 
energy loss is being used in generating gamma rays. It seems then, making 
use of Figure V-6, that such a power law extension must break down when 
the apparent age of the pulsar is around 2 - 6 X 10 6 years. 

The SAS-2 results indicate that, with reasonable assumptions of the galactic 
distribution of sources, the galaxy should be emitting IQ 42 photons s' 1 in 
photons above 35 MeV. What fraction of the galactic luminosity can we ac- 
count for by pulsars of all ages distributed throughout the galaxy? 

Let us first make some assumptions: 

1 . Pulsars are created during supernova explosions and the rate of generation 
is about 1 per 100 years in our galaxy. 



GAMMA RAYS FROM COMPACT OBJECTS 


121 


2. The apparent age P/2Pis a good approximation of the true age. 

3. Pulsars emit gamma rays as a function of time given by the power law 
relationship of Figure V-6: 

L,, (>35 MeV) ~ 4 X 10 37 "° 9 

4. This gamma ray emission continues until the apparent age is 5 X 10 6 years, 
and then turns off or decreases its luminosity very quickly. At this stage 
the gamma ray luminosity is 1 .5 X 10 35 photons s** 1 . 

Under these assumptions we can predict the Lq, the gamma ray luminosity 
of the galaxy. Writing the pulsar luminosity as a function of time as 

Ly - At- a , (I V-6) 

where A and a are constants appropriate for assumption 3. The pulsar gen- 
eration rate per unit time dn/dt can be expressed as 



(IV-7) 


The number of pulsars with gamma ray luminosity between L and L + dL is 


= ^ y ( » Al/aB L -a/a+n 
dL dt 7 \ dt / a t 


and the total galactic luminosity is 


r L~max , - . 

Lr - f — -L v dL v =— — t 1-1/a \ 

^ dL 7 ^> dL ^ a(l-l/a) rrmax-Smin )• 

•'L-min 


Using the constants 

A = 1 .6 X 10 41 photon s -1 (at t = 1 yr), 

B = 10“ 2 supernova yr -1 , 

L 7 min = 1.5 X 10 35 photon s’ 1 (t = 5 X 10 6 yrs), 


av-8) 


(IV-9) 
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Lyinax = 2.5 X 10 39 (t = 100 years). 


a = 0.9, 


we obtain 


Lq ^ 5 X 10 40 photons s -1 (>35 MeV), (TV-10) 

which is about 5% of the total measured galactic emission. Within just the 
errors of exponent a this value can change between 6 X 10 39 and 2 X lO 4 ^ 
photons s _1 . The additional uncertainties in the “typical” behavior of Crab 
and Vela pulsars increases the range of L G . 

In conclusion, let us note that pulsars can provide a measurable fraction of the 
total galactic gamma ray emission above 35 MeV. With the advent of more 
sophisticated and sensitive instrumentation, not only will the history of pulsar 
emission or its lifetime be better estimated, but also the details of gamma ray 
pulse profiles from these objects can enlighten our knowledge of the pulsars 
significantly. 

Recalling the spirit of the opening paragraph, describing the “underdeveloped” 
nature of the gamma ray region in terms of its “productivity” in the discovery 
of compact objects, let us remind ourselves that it may still hold very rich un- 
exploited reserves. 

Note added in Proof: A more refined analysis of 75 pulsars observed by SAS-2 
has been compiled (Ogelman et. ah, 1976) wherein PSR 1747-46 and PSR 1818- 
04 show positive effects at around 4a level. The rest of the conclusions of 
these lectures do not change significantly. 
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CHAPTER VI 

TEMPORAL ANALYSIS OF X-RAY DATA 

Richard Rothschild 

Goddard Space Flight Center 
Greenbelt, Maryland 

This chapter will be broken into three main categories. The first is “What do 
the data look like?” We shall look at the raw data from a GSFC rocket flight 
in October of 1973; it dramatically shows what we are up against when we 
come to the matter of interpretation. The second section will cover “What 
do we do with the data?” Manipulation of the actual data will be discussed 
and results shown. The final section will cover “What do the results mean?” 
This will include a discussion of how the data bias the results, and also how 
we can interpret the results. 


WHAT DO THE DATA LOOK LIKE? 

Figure VI-1 is a plot of integrated counts per 20.48 ms versus time for the 
exposure to Her X-l , Cyg X-3, Cyg X-l, and, naturally, the diffuse back- 
ground. The plot begins about 1 50 seconds after launch, when the doors of 
the rocket are open and viewing an area of the sky devoid of detectable 
sources. The count rate is about 75 per second. At about 160 seconds after 
launch, the source Her X-l comes into view of the detectors and the rocket 
points at it for fifty-odd seconds. During this time the intensity is seen to 
rise and fall with the 1 .24-second period characteristic of the pulsar. The 
mean rate is 4 to 5 hundred counts per second. The rocket then scanned 
over to Cyg X-3, viewing the diffuse background again and cutting across 
Cyg X-l on the way. When Cyg X-3 entered the field of view of the detectors 
the payload again stopped and pointed at the source for about 35 seconds. 
Cyg X-3 exhibited no discernible structure at this time, and in fact the vari- 
ations in count rate are consistent with those for a random source with con- 
stant mean rate. Next, the rocket began to scan over to Cyg X-l , and when 
the latter was in the field of view, once again the payload held its pointing 
for about a minute. The mean count rate here was about 1300 per second, 
the intensity is quite chaotic, and there seems to be flares of fraction of a 
second duration. Finally the doors closed, exposing a radioactive calibration 
source to the detectors. 

On this flight our finest temporal bin was 320 p%, and exposure to galactic 
x-ray sources totaled some three minutes. Thus we have about a half million 
time bins to work with. What we did with those counts is the subject of the 
next section. 
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WHAT DO WE DO WITH THE DATA? 

rhe analysis of temporal data can be naturally divided into studies of periodic 
and nonperiodic data. Periodic data can be studied by means of power spec- 
trum analysis, which includes the use of the Fast Fourier Transform (FFT) 
algorithm, and light curves. Autocorrelation functions can also be used on 
periodic data, but the discussion of their use will be presented with respect 
to nonperiodic data. Nonperiodic data analysis will include the use of auto- 
correlation functions, plots of variance/integration-time versus integration- 
time, searches for bursting beyond Poisson prediction, computation of burst 
widths, and searches for “ringing” of bursts. 


Power Spectrum 

The power spectrum of a process is the Fourier transform of the autocovari- 
ance function, and represents the distribution of the variance with frequency. 
If a periodic process is being measured, the variance will exceed the overall 
mean variance at a given frequency whenever the frequency is an integer 
times the fundamental frequency. First a few definitions: 

Mx - mean value of the process x(t), 

a^= variance of the process x(t) 

7 xx (r) = autocovariance of the process x(t) 

p xx (r) = autocorrelation of the process x(t) 

They are related as follows. If the process x(t) yields a time series of N counts 
Xj, then 

N 

^ = n J2 *» 

i = 1 
N 

°x ~ n ^ ' ^i “ ^x) 2 * 
i = 1 

N 

7xx( r ) ^ ^ C*i “ ^x) (*i+r * ^x)> 
i = 1 
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PxxOO = >xxW/ a x 
but note that = 7 XX (0). Hence 

^xx( t ) = 7 xx (t)/7 xx (0). 

The power spectrum is defined as 


r„(f)« f 7 xx (r)e-i2^dr, 

J-OO 

which is just the Fourier transform of 7 X x( r )* The i nverse transform states 
7 xx (r)= f r xx (0 e i2nfT df, 

•/_ oo 


If we set r = 0, we see that 




r xx (0df. 


This just confirms what was previously stated. The power spectrum at a given 
frequency r xx (f) df is the contribution to the variance from frequencies f to 
f + df. In order to compare various sets of data, it is necessary to normalize 
the power spectra. This is accomplished by dividing by the variance. This re- 
sult, r xx (0/d x ,is the spectral density function and can be written as 

R xx(0=r xx (f)/a2= f P xx (T)e- i2 * fT dT. 

J _ OO 

An alternate method for generating power spectra is the use of the Fourier 
coefficients of x(t). If x(t) is digitized on N intervals of width T, such that 
NT is the length of the exposure, the Fourier series approximating x(t) is de- 
noted by X'(t) and is given by 


n -1 


x(t) = 


L x - 


m = -n 


e i2?nnt/NT 
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where X m is the complex Fourier coefficient, n = 1/2 N, and the fundamental 
frequency is fj = 1/NT. The coefficients are then given by 


n -l 

x r e -i2,rmr/N , 

r = -n 

where x r are the digitized values of x(t). Then Parseval’s Theorem states that 
the average power of the signal Xj. is given by 


n - 1 n - 1 

x ' T= H iXmi2; 

r = n m = -n 

hence 

r X x (0 = lx m 12 

where 

f = mf x = m/NT. 

Consequently, we are dealing with Fourier transforms and their properties 
are well documented in the literature. We need not derive rigorously any of 
those results, but we shall use the graphic representation used by Brigham 
(1974) in order to illustrate some details that are important. 

Since our data are digitized, so necessarily must be anything calculated directly 
from the data, such as the autocovariance function and consequently the power 
spectrum. There is a theorem called the Sampling Theorem that states how the 
underlying continuous power spectrum or autocovariance can be calculated 
from the discrete ones. We will get to this later. The main use of power spec- 
tra or spectral density functions is to determine periodicities present in data, 
and these methods are very powerful in revealing even weak periodic signals 
in the presence of a strong random signal. Hence, one must be quite careful 
that the analysis does not introduce any false contributions to the spectra. 

If the signal were continuous and of infinite duration such that /fS o x(t) dt 
is finite, we would have no problem in calculating power spectra. Such is not 
the case, however. We will start with a continuous function 7(t) and its 
Fourier transform T(f), as illustrated below: 
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We have assumed that T(f) = 0 for | f | > f c ; that is, there is a critical frequency 
f c above which T(f) is zero. This says that r(f) is band-limited in frequency. 
Actually, x-ray data are band-limited in time (x(t) = 0, for 0 > t and t > T 0 ), 
i.e., there are only counts during the exposure. However, the results are 
analogous and are more easily depicted graphically for the frequency limited 
case. To continue, T(f) is the “pure” continuous power spectrum. It has no 
distortions. In order first to digitize 7 (t), we multiply by d(t), which is an 
infinite series of delta functions spaced at intervals T. The Fourier transform 
of d(t) is D(f) and is three delta functions of height 1/T at f = 0 , ± 1/T. This 
is illustrated below: 


d(tl 

1 


ilii 


LULL 



D(f) 


- 1 / T 


1 1 0 1/T 


The product 7 (t) * d(t) is then a series of delta functions whose amplitude is 
equal to the value of 7 (t) at that value of t. The Fourier transform of 7 (t) • 
d(t) is the convolution of T(f) and D(f) which turns out to be T(f) centered 
at each of the three delta functions D(f). The results are illustrated below; 
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As is shown in the T(f) * D(f) graph, the three T(f) functions do not overlap 
and thus affect each other’s shape. The function T(f) * D(f) centered at f = 0 
is identical to T(f). Hence we have not aliased the data through our digitaliza- 
tion process. This was only because it was drawn with 1/T > 2f c . It is easily 
seen that if 1/T < 2f c the central peak of the T(f) * D(f) plot would no longer 
be identical to T(f) since the overlapping with the side lobes would alter its 
shape. This is just the effect to be avoided. Hence, when picking the digiti- 
zation period T, we make sure that T < 1 / 2 f c . (With x-ray data which usually 
is time domain band-limited, the optimum procedure is to use the smallest T 
possible, so that the overlap is minimized. It cannot be eliminated: there 
will always be some aliasing due to the high frequency side lobes. Also, when 
f > f c T(0 * D(0 begins to repeat and no new information is generated.) This 
critical frequency, 1/2T, is called the Nyquist frequency. In order to avoid 
aliasing, or to minimize it if unavoidable, we always pick 1/T = 2f where T is 
the finest temporal bin possible. In order to remove the 1/T scale difference 
between T( f) and F(f) * D(f), |f I < f c , we multiply T(f) * D(f) by rectangular 
pulse Q(f) of height T and duration 2fc centered at f = 0. The inverse trans- 
form of this is q(t) and both are illustrated below 



-2T.-T 0 T 2T 



Q|f| 





T 



The result (r(f) * D(f) • Q(f) is then the same as the initial r(f). Hence 
(y(t) • d(t))* q(t) is the same as the original y(t). All of this yields the sampl- 
ing theorem. When stated for the usual case in x-ray astronomy (band-limited 
in time) we get that 


r *x(o 2 TV 



sin[27rT c (f-n/2T c )j 

7r(f-n/2T c ) 


hence we can reconstruct the best estimate of the continuous power spectrum 
r xx (0 fr° m ^e digitized power spectrum elements T xx (n/2T C ). 


Another potential problem arises from the finite length of the exposure. This 
is equivalent to multiplying 7(t) by a rectangle x(t) of length T 0 and height 1 . 
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This is illustrated below: 



The problem is obvious: Due to the infinite extent in the frequency domain 
of X(0, there will be a contribution to T xx (f) at all frequencies. This smear- 
ing is called leakage. It is due to the fact that an abrupt turn-on takes an infi- 
nite amount of power. That is, a rectangular pulse is made of contributions 
from cosines and sines of all frequencies. In order to minimize leakage, vari- 
ous experimenters use other “windows” than the rectangular pulse - ones 
that have less steep slopes. This lessens the effect of the data in the initial and 
final portions of the exposure, but for a long enough exposure it can be quite 
valuable. 


When one considers actually computing all the Fourier coefficients from the 
data, one is immediately presented with a problem of computer time. For 
determining the amplitude of N sinusoids the computer time is proportional 
to N 2 . For even a fraction of our 5 X 10 5 points, this is a large number. In 
the early 1960s an algorithm was developed that reduced this number to N 
Log 2 N. Hence there was a substantial savings in time. For 5 X 10 5 points, 
— 2 19 points) the ratio of (N log 2 N)/N 2 is ~4 X 10' 5 . This is the Cooley- 
Tukey algorithm. Later an algorithm based on (2 N )! experiments was de- 
veloped and called the fast Fourier fold. These are used today for computing 
the power spectra. Generally the mean value of x(t) is subtracted from each 
data point, in order to reduce the leakage from the spike at f = 0 which is 
due to the transform of ii. Then the power spectrum for random data has a 
mean of zero, instead of one. 

Figure VI-2 shows the power spectrum of 20 seconds of data from Her X-l . 
The four peaks are at f = .81062, and at two, three, and four times that. This 
frequency corresponds to a period of P = 1 .234 seconds which is near the 
quoted value of 1 .237 seconds. If we had used 40 seconds of exposure the 



TEMPORAL ANALYSIS OF X-RAY DATA 


133 


110 

1 

— I r 

1 

“T 1 

too 







FIRST 20 SECONOS OF 



EXPOSURE TO HER X-l 

90 

- 



- 

80 





70 

- 



- 

k 60 





CJ? 





50 

- 



- 

40 

- 



- 

30 

- 



- 

20 

- 



- 

10 ■ 

- 





t 


n 1 In 

■0 i « n 


1 lyjl 


.5 

1 1.5 

2.0 



MHZ) 

Figure VI-2. Power spectrum of 20 seconds of data from Her X-l. The 
four peaks are f = .81062, 2f, 3f, and 4f. 
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binning would be a factor of two finer and a better value of the period could 
be found. Note the low background in the spectrum at points away from the 
peaks. This is the case for proper use of the techniques explained above. In 
order to demonstrate the effects due to the finite length of the signal, the 
power spectrum was generated for data at the end of the Her X-l exposure. 

It includes about 14 seconds of Her X-l data, then 2 seconds of data when 
the payload was shifting off the source and about 4 seconds of background. 
Hence, we have included data whose average rate changed drastically during 
the exposure. We then would expect significant smearing of the spectrum. 

This spectrum is shown in Figure VI -3. Note the large low-frequency spike 
and the increased contribution away from the peaks. Notice also that the peaks 
are wider. All of this is due to smearing from the transition in the data, and 
points to the desirability of doing the power spectra away from the ends of the 
exposure where the rates are varying. 

Another example of a power spectrum is shown in Figure VI-4. The 55-second 
exposure to Cyg X-l , a non-periodic source on this time scale, was broken into 
five 10-1/2 second exposures, and the principal low frequency components 
are shown. Along the right edge is a histogram of the combined data. There 
is a lot of low frequency power, the vast majority of which does not repeat 
from one 10-second interval to another. This is due to the steep changes in 
intensity over short time intervals exhibited by this source. The early ob- 
servations of Cyg X-l claimed to have found periodicities, but this could not 
be confirmed on later observations. This is all due to the chaotic variations 
in intensity which can produce quasiperiodic temporal structure that last up 
to a minute or two. We now know this structure can be described by a non- 
periodic phenomenon - shot noise, which we shall discuss later. 


Light Curves 

Once the period has been specified, either from a power spectrum analysis or 
another independent experiment, the data can be folded modulo the period. 
This consists of calculating the phase with respect to the period of each tem- 
poral bin or event, and then binning the data by phase. Figure VI-5 shows 
seven Her X-l light curves binned by phase. The width of the bins is deter- 
mined from statistical considerations. The top plot contains data from the 
entire exposure and has not been corrected for background. The level of the 
background is indicated by the black arrow labeled B. The binning is 1% of 
the pulse period and hence this plot has 12.4 ms resolution. Three distinct 
peaks can be seen: one at phase ^.17, a large one at about .65, and a third 
at about .90. Note that the first and second pulse are nearly 180° out of 
phase. 
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Figure VI-3. Power spectrum of 14 seconds of data from Her X-l , 2 seconds of 
data when the payload was moving off of the source, and 4 seconds of diffuse 
background data. 
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10.486 SECONDS 

TIME 


Figure VI-4. Power spectrum of five 10-1/4-second intervals of data from 
Cyg X-l. The principal frequencies below 10 Hz are shown. A histogram 
of the combined data is along the right edge. (Source: Boldt et al., 1975). 
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The next two light curves show data from the first and last quarters of the 
exposure. The similarity of these two is borne out by x 2 tests: they are the 
equivalent. Hence no time variability in the pulse shape on a time scale of 
tens of seconds. The last four show the light curves for data in various energy 
ranges. Statistics prevent a quantitative statement, but it would seem that the 
peaks narrow with energy and might even shift position slightly. 


Autocorrelation Function 

The autocorrelation can be used to reveal periodicities in ddta, but it is far 
less accurate in determining periods, since a knowledge of the pulse profile 
is necessary to fit the data. In recent years it has become useful in the analy- 
sis of non-periodic data such as Cyg X-l . Even though the intensity of a source 
displays no periodicities, there still can be an underlying time structure, e.g., 
shot noise. The autocorrelation function can reveal such structure quite vividly. 

As noted previously the autocorrelation p(r), where r is the lag or shift param- 
eter, is defined as 


P(r) = 


r (T) 

7(0)’ 


where 7(7) is the autocovariance. 7(7) is calculated from the aspect-corrected 
data (corrected to counts seen at the center of the field of view, X i . If there 
are N values of Xj, then 


N — r 

7(r) = ^ -*o 

i - 1 

where r is digitized in units of temporal bins. Thus the autocorrelation func- 
tion can be calculated. This, however, is a composite of two sources of sig- 
nal - the x-ray source itself and the counting statistics. Some of the scatter 
due to statistics can be lessened by averaging p(r) and p(-7), since p(r) is an 
even function and they should be equal. Still, effects due to these two sources 
will still be present. 

Assume the signal z(t) = x(t) + y(t), where x(t) is the source signal and y(t) 
is the counting statistics effect. Then, as described by Weiskopf, et al (1975), 
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PzOO = 7 zz (r)/<^ 
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Thus we have a method for modifying the observed autocorrelation function 
p z (r) to get the x-ray source autocorrelation function p x (r). o z is just 7 ^ 
and is just p. 

As with the light curves, the autocorrelation functions can be performed for 
data in various energy ranges. The functions for various energy ranges can be 
compared in an effort to determine whether effects seen are energy dependent. 
Such an analysis was done on the data from Cyg X-l and Cyg X-3 and the re- 
sults are shown in Figure Vl- 6 . The Cyg X-3 autocorrelation function exhibits 
the characteristics of random noise p( 0) = 1 and p(r) = 0 for r =£ 0 . No struc- 
ture is apparent. On the other hand, the Cyg X-l autocorrelation function dis- 
plays considerable structure. Even though the power spectrum yielded no 
periodicity, there is a defmite exponential character in the data. This says 
that pulses from one time are correlated with pulses from another. We will 
discuss this in the interpretation section later, but suffice it to say that this 
exponential character is the shape of the shot noise pulses making up the Cyg 
X-l intensity. The wavy nature of the autocorrelation function is not due to 
other periodicities, but due to the finite interval being sampled and to the 
shot noise itself. 


Mean Variance 


Another method exists for displaying the data from a source in such a way as 
to allow the identification of the presence of shot noise processes. This em- 
ploys the computation of the mean variance of the data binned over intervals 
of duration AT. Then the plot of < ( 6 N ) 2 >/AT versus AT can be quite reveal- 
ing. Here ((5 N) 2 ) is defined as 


j 

<(5N) 2 ) = }^ 


[Xi-#i(AT)]2 


i = 1 


where J is the number of bins of width AT in the exposure being considered, 
Xj is the number of counts accumulated in bin i, and ju(AT) is the mean value 
of Xj for bins of width AT. Then as AT is varied, J varies, /t(AT) varies, and 
Xj varies, to yield <( 6 N) 2 ) as a function of AT. 

Appendix VI-B contains the derivation of the form that <(SN) 2 ) is expected 
to take as a function of AT for identical rectangular shot noise pulses of dura- 
tion r and frequency of occurrence X. The result is 
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OFFSET TAU (SECONDS) 

Figure VI-6. Autocorrelation p(r) versus lag parameter r for Cyg X-l and 
Cyg X-3 data. The upper plot is from Cyg X-l data and the lower plot is 
from Cyg X-3 data. 
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r r . rat - 

1 

1 rJi + 0 -e) 2 — 
\ 

AT « r ' 

Ir[1 + (1 -e) 2 R/X] 

| 

AT » t J 


where R is the mean count rate per second on source and e is the background 
signal to total signal ratio. Hence, if shot noise exists at the source, for bins 
small with respect to shot pulse width, the ratio MSD should increase linearly 
with AT from the value MSD = R. However, when AT is large with respect to 
the pulse width, MSD is expected to be independent of AT and elevated above 
MSD = R by the amount (1 - e) 2 R 2 /X. 

If shot noise were not present, then r = 0. But to give any count rate at all, 
the rate of such zero-width pulses would have to become infinite in such a 
way as to have Xr be a constant. In this case we only have AT » r , and 
since 1/X -+ 0, 


<(6N) 2 ) 

AT 


= R. 


This is the result to be expected for Poisson-distributed counts: the variance 
equals the mean. Hence the reason for comparing to the MSD = R level. 

Figure VI-7 shows the plot of MSD versus AT for Cyg X-l and Cyg X-3. The 
points for Cyg X-3 scatter about its mean rate with increased scatter as ex- 
pected as increasing AT reduces the number of terms involved in the calcula- 
tion of MSD. On the other hand, the points for Cyg X-l increase with bin 
width and appear to level off for large AT. This is solid evidence for the shot 
noise process occurring at Cyg X-l . Comparing the MSD plot to our theory 
gives values for X and t in agreement (at least for r) with the autocorrelation 
function shape. 


Bursts 

Definition 

Since the discovery of Cyg X-l and its chaotic temporal behavior, searches 
have been made for very short time-scale fluctuations that would be evidence 
for the compact nature of this object. The exposure to Cyg X-l we have been 
using as an example contains several 1-ms bursts. Figure VI-8 shows 80-ms 
of data binned every 640 /is. Quite obvious in this data are at least three 
bursts of 1 ms duration, which implies an emission region of radius less than 
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Figure VI-8. Eighty-two milliseconds of scan temporal data from the 
October 4, 1973 exposure to Cyg X-l , containing 3-millisecond bursts 
(shaded) (Source: Rothschild et al., 1974). 


100 km for a 10M® black hole. Hence, proving that these bursts are statis- 
tically significant is of utmost importance. 

In computing statistical significance, Poisson statistics are employed. One of 
the parameters in the Poisson probability function is the mean rate X. Just 
looking at the data will show that no single X will do (in the Poisson sense 
of X = a 2 ) for describing the entire exposure. Hence the exposure must be 
broken into intervals where the data do not vary radically from the local 
mean. This implies using as short an interval as possible. On the other hand 
the shorter the interval the less the counts and, consequently, the larger the 
statistical uncertainties. A compromise must be reached between the two 
factors. Since our data could conveniently be binned every .4096 second, 
and since the flare that contained the three bursts was of approximately the 
same duration, the exposure was divided into 120 contiguous intervals of 
.4096 second, and then each interval could be divided into 320 bins of 1 .28- 
millisecond bins. This implies that the statistical uncertainty on the predicted 
number of events is about 6%. In each interval the mean number of counts 
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per bin was calculated and the number of counts for each bin was noted. The 
expectation value N(n) for the number of 1 .28 -ms bins containing n counts 
was determined for each interval by using the local mean count per bin and 
Poisson statistics. In each interval a count n* was picked such that the expec- 
tation value for the number of bursts N* = 



n^n* 


was less than or equal to .01, where n* was the lowest possible integer. All 
bins in the interval were then examined to see how many had n > n*. The 
limit .01 was chosen such that when the results from all 120 intervals were 
combined, the expected number of bursts due to statistics would be about 1. 
The result of this analysis, summed over all 120 intervals, was that in six in- 
tervals there existed at least one burst for which n > n*. These six intervals 
contained eight bursts whereas the expectation value for the entire exposure 
is .73. Hence we have shown that on eight occasions throughout the exposure, 
including the three seen in Figure VI-8, bursts occurred when Poisson statistics 
would predict less than 1 . This establishes the fact that Cyg X-l occasionally 
bursts for a millisecond or thereabouts, and incidentally the luminosity can 
increase by a factor of from 4 to 7 over the mean luminosity at these times. 

In order to search for bursts of other durations the bin width was varied, and 
the result was that there were no significant new bursts of shorter or longer 
duration up to about 5 ms. Similarly the duration of the interval over which 
the local mean rate was determined was varied from .1 to .4 second without 
a change in the number of bursts. Hence the interval duration did not affect 
the determination of bursts in the range mentioned above. One thing further 
should be pointed out. Since the probability that a large number of counts 
(a burst) occurs is directly dependent on the local mean rate, the sensitivity 
for revealing bursts drops with increased local mean rate. Thus this method 
of determining bursts is biased against finding bursts in flares or times of en- 
hanced activity as opposed to lower-intensity periods. This makes the ob- 
servation of the three bursts during the flare at =*318 seconds after launch 
even more significant. 


Width 


The next question to be resolved is the pulse width. Our binning in the burst 
analysis was restricted by the telemetry format. Can we get the true burst 
width from the data? If this were a periodic phenomenon there would be 
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no problem - we would just make a light curve phased to the period, as 
for Her X-l , and the mean pulse shape would appear. But here there is no 
periodicity, nothing to synchronize one burst with the next. Since a given 
burst contains about 10 counts in four bins, the statistics prevents us from 
getting a reliable pulse shape by synchronizing on, say, the largest bin. Hence, 
we have chosen to synchronize on the centroid of each burst. The eight bursts 
were binned as finely as possible (320 jts) with the centroids aligned. This is 
shown in Figure VI-9. 



Figure Vl-9. Centroid-aligned burst profile histogram. Error 
bars show statistical error only. 
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The full width at half maximum, corrected for background, is - to the near- 
est channel — three channels, or .96 ms. With the small number of counts com- 
prising the graph, not much can be said for structure within the bursts or really 
about pulse shape. We can say that it is consistent with a rectangular pulse, but 
cannot rule out absolutely a triangular or exponential shape. What is needed is 
finer temporal resolution and more bursts and a larger area of detectors. This 
was partially accomplished on the next rocket flight a year later, when the 
resolution was 160 /as so that the bursts would cover six or seven bins, even 
though the counts per bin would be halved. The bursts were seen again, but 
in lesser numbers than expected on the basis of the first flight. Also, no two 
or more occurred near each other. In any case the centroid-aligned pulse pro- 
file was calculated. Since the error bars are quite large the pulse profile is still 
consistent with a rectangular pulse of 1-ms duration. Even with the finer 
binning, the width, full width at half maximum, is — to the nearest channel — 
1.12 ms. Combining the results of the two flights on 320 gs binning gives 
again ^1 ms as the best guess width of the rectangular pulses. 

Obviously a long exposure with a detector of large area is needed to get a good 
pulse profile (and also to get a good spectrum). The main concern here is 
counts above background. As opposed to radio pulsar studies, for which long 
times are needed, here we need area, i.e., a large number of counts per burst 
with respect to background, not a large number of bursts (even though that 
would not be unwelcome). 


Ringing 

The next question to ask of the data is, “Do the bursts come in bunches or 
alone?” Except for the three that occurred within 20 ms in the first flight, 
all are singles if we just inspect the data by eye. But what if there is an en- 
hancement at some interval before or after each burst that is of the same or- 
der of magnitude as the background? 

This effect will be enhanced, if it exists, by combining all of the bursts as we 
did to get the pulse profile. We can then display this plot for hundreds of 
milliseconds (or any arbitrary time) before and after the bursts to look for 
structure. If the limits are ±T seconds about the burst and two bursts are 
separated by less than T seconds, then the second of the two should be in- 
cluded in the background and not included in the composite burst. This is 
to insure that all bins of the background come frm the same number of bursts. 
This was done for both flights, and the result for the first flight is displayed 
in Figure VI-10, which shows plus and minus 16 ms centered on the aligned 
bursts. The large scatter due to statistics can lead the eye to see a variety of 
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things, but a x 2 test shows the background on either side of the bursts to be 
consistent with a constant. 

Since the scatter in the points goes as \/N and the number of counts in the 
peak goes as N, it is advantageous to gather more bursts. When the bursts of 
the second flight are added in, the statistics only change by about a factor of 
y/T. Hence, the scatter is still large. 


WHAT DO THE RESULTS MEAN? 

The answer to this question could be as detailed and theoretical as desired. 

The intent of this section is to try to describe generally the situation at these 
three sources, and how various types of analysis lead us to those conclusions. 

For Her X-l the situation is basically clear, but muddied by details. The power 
spectrum gave us the period of the pulsar. The light curve gave us the mean 
pulse shape. From the earlier chapters we have learned that pulsars have a beam 
of sorts that sweeps past us to create the pulse. Since two of the bumps are 
180° apart in the pulse profile, the “light house” must have two beams shining 
in opposite directions. Since there are generally two magnetic poles at a mag- 
netized object, the fact that the light curve is double peaked is not surprising. 
The third bump indicates either another beam or some structure to one of the 
two polar beams. 

If we could watch Her X-l for days on end, we would note a 1.7 day period 
present in both the rate data and in the measured pulsar period. This is just 
the binary motion of the stars. We would also notice that this source had 
another periodicity or two, and then the situation becomes complicated. The 
pulse shape changes with binary phase and the emission is only present for 9 
out of 35 days. Hence this “well understood source is not so well understood 
after all. 

From all our tests Cyg X-3 is a random source of x-rays. -Its autocorrelation 
function is zero for r =£ 0. It contains no bursts that cannot be accounted for 
by the Poisson probability distribution. Its power spectrum is featureless on 
our time scale. Longer observations show it to exhibit a 4.8 hour sinusoidal 
variation. But this apparent simplicity is deceptive. Its spectrum varies dra- 
matically, its mean intensity seems to have various levels. And intense radio 
flares, so far uncorrelated with its x-ray or infrared emission, are seen. 

And then there is Cyg X-l . Our observation of bursts indicates a radius of 
less than 100 km. Optical observations predict a mass in excess of 5 M®. It 
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certainly appears to be a black hole. If so, the x-rays we see are coming 
from the inner portions of the accretion disk. Thus the shot noise character 
of its intensity is telling us something about the processes occurring there. 

In Appendix VI-A, the mean, variance, power spectrum and autocorrelation 
for shot noise are derived as functions of the pulse shape. It is seen that the 
observed autocorrelation function shape can be predicted from the assump- 
tion of exponentially shaped pulses. We can get other parameters of the sys- 
tem also. The fit to the autocorrelation function gives the exponential time 
constant. This, combined with the mean (XAr 0 ) and the variance (XA 2 Tq/ 2) 
give the pulse rate X and pulse height A. These turn out to be about 20 per 
second with fall time about 1/3 - 1/2 second, and with about 65 x-ray counts 
detected per pulse. This is true only for the regime AT « r. 

The process of aligning centroids to give a light curve for bursts yielded no 
finer temporal structure. This is only based on a dozen or so bursts, and 
statistical uncertainties are large. The different metrics (Schwartzschild and 
Kerr) predict minimum burst durations. As yet, no distinction between 
metrics can be made. If and when shorter duration time structure is seen, 
perhaps the parameters of the metric can be found. 



APPENDIX VI A 


Derivation of autocorrelation function, spectral density, and probability den- 
sity functions for shot noise processes. 

This derivation will proceed in three steps. First is the derivation of the auto- 
correlation function for Poisson impulses. Second we will convolve this re- 
sult with the shape of an individual shot pulse to reveal both the autocorrela- 
tion function and the spectral density of shot noise. Finally, the probability 
density function will be derived. 

First, a few definitions. 

M x = E[x(t)] = mean of x(t), 

= variance of x(t), 

f(x) = probability density function of x(t), 

PxxOl - ‘2>j 

or > = autocorrelation function of x(t), 

Pxx( T ) J 

TxxOl > h)] 

or > = autocovariance function of x(t) 

TxxOO J 


r xx (0 ~ power spectrum of x(t) 


R* x (0 = spectral density of x(t) 

We know from probability and statistics (e.g., see Papoulis 1965, or Jenkins 
and Watts 1 969) that 

f(x) dx = P[x < x(t) < x + dx], 

N 

Mx =E[x(t)]=l Xj 

i = 1 
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N 

ct x = n F (Xi ~ Mx)2 

i = 1 

The mean and variance characterize a set of data when all the points, are 
independent of each other. This generally is not the case in nature. Two 
other quantities are needed: the autocovariance and autocorrelation func- 
tions. They are defined as: 

TxxOi , t 2 ) = E {[<h)-n x ) -(x(t 2 ) -Mj} , 

PxxOi - h) = TxxOi . t2V°x( t l) o x02) 

Some temporal features are more easily recognized through the power spectrum. 


r xxM = 



Txx( T ) dT - 


where w = 2 tt/t and 7 xx (r) will be defined later, but is just Txx^l > x 2> for 
t j = t, t 2 = t + r. The above can be inverted, if 7 xx (t) is measurable, to give 


?xx( T ) = 


1 

27T 



r xx (w) e itJT du>. 


Then for stationary processes, /u # 


and 


Pxx( T ) 


Txx( T ) 



Rxx(«) : 


£ 


p xx (r)dT, 


So much for definitions. Much of what follows is quite general and has wide 
applications, but the underlying goal is to describe successfully the temporal 
character of Cyg X-l and other sources that display a shot noise type of count 
rate. Hence we shall consider a process whose intensity profile is a composite 
of overlapping and randomly occurring pulses of defined shape. 
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In order to achieve this result we start with the Poisson process that governs 
the number of shot pulses in any time interval. Consider a group of points tj 
randomly distributed on the interval (0, t) with mean rate of X. The prob- 
ability of finding k points in this interval is just 

r ». ag 

k! 

Form the random variable x(t), which is defined as the number of points in 
the interval (0, t). Then the difference x(t a ) - x(t b ) is also a random variable 
and is the number of counts in (t b , t a ). 

Temporarily, let us denote x(t a ) - x(t b ) by co(t a ,t b ). Then it is easy to show 
that 


E I C0 ( t a> tb» = E [ x (t a ) - x 0b)] = * ’ (t a - t b ) = M w (t a , t b ). 


and similarly 


E[w 2 (t a , t b )] = \2(t a - t b )2 + \(t a - t b ); 

hence 

<£(ta.t b ) = X(t a -t b ). 

Thus, this new random variable is also Poisson distributed. Pick four times 
such that t a > t b > t c > t d . The intervals (t b , t a ) and (t d , t c ) are non- 
overlapping and consequently independent. The autocovariance can then be 
written as 

T'xxOa " t b . t c - td) = E {MV t b ) t b )J • [<o(t c , t d ) - t d ]} 

= E[a)(t a , t b ) • co(t c , t d )] - /i u (t a , t b ) (t c , t d ); 

thus 

Txx^a “ *c “ t d) = E[co(t a , t b ) * co(t c , t d )] ~7 2 (t a - t b ) (t c - t d ), 

Let us now evaluate the expectation value: 

E {[x(t a ) - x(t b )] [x(t c ) - x(t d )] } = \2(t a - t b ) (t c - t d ), 

since for two independent random variables the mean of the product is the 
product of the means. Thus in the non-overlap regime 
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T'xx^a " M)> “ *d) 


which makes sense. If, however, the intervals overlap, i.e., t a > t c > t b > t d , 
then intervals (t b , t a ) and (t d , t c ) overlap and are not independent. Adding 
and subtracting x(t c ) and x(t b ) appropriately and regrouping terms gives 

E{[x(t a )-x(t b )] [x(t c )-x(t d )]} =E[u(t a ,t c )u(t c ,t b )| 

- E [co(t a , t c )co(t b f t d )] + E [w 2 (t c , t b )] 

+ E[co(t c , t b )w(t b ,t d )]. 

These are now nonoverlapping intervals and we can apply the formulas dis- 
cussed above: 


E [«(t a , t b ) co(t c , t d )] = X2(t c - t d ) • (t a - t b ) + X(t c - t b ). 

This is identical to the non-overlap case plus a term dependent upon the 
overlap. Then 


Txx(*a ~ ^b’ ~ *d) ^b)* 

If this is generalized by setting t a = ^ , t b = 0,t c =t2>andt d = 0, we get 


E[x(t)] = Xt, where t = t x or t 2 ; 
Xt-i 


Txx^l*^) 


Xt 


1 


* 1^*2 




Now form a new random variable that eventually will be the derivative of x(t) 
with respect to t: 

y(t) = — [x(t + e) - x(t)] , e>0. 


y(t) is Poisson distributed and represents 1/e times the number of points in 
(t, t + e). Let us now calculate the mean and autocovariance of y(t). 
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E[y(t)] = j E[x(t + e)] -I E[x(t)] 

= I\(t + e )-ixt 

= X. 

In order to calculate y yy (tj ,* t 2 ), we must first determine the amount of 
overlap. 

(1) For t^ > t 2 + e, (non-overlap) 


7yy(*l » ^ 2 ) 


(2) For t 2 + e > tj > t 2 , (overlap) 

1 

T'yyOl - * 2 ) = “T TxxOi + e - tj , t 2 + e - 1 2 ) 


Generalizing this gives 


[X(t 2 + e)-Xt 1 ] 


1 . 1 w 
-X-— xOi -t 2 ). 


7yy(tj , t 2 ) 


0 Iti -t 2 |>e~ 

g ^ “3 X |ti - 1 2 | |tj-t 2 |<e 


Finally, form the random variable z(t) as 
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Performing this limiting process yields 
E[z(t)l = X 


TzzOl* t 2> = X6 ( t l 


since |tj - t 2 1 must also go to zero when e -► 0; hence the introduction of the 
delta function. 


This result can be put in a form more readily utilized. Most autocovariances 
are a function of a shift parameter r. Let t 2 = t and tj = t + r. Then 

7 zz ( r ) “ Tzz(^l > ^2) ” Tzz(*» 1 + r ) 

thus 

7 ZZ 00 = X8(t). 


From this we immediately get the variance a \ by setting t - 0: 


a 


2 

z 


= X 


Hence z(t) is also Poisson distributed. The autocorrelation function is then 


Pzz<7) = *(t). 


The power spectrum becomes 


r zz (w) = x 



6(r) e -iuJT dr = X. 


Now we need to fold in the actual pulse profile. Let h(t) be this pulse shape 
and let s(t) be the cumulative effect of shot pulses that occur randomly. Then 
we say 




*)»/ hQ-tj), 


where t 4 are randomly distributed points with uniform density (in time) X. 
The tj are then described by z(t): 
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z(t)= Jj>(t - t j) 
i 

How do we relate the two? They are related mathematically through a con- 
volution: 


i 


's(t) =1 (t - a) h (a) da. 


Assuming h(t) is a real process and obeys casualty (h(t) = 0 for t < 0), this 
becomes 


r 


s(t) = (t - ai h(a) da. 

A quantity we will have need for later is the Fourier transform of h(t). 

r 

H(icu) = h(t) e~ lWT dt. 


Note that 


1 

= f h(t)dt. 

•/- OO 


First, calculate the mean of s(t), assuming X is a constant in time. 


E[s(t)] = 



E [z(t - a)] h(a) da. 


The expectation can go inside the integral since it is performed over t, not a. 
But expectations are constants in time, or so we assumed when we assumed 
X constant; thus 


E[z(t -o)] =E[z(t)] =n z . 
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and so we have 


E[s(t)] 


-*r 


h(a) da 


= M Z H(0). 


What about the autocovariance? By expanding its definition we see that for 
an arbitrary function f(t) 

T ff ( t) = E [f(t) • f(t + t) - /if(t) - M f(t + t) + M 2 ] 

= E [f(t) • f(t + r)]-fi E [f(t)] + E [f(t + t)] + m 2 

= E[f(t) • f(t + t)] -m 2 ; 

hence multiply s(t) by z(t - r). Thus, since z(t - r) is not a function of a, 


s(t) z(t - r) = z(t - a) z(t - r) hfa) da. 

Jo 

Take the expectation of both sides: 

E [s(t) z(t - r)] = f E [z(t - a) z(t - r)] h(a) da. 


Change variables on the right side: t' = t - a, t - r = t f + (a - r). Then the 
expectation value inside the integral becomes 


E[z(t - a) z(t - t)] = 1e[z( 0 z(t' + (a - r))] ; 


but this yields 


E [z(t - a) z(t - r)] = 7 zz (a - r) + 


and since 7 zz (a - r) = 7 ZZ (^ - a )> we g et 


■i. 


E[s(t) z(t - r)] = [7 zz (r - a) + h(a) da. 
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But the left hand side is the cross covariance between s(t) and z(t). Then 
E [s(t) z(t - t )] = y sz (r) + M s M z , 

so we get 

7 sz (r ) = . 7 ZZ (r - a) h(a) da + \i\ H(0) - M s M z • 

•'0 

Or, using the symbol * to denote convolution, 

TszOO = 7 zz (j) * h(T) + nl H(0) -M s M z - 

But we know that 

Ms =M z H(0); 

thus the last two terms cancel to give 

7 s z( r ) = 7 ZZ ( T ) * K t )- 

Similarly, we can show that 

7 ss ( t ) = T sz ( r ) * h(-r). 

Combining these two results yields the desired relation between the autoco- 
variance of s(t) and that for z(t), which has already been calculated. We have 

T'ssO) = 7 zz ( T ) * h(-r) * h(r). 

Now let us proceed to the power spectrum. By definition 


r zz (w) - 



e“ 1WT 7 Z z(r)dr. 
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Similarly 


r zs M ■ 


f- 


- ,wr 7 zs (r)dr 


£ 


7 ZZ (T - a) h(a) da 


dr 


r°° 

r°° 

h( a ) 

I e' iwT 7 ZZ (r - a) dr 

/ -OO 

_ J -.OO 


da 


For the integral in brackets, let 0 = r - a; then 


zsv 


r°° 

r°° 

h(a) 

e-iwW+a) 7zz (j3) <j 0 

/_oo 

•/ -OO 


da 


r 00 

r°° 

I h(a)e -iu,a 

e~^y zz (p)d0 

J - OO 

OO 


da 




h(a) e~ lu,a r zz (to) da 


= r zz (w) H(iw). 

Similarly, we can show that 

r ss (w) = r zs (co)H*(ito). 

Combining these results gives 

r ss(w) = r zz (o;)|H(iro)i2. 
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We now have all we need to compute the various functions associated with 

<t). 


M s = XH(0), 

r ss (co) = X|H(ico)|2, 


1 f 

7ss(r } ~ 2n 

•J —O 


X|H(i«)| 2 e itJT dco. 


This integral is the inverse transform of | H(ico) J 2 which is just h(t) * h (- t), 
since the transform changes multiplication to convolution. Then 


-r 

•'n 


y ss(0 = X h(r + 0)h(/3)d0. 

'o 


This gives us 

o] = x h2(j3) d/?. 

•'O 

The expressions for n s and form Campbell’s theorem. 

f 


P S *(T) 


h(f + P) h(/3) dp 

f h2(p)d0 

^0 


Now we need just supply the proper form of h(t) that describes the shot noise 
process. 


Example: rectangular pulse h(t) 


M 0 _ 


T o 

M 0 
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M 0 /V 0<t<r 0 


0 otherwise 


JU S = XM 0 , 




M 0 \ |1 -e“ iwT 0 


H(icj) = 


, /MqV 1-coswt 0 
|H(ico)| 2 =2 — 0 


/M n \ 2 fl-coswT 0 

r (w) = 2 X — - 

ssl ’ Uo ) CO 2 


We must determine the proper integration limits for calculating 7 SS (T)- For 
h(/3) to be non zero, 0 < 0 < t 0 . For h(r + 0) to be non zero, -r < 0 < r 0 - r. 
Hence the integral is nonzero for 0 < (3 < r 0 - r. Hence 


”l , ss(’ r ) = 


M— J (T °’ T) 0<T<T ° 

0 otherwise 


and finally we get 


1 - r/r n 0 < r < r ( 


0 otherwise 


Example: exponential pulse 


Ae" t/r 0 0 < t 


0 t<0 
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M s = XAr 0 , 

°l = J*A 2 t 0 , 

T 0 

r ss (co) = XA 2 — - 2 — 

1 + “ 2 T l 

7 ss (r) = j\A 2 r 0 e-'/T 0 , 

P„(r) = e-^0, 

which is just the shape of the autocorrelation plot shown for Cyg X-l. 

Can we determine the probability density function and its Fourier transform, 
the characteristic function? First of all, by definition of stationary processes, 

P{s < s(t)< s + ds| = f(s) ds, 

and this is time inedpendent. Let us assume the process h(t) only lasts for a 
time T. Then 

h(t) = 0 for t < 0 and for t > T, 

In the exponential of the second example, we could let T = 10 Tq or 100 r 0 
and be fine. Then let 

x 1 (t) = h(t-t 1 ) > 

where t 1 is a random variable uniformly distributed on (t - T, t). (Remember 
we are looking for effects of things that have already happened; hence t - T.) 
The probability density of X l (t) is g, (x), where 

gj(x) dx = p{x<Xj(t) < x + dx} . 

Now consider 


x 2 (t) = h(t-t 1 )+h(t-t 2 ). 
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The probability density of X 2 (t) is g 2 (x), which is the convolution of gj (x) 
with itself, i.e. 

g2( x ) = gl( x )*8l( x >- 

Extending this to k terms. 


and 



gk( x ) = 8k- 1 ( x ) * 8i ( x ) 


= g 1 (x) * • • • *gj(x) k terms. 

Just for the sake of completeness, define 

g 0 (x) = 5(x) 

Thus if there are k shot pulses in the interval (t - T, t), g k (s) is the probability 
density of s(t). This can be written as 

f[s(t)|n(t-T,t) = k]=g k (s) 

where n(a, b) is a random variable equal to the number of points in (a, b). 
Since everything is randomly distributed, 

s(t) = X k (t), 

Then the probability density function for shot noise, f(s), is given by 
f(s) = f [s| n = 0] P {n = 0 } + f [s| n = 1 ] P{n = l} + ••• + f [s|n = k] P|n = k| 
where P|n = k| is the Poisson probability, i.e. 


So we get 
f(s) = e -xT 


go(s) + gj(s) ■ 


62 w 
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All we need now are the gj(s). If we change variables such that rj = t - tj , 
then Xj(t) = h(rj) where 17 is a random variable on (0, T). Thus we know the 
form of Xj(t). We also know 

gj (x) dx = P jx < Xj (t) < x + dx } , 

and since tj is uniformly distributed over (0, T), the probability of Xj(t) being 
in a band of width dx is just dt/T where dt is the interval corresponding to dx. 
Hence we can say 


gj (x) dx = dt/T 
or 


gl ( X > dT dt = dt / T ’ 

which implies that 

( ^ = I._L- 1 1 

gl W T ' dx/dt T ' dh((J) ' 

d0 

It can also be shown that (Papoules 1965, p. 566) the characteristic function 
< l > (w) for shot noise is 


<h(co) 


Jo 


e iwh(a)_! da 


Then we can get f(s) from the inversion 


f(s) 


■*£ 


4>(w) e- iws dw. 


The integrals start to get troublesome about here, so we shall end. 
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APPENDIX VI-B 

DERIVATION OF VARIANCE AS A FUNCTION 
OF BIN WIDTH 


In order to derive the variance as a function of bin width the following termi- 
nology will be used : 

M 0 = mean counts per burst of duration r 
AT = temporal width of data bin 
X = mean shot burst rate 
N A = total counts in data bin i 
My = counts from burst j in data bin i 
- number of bursts in bin i 
Bj = background counts in bin i 
R = mean on-source count rate 
R b = mean background count rate. 

First of all, 


R = <N>/AT, R b =<B>/AT, e = R„/R. 
Then summing all the contributions to data bin i, 


N i=2^ M ij +B i- 


(Bl) 


j= 1 


The variance is then given by 



(B2) 
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Let us evaluate the partial derivatives: 

9Nj 9Nj 9Nj 

9m^ = 1, ai[ = 1, 9^ =Mi 

Let us show this to be true. The first is obviously true, since the coefficient 
of all My in the sum is one. Hence, the partial derivative with respect to it is 
just one. The second partial is one trivially. The third is not so obvious. How- 
ever, we can treat the partial derivative of a sum similarly to the partial of an 
integral between fixed limits, i.e. 


3Ni 

dn i 


_d__ 

dn i 



Mj(y) dy, 


where Mj(y) is the contribution to the i th data bin from the y th burst. This 
yields 


dNj 

^=M i (n i )-M i (0) = M i , (B3) 

Since Mj(0) = 0. The interpretation of M| is that this is the mean counts per 
burst for the i th data bin. This is not to be confused with the overall mean 
counts per burst <M>. We immediately can write 

Nj = Mj + Bj. 

The variances of My and Bj are given by Poisson statistics as 
(6M ij )2 = M ij and(6B i )=B i . 

We then can write the variance as 

n i 

(5Nj)2 = J2 My + Bj + Mj ( 61^)2 
j = 1 

= nj Mj + Bj + Mj (Snj) 2 


= Nj + Mj (6rij)2. 
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If we average this over the entire exposure, we get 
((5N) 2 > = <N> + <M j (6n i ) 2 >. 

If we assume the burst amplitude to be independent of burst frequency, we 
get 


<(6N) 2 > = <N> + (M(6n) 2 >. 

Here we must enter some assumptions as to the process involved. We choose 
the interpretation that the overall count rate is composed of randomly occur- 
ring, overlapping pulses whose cumulative effect exceeds that of the background 
and other processes that generate structureless counts. This can be quantified 
as follows: 


(a) <(Sn) 2 > = (n>, 

(b) <(5M) 2 >«<M> 2 , 

(c) Mj = <M>+ 5M. 

The first of these just states that the bursts are Poisson distributed; i.e., shot 
noise-like. The second requires that the shot noise process is the dominate 
process in generating the data. The last statement just shows how Mj varies 
from the mean <M>. Incorporating this into Equation B3 yields 

<(6n) 2 > = <N>+<((M> + SM) 2 (6n) 2 > 

= < N > + < (< M > 2 (5 n) 2 + 2< M > 5M(Sn) 2 + (5M) 2 (5 n) 2 

= (N> + <M> 2 (6n) 2 

^<N>+ <M> 2 <n>. (B4) 

Let us now further assume rectangular pulses of width r and total counts M 0 . 
There are two temporal regions in which we can make definite statements 
AT « t and AT » t. 


AT « r : 

Here the data bin width is small compared to the burst width. Then assum- 
ing stationary mean values, we have 
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<M> = M 0 — 

<n> = Ar 

Plugging this into Equation (4) gives us 
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(B5) 

(B6) 


<(6n) 2 > = <N> + M| 




•(Xr) 


since 


M i=- 


Nj ■ Bj 


= RAT + 


and MqAT : 


(R-R b ) 2 AT /at\ 2 


(At ) 2 
RAT - R b AT 




'(Xr), 


Thus we get, after some algebra, 


or 


<(5N) 2 > = RAT 


<(6N) 2 > 

AT 




(B7) 


The other temporal region is for the data bin width large compared to the 
pulse width. Then 


A » T 


<M) = M 0 
< n > = A AT 


(B8) 

(B9) 


This yields, after similar manipulation of equations, 


<(SN ) 2 > 
AT 


= R 


1 + (1 -<0 2 y- 


(BIO) 


The obvious features of these two regions is that for data bins small with re- 
spect to pulse width, the mean square deviation over bin width increases with 
bin width from the Poisson prediction of <N>. For data bins large with respect 
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to pulse width, the mean square deviation over bin width is constant but ele- 
vated above the Poisson prediction by (1 - e) 2 r 2/\. This gives us a direct 
means for measuring the shot noise parameters of an x-ray source. If the data 
are not continuous or cannot readily be put in the form of pure counts, we 
can use intensities as long as we use the proper variance. Let I be the intensity, 
analogous to R above, and let (var) s be the expected statistical variance. Then 
instead of considering mean square deviation over width, we can use observed 
variance over expected variance from Poisson statistics. This gives the added 
advantage that truly random sources will yield unity for this ratio. Then 


Ii = Nj/AT, 


((5N) 2 ) ((SI) 2 ) 

<N> [(var) s I] * 


and for the two regimes 


where 


((51)2) 

(var) s I 


f (I)AT 
1 + — n ; — 





AT «t 


AT»r 


X' = 


X 

(1 -e) 2 
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CHAPTER VII 

RECOGNITION OF COMPACT ASTROPHYSICAL OBJECTS 
THE BINARY PULSAR PSR 1913 + 16 

Steven Detweiler 

University of Maryland 
College Park , Maryland 

In the summer of 1974 Hulse and Taylor (1975) from the University of Mass- 
achusetts made the important discovery of a radio pulsar (PSR 1913+16) in a 
binary system. Hopefully, in the near future this system will allow unprece- 
dented tests of general relativity and will afford new insight into the evolution 
of close binary systems. 


THE PARAMETERS OF PSR 1913+16 

The important parameters of the pulsar and its binary system are given in 
Table VII- 1, and a schematic diagram of the velocity curve is given in Fig- 
ure VII- 1. 


Table VII-1 

Parameters of PSR 1913+16 


p 

Period of Pulsar 

.05902995 S 

IP 1 

Absolute value of time 
derivative of P 

CIO" 16 

DM 

Dispersion measure 

167 ± 5 cm” ^ pc 


Variation with phase 

<10 cm - ^ p C 

P 

Period of binary 

27906.981 ± .002 S 

€ 

Eccentricity 

0.61723 ± .00003 S 

CO 

Longitude of periastron 

179° 

CO 

Precession 

4.24 ± .04 °/year 

a sin i 

Projected semi-major axis 

7.004 X 10 5 km- 1R^ 

(M c sin i) 3 



f ‘(M p + M c )2 

Mass function 

.1313 M® 

- 

No eclipse observed 

- 
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Figure VII-1. The velocity of the binary pulsar as a function of the 
orbital phase (Source: Hulse and Taylor, 1975). 

Two bits of qualitative information are immediately apparent from the velocity 
curve: (1) The eccentricity of the orbit must be large, as is implied by the sharp 
dip in the velocity curve; the pulsar is moving toward us at a high speed for a 
short amount of time. (2) The longitude of the periastron, must be close 
to 180°; this may be inferred from the high degree of symmetry about the 
dip in the velocity curve. 

The period of the pulsar, P, is only .06 s, making it faster^than all other known 
pulsars but the Crab. The time derivative of the period, P, has not yet been 
determined but the upper limit of 10" 16 is usually small. 

Two facts that will be useful in the discussion of the nature of the companion 
of PSR 1913+16 are that no eclipses are observed and that there is no signifi- 
cant correlation between variations in the dispersion measure and the orbital 
phase. 

With an orbital period of only 7.75 hours and a projected semi-major axis of 
the orbit (a sin i), of one solar radius, this binary system will clearly exhibit 
measurable general relativistic effects. In addition, the evidence implies that 
the companion to the pulsar is also a collapsed object ; thus, the effects of the 
tidal interaction are likely to be at least an order of magnitude smaller than 
the general relativistic effects. 
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THE NATURE OF THE COMPANION 1 

We temporarily assume that the companion is a main sequence star; this 
section then shows that such an assumption is untenable. 

It is obvious that the radius of the companion, R c , must be less than the sepa- 
ration of the two components at periastron, r $ . Thus 

( M r sin i\ 

R c <r s = (asini)(l - e) ^ c f -j , (1) 

when f is the observed mass function of the system (cf. Aiken 1964, for a de- 
scription of the classical problems concerning the determination of the param- 
eters of a binary system). In addition, no eclipse of the system is observed; 
thus, R c must be less than the projected stellar separation, r e , when the pulsar 
is at superior conjunction, i.e. 

Re < r e = [(1 + e)/(l + e sin co)] r s cos i. (2) 

With the assumption that the companion is a main sequence star the two in- 
equalities on Rg may be converted to inequalities on M c by using a known 
mass-radius relationship for main sequence stars, 

log (Re/Ro) = log A + B log (M c /M®). (3) 

The constants A and B are determined by the specific model chosen for the 
main sequence stars. Now Inequalities 1 and 2, treated as inequalities on M c , 
may be used with the measured values of a sin i and f to put limits on the four 
quantities M p , M c , a, and sin i. (For details see Masters and Roberts, 1975). 

It is found that for all acceptable values of these parameters there would be 
a classical advance of the periastron of at least 3°/orbit. This shift of the 
periastron would be due to the quadrupole moment of the companion in- 
duced by the tidal interaction with the pulsar. The observed value is about 
three orders of magnitude smaller than this. Table VII-2 gives the relative 
magnitudes for different causes of the apsidal advance. Thus the advance of 
the periastron gives us strong evidence that the companion is not on the main 
sequence. 

Inequality 1 could even be made stronger. The lack of x-rays and the lack of 
correlation between dispersion measure and orbital phase shows that there is 
no evidence of mass transfer or accretion; thus the companion must be inside 
its Roche lobe. 


1 This section is taken mainly from Masters and Roberts (1975). 
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Table VII-2 


Order of Magnitude of Shift in the Periastron 


Quadrupole moment of main sequence 
companion 

^>3°/orbit 

General relativistic effect 

~2°/year 

Quadrupole moment of fast rotating 
white dwarf companion 

^2°/year 

Quadrupole moment of neutron star 
companion 

<1q- 3 °/year 

Quadrupole moment of Helium 
Main sequence companion 

~2 °/year 


There does exist one difficulty with the above conclusion: If the companion 
is a helium burning main sequence star (large centeal density and small radius) 
it is possible that the classical apsidal advance may be as small as the measured 
4 °/ year. While such a possibility cannot yet be ruled out on observational 
evidence, the experts 2 consider it unlikely on evolutionary grounds. 

THE DETERMINATION OF THE MASSES OF THE BINARY SYSTEM 
The Advance of the Periastron 

The measured values of a sin i and f yield two equations for the four unknown 
parameters a, sin i, M p and M c . In order to determine the two masses explicitly 
it is necessary to find at least two more equations governing the four unknowns. 
For these purposes, we make the likely assumptions that the companion is a 
compact object and that general relativity is the correct theory of gravity in 
this regime of physics. 

The general relativistic shift in the periastron is given by 

. GM C 

5w = — — per orbit. (4) 

(1 - e 2 ) ac 2 


2 Informal discussion at the 1975 Enrico Fermi Summer School 
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If the measured 4 °/year is caused solely by the general relativistic effect then 
Equation 4 yields a third relationship between the four unknowns. If a and 
sin i are eliminated from the system it follows that the total mass of the sys- 
tem is 

Mp + M c = 2.87 M@ . (5) 

If a and M p are eliminated, then 

M c sin i = 1 .03 M@ . (6) 

The possible relationships between the masses and sin i are given in Figure 
VII -2. 



Figure VII-2. The masses of components of the binary system as a 
function of i (Source: notes from the Enrico Fermi Summer School, 

1975). 

The assumption that the entire shift in the periastron is due to the general 
relativistic effect does not necessarily follow from the assumption that the 
companion is compact. If it is a rapidly rotating white dwarf (one revolution 
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per 100 seconds or faster) then its rotationally induced quadrupole moment 
may be large enough to give an apsidal advance comparable to the general rela- 
tivistic advance. We will assume that this is not the case — though it is a dis- 
tinct possibility. 


Second Order Doppler Shift and the Gravitational Redshift 

The ratio of the observed period of the pulsar, P ob , to the emitted period, P, 
is given by 

p u v i v2 GM r 

_2^ = i + — + i+ — + + higher order terms, (7) 

P c 2 c rc 2 

where r is the separation of the stars, V is the orbital velocity of the pulsar, 
and ri is a unit vector pointing from the observer to the emitter. 

The term proportional toV/c is just the first order Doppler shift and is suffi- 
cient for determining the mass function of the system. The last two terms 
are the second order Doppler shift and the gravitational redshift which are on 
the order of 10 -6 . 

When the equations for elliptic motion (Aiken, 1964; Blandford and Teukolsky, 
1975) are assumed, Equation 7 takes the form 

p 

—■ = 1 + K[cos (to + 0) + e cos to] + 0 cos 0 + y, (8) 

where 0 is the true anomaly and the constants K, p and y are given by 


and 


2tt a sin i 
cpO-e*) 1 ' 2 ’ 

GM 2 (M p + 2M c )e 
c 2 (M p + M c ) 2 a(l - e 2 ) 


e ’ 2c 2 (Mp + M c ) 2 a ' 


( 9 ) 


( 10 ) 


GO 
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K comes from the first order effect and y cannot be determined without know- 
ing the true period of the pulsar. Thus the next observational goal will be to 
determine the constant j3, yielding a fourth equation for the four unknowns 
and thus finally determining the masses of the system. 

Unfortunately present estimates (Blandford and Teukolsky, 1975) are that it 
will take about five more years of observations to determine 0. 


EFFECTS OF GENERAL RELATIVITY 
Gravitational Radiation 

Once the masses are determined the effects of the gravitational radiation from 
the binary system may be looked for. 

The losses of energy and angular momentum to gravitational radiation will in- 
duce a small change in the period of the binary given by (Peters, 1964; Wagoner, 
1975) 


where 


1 dp _ 96 g3 M p M c (M p + M c ) 
p dt 5 c 5 a 4 


F(e) = (1 - e 2 )- 7/2 + § e2 + ^ e 4 ) • 


For our system this reduces to 


1 dp 
p dt 


-6.3 X 10' 11 


X(1 + X) 3 
sin 5 i 


where 


X = Mp/M c 


( 12 ) 


(13) 


(14) 


(15) 


or approximately one part in 10 10 per year, corresponding to a net loss in 
orbital energy of about 10 30 erg s _1 . It is estimated that this effect will be 
measurable in 5-10 years (Blandford and Teukolsky, 1975). 


The measurement of 1/p dp/dt yields a fifth equation for the unknowns a, 
sin i, M p and M c and thus yields an important check on the consistency of 
our analysis. If the consistency is upheld, then this system will provide 
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verification of three widely different relativistic effects: (A) precession of the 
periastron, Equation (4); (B) gravitational red shift and second order Doppler 
shift, Equation (8); and (C) gravitational radiation, Equation (12). 

If the system of equations is found to be inconsistent, then either some aspect 
of general relativity may be suspect or the assumption of the nature of the 
companion may be wrong (it may be a rapidly rotating white dwarf). 


Precession of the Spin Axis of the Pulsar 

One final general relativistic effect is the coupling of the spin with the orbital 
angular momentum. This coupling may give a precession of the spin axis of 
the pulsar on the order of 2 °/year (Esposito and Harrison, 1975). This may 
cause the plane of the pulsar beam to precess off of the earth, effectively 
shutting off the observations of the pulsars in a time scale of ten years - 
although it might reappear 200 years later. 
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CHAPTER VIII 
PULSAR THEORY 
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University of Maryland 
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It is generally accepted that pulsars are spinning, magnetized neutron stars 
with surface magnetic fields of approximately 10 12 gauss. The mean pulse 
envelope of NP 0532 at various wavelengths is shown in Figure VIII- 1 . The 
simplest physical model consistent with this interpretation of the observed 
properties of pulsars assumes that a neutron star has a dipole magnetic field 
and is surrounded by a vacuum. However, as will be shown below, this model 
is not physically self-consistent since the vacuum solutions of Maxwell’s equa- 
tions indicate that there is an electric field sufficiently strong to pull charged 
particles from the surface of the neutron star and thereby produce a plasma. 

The electromagnetic field on a charge inside the highly conducting neutron 
star must vanish. This requirement implies 


E + 


(fl X r)X B 

~ - 


= 0 , 


0) 


where f2 is the angular velocity. From Equation 1 it follows that the compo- 
nents of E are 


E =-J^sin 2 0, 

* r 


-\ 


~w 


BrJ2 Brf2 3P 7 (cos0) 

E„ = cos 0 sin 0 = 2V ’ 

° c 3c 


Erf, = 0, 


r ( 2 ) 


- s 

where we have assumed that the magnetic and rotational axes are aligned. 


Next we solve for a time-independent solution outside the neutron star and 
then match boundary conditions at the surface. If the field is time-independent, 
it follows from Maxwell’s equations that 


VX E = 0. 


( 3 ) 
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Figure VIII-1. Mean pulse envelope of the Crab pulsar at 
radio, optical, and x-ray wavelengths (Source: Hewish, 
1970), 
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and therefore 


E = - V<b. (4) 

The assumption that the neutron star is surrounded by a vacuum implies 

V • E = V 2 = 0, (5) 

and therefore we have a standard problem in electrostatics. The boundary 
condition at infinity is that <p vanish as r -*■ oo , This condition implies 


$ - ^ ' Cgr-( c+ D Pj (cos 0). 
8=0 


( 6 ) 


The corresponding boundary condition at the surface of the neutron star is 
that the tangential components of E be continuous. This condition is satis- 
fied in a trivial manner for E 0 . The continuity condition for E g implies 


B s Rf2 3P 2 (cos0) 
3c 30 


i 3$ 
R30 



3Pg(cos 0) 
30 


8=o 


(7) 


where R is the radius of the neutron star. The coefficient C 0 = 0 since the 
neutron star is assumed to have no net charge. From Equation 7 it follows 
that 


C 2 = - 


B s R 5 fi 

3c 


Cg = 0 8 = 0, 1 , 8 > 2. 


( 8 ) 


It follows from equations (6) and (8) that the vacuum exterior electric field is 


9$ R 5 B s 

E. — = (3 cos 2 0 - 1), 


dr 


2cr* 


1 3<b r5 B s 


> 


F = _ i __ - 
r dd " 


cr 4 


cos 0 sin 6 


(9) 


184 


RECOGNITION OF COMPACT ASTROPHYSICAL OBJECTS 


Since V X B = 0 and V • B = 0 outside the neutron star, we can proceed 
as above solving for the external magnetic field. It follows that 


B r = B s p cose > 

1 R 3 

B e = 2 Bs sin ^ ’ 


( 10 ) 


At the neutron star’s surface, the gravitational force in the same direction is 


ZeE • B s f2R^B s £e 
mg • B s me GM 


( 11 ) 


which, for M = 1 M@ and B s = 10 12 gauss is ~ 10 7 for iron nuclei and — 10 12 
for electrons. Therefore, charged particles are likely to be pulled from the 
surface of a neutron star, unless strongly bound to the surface by other forces. 


The above discussion implies that it is plausible to assume that there is a plasma 
surrounding a neutron star and that charged particles generally experience no 
net force in this circumstellar plasma. If this assumption is made, then it fol- 
lows from Equation 1 that 


E • B = 0 (12) 

in the region surrounding a neutron star as well as its interior. Equation 12 
is equivalent to assuming that E = 0 in the corotating frame of reference. In 
order to satisfy Equation 1 the charges must redistribute themselves so that 

V • £x B^ = 4?re(n_ - n + ), 


or 


n_ - n + 


a -b 

27rec ’ 


(13) 


where we have assumed that the positively charged particles have unit charge. 
Close to the light cylinder (R L = r sin 9 = c/r) relativistic effects are important 
and it can be shown that the local charge density is 
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n_ - n + 


> B 

2vec 


1 

ft 2 r 2 sin 2 0 
c 2 


(14) 


The electric field produced by the above distribution of charge density is a 
quadrupole electric field since £2 * B changes sign (see Figure VIII-2). 



Figure VIIJ-2. A neutron star magnetosphere for the case where the 
magnetic and rotation axes are aligned 

It is convenient to divide the magnetic field lines outside the neutron star 
into those that pass through the light cylinder and those that close within it. 
The co-rotating magnetosphere is defined to include those magnetic field 
lines that close within the light cylinder. To estimate the boundary of the 
magnetosphere we assume that the poloidal magnetic field is a dipole field. 
For this case the field lines are described by the equation 

r = const. X sin 2 0, (15) 

and therefore the boundary of the magnetosphere includes magnetic field 
lines such that 
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Magnetic fields with directions <0 S at the surface of the neutron star are 
open field lines that extend out to infinity (i.e., to large distances from the 
neutron star). 

Since the electric field on the surface of the neutron star is known (see Equa- 
tion 2), we can calculate the potential difference between any two points on 
the surface. We find 


B R12 

V B- V A = — (sin 2 0 B -sin 2 0 A ). (17) 

The potential difference between the magnetic pole and 0 S is 



Since the potential varies along the surface of the neutron star, the potential 
of a point on the surface cannot in general be the same as the potential at 
large distances, and consequently charges will move along magnetic lines of 
force as they do along wires. In order to maintain charge neutrality there 
must exist a critical angle 0 S between the magnetic pole 0 S at which the po- 
tential is the same as at large distances. It follows that the current flow out 
along open field lines has one sign in the polar regions (0 < 0 S ) and the oppo- 
site sign in the surrounding sheath (0 C < 0 < 0 S ). The most energetic particles 
may attain energies as high as 

1 R ? B 1? 

Emax^ 6 ^^ 1012 -^’ 0*0 

where is the radius of the neutron star in units of 10 6 cm, B 12 is the sur 
face magnetic field in units of 10^ gauss and P is the period of the pulsar in 
seconds. We note that the limiting energy is not strongly dependent on the 
model details such as the relative orientation of B and 12. 

The total magnetic flux that passes beyond the light cylinder is 

<I> = ntp B s • e r , (20) 

where B s is the surface magnetic field and 7rrp is the area of the polar caps 
from which magnetic lines of force penetrate the light cylinder. It follows 
from Equations 13 and 20 that the maximum current from a polar cap is 
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II 


max 



n ' B s 

— i __ 

2itc 2tt 


( 21 ) 


As the above current flows from negatively charged to positively charged 
regions, it crosses magnetic field lines and consequently a torque is exerted 
on the rotating neutron star. This torque acts to lower the rotational energy 
at the rate 



Equations 15, 16 and 20 allow us to rewrite this expression as 

dE B2J24R6 

df < — 7 (23) 

We note that other modes of energy loss such as magnetic dipole radiation 
may be responsible for the loss of a pulsar’s rotational energy. 

In our above discussion we have assumed that the presence of electric forces 
greater than gravity will cause charged particles to be removed from the sur- 
face of a neutron star. However, the physical properties of matter are modi- 
fied in a very significant manner when the magnetic field is as strong as 10 12 
gauss. In such magnetic fields matter will distribute itself in long molecular 
chains along the magnetic fines of force. The ions, which consist of nuclei 
(probably Fe 56 ) and some core electrons, form a one dimensional lattice 
along the magnetic field while the remaining electrons distribute themselves 
in an outer sheath. These molecular chains are very strongly bound to each 
other as a result of strong electric fields. For the case of Fe 56 the binding 
energy per ion is = 14 keV in a magnetic field of=2X 10 12 gauss. The 
Fermi energy of the sheath electrons is =*750 eV and the lattice spacing 
~ 10 cm - The electric field that is required to pull ions from such a chain 
is 14 ke V/Ze Ax where Z — 26 for iron and Ax is the spacing of the one 
dimensional lattice. It follows that the electric field necessary to unbind ions 
is -1° 12 V cm" 1 , which is generally less than the electric field that is predicted 
to exist on the surface of a neutron star, 

nRB s 10 11 

^surface — ” — p volts cm - * , ( 24 ) 
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where P is the rotation period of the pulsar. We can conclude that unless the 
pulsar surface temperature is relatively high kT surface > 2 keV), ions should 
remain bound to the surface. On the other hand, electrons can be much more 
easily pulled from the surface of a neutron star. This circumstance implies 
that electrical discharges will occur on the surface of a pulsar. Such electrical 
discharges are likely to be the cause of the observed pulsed radiation from 
pulsars. 

Let us consider a neutron star whose magnetic and spin axes are oriented so 
that the induced electric charge is positive over its polar caps. This distribu- 
tion of charge is realized when the magnetic axis is antiparallel to the spin 
axis. In this case any positive charge flowing away from the polar regions 
will not be replaced, since the ions are strongly bound to the surface of the 
neutron star. For this reason, a potential gap will develop at the poles. From 
Poisson’s equation, it follows that the potential differences across the gap is 

AV — — — h 2 , (25) 

c 

where h is the gap thickness. Equation 25 remains valid as long as the gap 
thickness is appreciably less than the radius of the polar cap. The maximum 
potential that can exist across this gap is given by Equation 18. 

If a gap develops, its dimension will very rapidly increase with a speed close 
to that of light. However, since the potential difference across the gap varies 
as h 2 (see Equation 25), the gap will soon be discharged by the formation of 
electron-positron pairs. Let us examine more closely how this electron- 
positron discharge might be produced. Electrons accelerated in the gap 
above the polar caps will produce curvature radiation since they are constrained 
to move along magnetic lines of force. The characteristic energy of photons 
produced by such curvature radiation is 

hw — — , (26) 

where r c is the radius of curvature of the magnetic field. Close to the poles 
of a neutron star r c should be 10^ cm. When the gap potential is such that 
eAV > 10 11 eV, the photons that are generated by means of curvature radi- 
ation will be sufficiently energetic that they will create new electron-positron 
pairs on nearby magnetic lines of force. If the discharge is sufficiently well 
developed, these secondary electron-positron pairs will produce new pairs 
and an electron-positron avalanche will be created. The discharge will eventu- 
ally become self-limiting after a sufficient number of electrons are pulled back 
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onto the positively charged polar caps. The gap voltage that is sufficient to 
produce an electron-positron avalanche is ^10 12 eV. Higher voltage drops 
may arise, since the gap thickness increases at nearly the speed of light. 

A crucial physical parameter in determining the mean free path of a gamma 
ray in a magnetic field is 


hco Bf 
2mc 2 ®crit ’ 


(27) 


where B crit = m 2 c 3 /eh = 4.4 X 10 13 gauss and B x = B sin 0 is the component 
of the magnetic field perpendicular to the motion of the gamma ray. Since 
the initial direction of the gamma rays is nearly parallel to the magnetic field 
(i.e., B x < B/y), the parameter x will be « 1 and the mean free path of the 
gamma ray is (Erber 1966). 


4.4 _h Bcrit /_4 \ 
(e 2 /hc) mc B ± CXp \3x / 


(28) 


It is plausible to assume that the mean free path X is comparable to the gap 
thickness h at the time an electron-positron avalanche is produced. We note 
that since the gamma rays are emitted almost parallel to magnetic lines of 
force, they must travel a distance comparable to the radius of curvature r 
before becomes comparable to B. 


Charged particles will achieve energies of y ~ 10 6 inside the gap. Once be- 
yond the gap they will experience no acceleration (i.e., E • B'= 0). Then, how- 
ever, curvature radiation will produce gamma rays and electron-positron pairs 
with energies y ~ 10 2 - 10 3 (i.e., =*50 - 500 MeV). It is possible that the ob- 
served pulsed gamma ray emission at energies >35 MeV is caused by this proc- 
ess. We note that radiative energy losses from curvature radiation will prevent 
electrons or positrons from attaining energies much greater than y ~ 10 6 . 

Our discussion above has led us to conclude that relativistic beams of electrons 
and positrons should stream out from the polar regions of pulsars. It is plaus- 
ible that the observed emission from pulsars should be produced by such beams 
of particles. Moreover, since collective effects are important in understanding 
laboratory plasmas, it is quite reasonable to suppose that they also play a 
fundamental role in determining the observed properties of pulsars. 

At radio wavelengths, pulsars are observed to have very high radio brightness 
temperatures (i.e., T b 10 25 - 10 30 K). A coherent radiation process is 
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required to attain such high brightness temperatures. The bunching of charged 
particles (electrons and positrons) is a plausible mechanism for producing co- 
herent radio emission. It has been suggested that the two-stream instability is 
responsible for the bunching of charged particles. The two-stream instability 
is excited by a charged beam of particles passing through a stationary plasma. 
In our case the plasma is an electron-positron plasma. In order for the two 
stream instability to cause bunching, the growth rate of the instability must 
be greater than the damping caused by coherent curvature radiation. 


The relativistic dispersion relation for the two-stream instability can be ex- 
pressed as (Montgomery and Tidman, 1964) 


where 






7p ( w “ kU P ) 2 7^(0 ) -kU b ) 2 


= 1 , 


47rn P e 2 47rn b e 2 

(Op = — - — and = — — 

Pm b m 


(29) 


are the plasma frequencies of the electron-positron plasma and the ultra- 
relativistic beam respectively. It can be shown that the maximum linear 
growth rate for the two-stream instability is 

3I/2 /n b \l/3 7 

-<v <*» 

Equation 30 implies that the rate of energy gain per cm 3 caused by the two- 
stream instability is 




#here 6E 2 /87t is the electrostatic energy density in the bunches. 


(31) 


In order to compute the corresponding energy loss caused by coherent curva- 
ture radiation we must determine the wave number and phase velocity of the 
fastest growing wave. The wave number and phase velocity that correspond 
to the wave with the maximum growth rate given in Equation 30 are 


k = 2tt/X = 




Au7^ /2 


(32) 
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and 


u h tl-2 


-4/3 




respectively. The quantity Au is defined as 


Au = u b - v 


(33) 


(34) 


The longitudinal waves excited by the two-stream instability will cause bunch- 
ing of charged particles. Since all of the above expressions have been evaluated 
in the observer’s frame of reference, the bunches have a characteristic dimen- 
sion ^ XI 2 and move with the velocity u 0 in this frame of reference. It follows 
that the rate of energy loss per cm 3 due to coherent curvature radiation from 
a bunch of electrons (or positrons) is 



(35) 


where 



r c is the radius of curvature, and 5 n is the number of electrons (or positrons) 
per cm 3 in a bunch. 


From Poisson’s equation we find 


k5E = 47re6n. ( 36 ) 

It follows from Equations 35 and 36 that the rate of energy loss from curva- 
ture radiation is 


('£') =„!V — 

\<li /loss ,2 * 8» • 
c 


(37) 


Since both the rate of energy loss and energy gain are proportional to the 
electrostatic energy density, we can define a damping rate by the expression 
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^damping “ 2 ~^<!> ' 
r c 

The condition that the two stream instability develop becomes 

^damping- 


( 38 ) 


(39) 


This condition for bunching makes it implausible for a beam of positrons (or 
electrons) as energetic as y ~ 10 6 to excite charge bunches. However, it does 
not rule out the possibility that electrons (and/or positrons) with 7 ~ 10 3 can 
excite bunches as they penetrate a lower energy electron-positron plasma. One 
additional condition necessary for the two-stream instability to be excited is 
the requirement that the electron and/or positron beams be separated in energy. 
If the energy distribution of particles in the beam is monotonically decreasing, 
then no two-stream instability will develop. 

The two-stream instability is only one of a number of plasma instabilities that 
may be excited by charged relativistic particles moving through a plasma. The 
transverse electromagnetic instability (ordinary mode instability) is an example 
of a plasma instability that may have important application to pulsar theory . 
Recent studies of relativistic electron beams have shown that the ordinary mode 
instability can cause the filamentation of electron beams moving through a 
plasma. The filamentation of the electron beam is caused by a perturbed mag- 
netic field that is perpendicular to the streaming motion of the beam. 

It can be shown (e.g., Hardee and Rose, 1974) that the condition for the ex- 
citation of those transverse electromagnetic waves is 


cu 


?b 


b > & 



(40) 


for the case of an electron (or positron) beam moving along a magnetic field 
embedded in an electron-positron plasma. The ordinary mode instability pro- 
duces a purely growing mode whose maximum growth can be shown to be 


co? 

•max 7 b 


(41) 


The maximum growth rate is achieved for beam and background plasma of 
equal density. If we compare the growth rate for the ordinary mode instabil- 
ity with that for the two-stream instability given in Equation 30, we see that 
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the ordinary mode instability, if the necessary condition for its excitation 
(Equation 40) is realized, can develop more rapidly for a relativistic beam 
than can the two-stream instability. 

Close to the surface of a neutron star the magnetic field strength is very 
strong and consequently the plasma should be stable against the excitation 
of the ordinary mode instability. However, since the magnetic field outside 
the surface of a neutron star should scale as B = B s (R/r) 3 and the particle 
density as n = n s (R/r) 3 , particle densities sufficiently high to drive the ordi- 
nary mode instability may arise inside the light cylinder (i.e., R L — 1.5 X 10 8 
cm for the crab pulsar). The excitation of the ordinary mode instability im- 
plies that a particle moving along a magnetic field will see a transverse mag- 
netic field excited by the instability. The presence of a magnetic field trans- 
verse to the direction of motion of the particle (electron or positron) will 
cause it to radiate by means of the synchrotron process. It is plausible that 
the x-ray and optical pulses observed from the crab pulsar are produced in 
this manner. 

If the initial particle distribution is of a two-stream type, then in the regime 
where 


<4h h «V?hl ( 42 ) 

there exists an electrostatic two-stream instability and an associated low- 
frequency transverse wave. The phase propagation of the transverse wave 
occurs at any angle (other than zero) with respect to the ambient magnetic 
field. The group velocity, however, is along the ambient magnetic field. It 
has been suggested (Hardee and Rose, 1976) that the observed radio pulses 
from pulsars may be caused by the focusing of these transverse waves along 
the open field lines extending from the polar cap to the light cylinder. 
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